mirror of
https://asciireactor.com/otho/psdlag-agn.git
synced 2024-11-22 08:35:06 +00:00
updated log, setting up analysis script with proper initial parameters.
This commit is contained in:
parent
0f2af5fa04
commit
2981743681
651
analyses/all
651
analyses/all
@ -1,651 +0,0 @@
|
||||
==> g(𝛌=4775)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
-9.745e-01 2.756e-01
|
||||
-1.384e+00 2.128e-01
|
||||
-2.748e+00 2.454e-01
|
||||
-3.305e+00 2.335e-01
|
||||
-3.314e+00 1.777e-01
|
||||
-3.389e+00 1.613e-01
|
||||
-4.198e+00 1.693e-01
|
||||
-4.465e+00 1.915e-01
|
||||
-4.700e-01 1.237e+00
|
||||
-7.487e-01 1.302e+00
|
||||
-2.046e+00 1.348e+00
|
||||
-2.428e+00 1.335e+00
|
||||
-2.953e+00 1.392e+00
|
||||
-3.086e+00 1.396e+00
|
||||
-3.761e+00 1.443e+00
|
||||
-4.290e+00 1.408e+00
|
||||
9.862e-02 3.790e-01
|
||||
3.899e-01 4.121e-01
|
||||
8.650e-01 4.290e-01
|
||||
5.516e-01 4.096e-01
|
||||
2.228e-01 4.090e-01
|
||||
9.508e-01 4.215e-01
|
||||
-2.872e-01 3.988e-01
|
||||
9.059e-02 3.866e-01
|
||||
***************************************
|
||||
|
||||
==> HST(𝛌=1158)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
2.402e-01 2.838e-01
|
||||
1.329e-01 2.104e-01
|
||||
-8.990e-01 2.184e-01
|
||||
-1.229e+00 1.729e-01
|
||||
-1.825e+00 1.514e-01
|
||||
-2.236e+00 1.569e-01
|
||||
-2.758e+00 1.438e-01
|
||||
-4.749e+00 3.931e+00
|
||||
1.518e-01 3.556e-02
|
||||
1.794e-02 4.241e-03
|
||||
-1.040e+00 2.554e-01
|
||||
-1.377e+00 3.397e-01
|
||||
-1.989e+00 4.898e-01
|
||||
-2.375e+00 5.854e-01
|
||||
-2.886e+00 7.127e-01
|
||||
-3.737e+00 9.211e-01
|
||||
-3.034e-02 7.945e-03
|
||||
-6.855e-02 1.742e-02
|
||||
-6.829e-02 1.598e-02
|
||||
-4.907e-02 1.126e-02
|
||||
-5.681e-02 1.297e-02
|
||||
-1.644e-01 3.816e-02
|
||||
-3.273e-01 8.072e-02
|
||||
1.333e-01 3.447e-02
|
||||
***************************************
|
||||
|
||||
==> HST(𝛌=1367)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
6.534e-02 2.733e-01
|
||||
-9.689e-02 2.184e-01
|
||||
-1.175e+00 2.494e-01
|
||||
-1.525e+00 1.446e-01
|
||||
-2.167e+00 1.617e-01
|
||||
-2.492e+00 1.545e-01
|
||||
-3.258e+00 1.095e-01
|
||||
-1.115e+01 5.296e+03
|
||||
6.378e-02 1.554e-02
|
||||
-9.676e-02 2.400e-02
|
||||
-1.181e+00 2.892e-01
|
||||
-1.520e+00 3.738e-01
|
||||
-2.180e+00 5.369e-01
|
||||
-2.532e+00 6.061e-01
|
||||
-3.129e+00 7.691e-01
|
||||
-3.794e+00 9.654e-01
|
||||
-8.535e-05 2.174e-05
|
||||
-5.126e-04 1.233e-04
|
||||
-4.066e-04 9.595e-05
|
||||
1.080e-03 2.497e-04
|
||||
-2.682e-03 6.077e-04
|
||||
-5.462e-03 1.334e-03
|
||||
-2.209e-02 5.376e-03
|
||||
1.564e-01 3.945e-02
|
||||
***************************************
|
||||
|
||||
==> HST(𝛌=1479)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
1.224e-01 2.666e-01
|
||||
-7.607e-02 2.200e-01
|
||||
-1.174e+00 2.513e-01
|
||||
-1.500e+00 1.446e-01
|
||||
-2.075e+00 1.576e-01
|
||||
-2.554e+00 1.536e-01
|
||||
-3.337e+00 1.052e-01
|
||||
-1.152e+01 9.495e+03
|
||||
9.340e-02 2.202e-02
|
||||
-8.842e-02 2.148e-02
|
||||
-1.173e+00 2.885e-01
|
||||
-1.512e+00 3.787e-01
|
||||
-2.113e+00 5.217e-01
|
||||
-2.547e+00 6.151e-01
|
||||
-3.178e+00 7.800e-01
|
||||
-7.000e+00 1.520e+00
|
||||
-6.035e-02 1.487e-02
|
||||
2.044e-02 5.014e-03
|
||||
1.258e-01 3.039e-02
|
||||
5.764e-02 1.497e-02
|
||||
2.580e-02 6.105e-03
|
||||
-4.737e-02 1.144e-02
|
||||
9.579e-02 2.320e-02
|
||||
6.411e-01 1.538e-01
|
||||
***************************************
|
||||
|
||||
==> HST(𝛌=1746)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
3.934e-02 2.650e-01
|
||||
-1.638e-01 2.223e-01
|
||||
-1.339e+00 2.641e-01
|
||||
-1.655e+00 1.414e-01
|
||||
-2.232e+00 1.613e-01
|
||||
-2.793e+00 1.845e-01
|
||||
-4.025e+00 1.005e-01
|
||||
-1.072e+01 3.151e+03
|
||||
5.043e-02 1.188e-02
|
||||
-1.304e-01 3.177e-02
|
||||
-1.254e+00 3.099e-01
|
||||
-1.587e+00 4.030e-01
|
||||
-2.207e+00 5.484e-01
|
||||
-2.693e+00 6.488e-01
|
||||
-3.209e+00 8.109e-01
|
||||
-3.820e+00 9.624e-01
|
||||
3.040e-02 8.141e-03
|
||||
2.438e-02 6.215e-03
|
||||
1.450e-01 3.324e-02
|
||||
1.458e-01 3.170e-02
|
||||
1.772e-01 3.882e-02
|
||||
2.076e-01 4.966e-02
|
||||
-1.390e-01 3.462e-02
|
||||
3.929e-01 1.097e-01
|
||||
***************************************
|
||||
|
||||
==> i(𝛌=7647)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
-1.510e+00 2.756e-01
|
||||
-1.872e+00 2.087e-01
|
||||
-3.441e+00 2.422e-01
|
||||
-3.722e+00 1.931e-01
|
||||
-4.794e+00 3.059e-01
|
||||
-4.336e+00 1.760e-01
|
||||
-4.644e+00 1.360e-01
|
||||
-5.061e+00 1.844e-01
|
||||
-7.402e-01 1.117e+00
|
||||
-9.967e-01 1.174e+00
|
||||
-2.393e+00 1.214e+00
|
||||
-2.723e+00 1.271e+00
|
||||
-3.701e+00 1.349e+00
|
||||
-3.742e+00 1.343e+00
|
||||
-4.317e+00 1.393e+00
|
||||
-4.979e+00 1.451e+00
|
||||
2.638e-01 3.565e-01
|
||||
7.095e-01 3.780e-01
|
||||
1.297e+00 4.108e-01
|
||||
1.035e+00 3.960e-01
|
||||
9.267e-01 3.966e-01
|
||||
8.982e-01 4.049e-01
|
||||
7.318e-01 3.899e-01
|
||||
9.092e-01 3.982e-01
|
||||
***************************************
|
||||
|
||||
==> I(𝛌=8560)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
-1.497e+00 3.256e-01
|
||||
-1.786e+00 2.808e-01
|
||||
-3.555e+00 5.022e-01
|
||||
-4.203e+00 1.412e-01
|
||||
-3.826e+00 2.505e-01
|
||||
-4.096e+00 5.759e-01
|
||||
-4.304e+00 1.379e-01
|
||||
-1.109e+01 5.792e+03
|
||||
-7.542e-01 1.301e+00
|
||||
-9.532e-01 1.415e+00
|
||||
-2.397e+00 1.461e+00
|
||||
-3.214e+00 1.450e+00
|
||||
-3.008e+00 1.544e+00
|
||||
-3.489e+00 1.514e+00
|
||||
-4.353e+00 1.499e+00
|
||||
-4.061e+00 1.520e+00
|
||||
1.675e-01 4.363e-01
|
||||
7.660e-01 4.491e-01
|
||||
1.781e+00 4.834e-01
|
||||
1.971e+00 4.573e-01
|
||||
6.572e-01 4.531e-01
|
||||
8.520e-01 4.811e-01
|
||||
2.405e+00 4.975e-01
|
||||
-1.573e-01 4.341e-01
|
||||
***************************************
|
||||
|
||||
==> plots <==
|
||||
|
||||
==> r(𝛌=6175)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
-1.372e+00 2.752e-01
|
||||
-1.786e+00 2.125e-01
|
||||
-3.214e+00 2.522e-01
|
||||
-3.558e+00 2.087e-01
|
||||
-4.416e+00 3.027e-01
|
||||
-4.366e+00 2.391e-01
|
||||
-4.559e+00 1.732e-01
|
||||
-4.586e+00 1.687e-01
|
||||
-6.693e-01 4.559e-01
|
||||
-5.998e+00 1.480e+00
|
||||
-2.390e+00 9.034e-01
|
||||
-2.927e+00 1.017e+00
|
||||
-3.411e+00 1.073e+00
|
||||
-3.731e+00 1.204e+00
|
||||
-3.930e+00 1.152e+00
|
||||
-7.000e+00 1.390e+00
|
||||
2.905e-01 1.874e-01
|
||||
1.835e+00 4.578e-01
|
||||
8.760e-01 3.014e-01
|
||||
1.476e+00 3.951e-01
|
||||
-3.280e-01 2.003e-01
|
||||
1.439e+00 4.270e-01
|
||||
9.464e-02 1.784e-01
|
||||
1.532e+00 4.331e-01
|
||||
***************************************
|
||||
|
||||
==> R(𝛌=6439)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
-1.287e+00 2.795e-01
|
||||
-1.580e+00 2.127e-01
|
||||
-2.879e+00 2.577e-01
|
||||
-3.569e+00 2.889e-01
|
||||
-3.661e+00 2.296e-01
|
||||
-4.471e+00 5.134e-01
|
||||
-3.998e+00 1.564e-01
|
||||
-5.001e+00 4.661e-01
|
||||
-6.583e-01 7.106e-01
|
||||
-6.900e+00 1.424e+00
|
||||
-2.060e+00 9.613e-01
|
||||
-2.727e+00 1.047e+00
|
||||
-2.954e+00 1.131e+00
|
||||
-3.654e+00 1.178e+00
|
||||
-4.406e+00 1.288e+00
|
||||
-5.203e+00 1.386e+00
|
||||
1.448e-01 2.645e-01
|
||||
-1.004e-01 2.434e-01
|
||||
9.754e-01 3.634e-01
|
||||
1.140e+00 3.905e-01
|
||||
9.284e-01 3.632e-01
|
||||
1.116e+00 3.672e-01
|
||||
2.092e+00 4.941e-01
|
||||
8.634e-01 3.384e-01
|
||||
***************************************
|
||||
|
||||
==> Swiftb(𝛌=4392)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
-5.521e-01 3.134e-01
|
||||
-1.135e+00 2.562e-01
|
||||
-2.691e+00 4.640e-01
|
||||
-3.626e+00 1.546e-01
|
||||
-4.510e+00 2.215e-01
|
||||
-1.036e+01 5.116e-01
|
||||
-1.076e+01 1.445e-01
|
||||
-1.213e+01 8.139e+04
|
||||
-2.515e-01 3.686e-01
|
||||
-6.156e-01 4.904e-01
|
||||
-1.924e+00 8.403e-01
|
||||
-2.417e+00 8.016e-01
|
||||
-3.053e+00 1.013e+00
|
||||
-6.570e+00 1.397e+00
|
||||
-3.450e+00 1.038e+00
|
||||
-6.333e+00 1.336e+00
|
||||
7.408e-02 1.220e-01
|
||||
3.089e-01 2.006e-01
|
||||
5.057e-01 2.378e-01
|
||||
4.831e-01 1.961e-01
|
||||
5.779e-01 2.557e-01
|
||||
6.277e-01 2.966e-01
|
||||
4.823e-01 2.573e-01
|
||||
6.110e-01 2.531e-01
|
||||
***************************************
|
||||
|
||||
==> Swiftu(𝛌=3465)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
-3.178e-01 3.178e-01
|
||||
-7.537e-01 2.541e-01
|
||||
-2.402e+00 4.048e-01
|
||||
-2.624e+00 1.540e-01
|
||||
-3.413e+00 1.944e-01
|
||||
-3.834e+00 2.633e-01
|
||||
-1.236e+01 1.174e-01
|
||||
-1.071e+01 3.318e+03
|
||||
-1.376e-01 3.796e-01
|
||||
-4.256e-01 4.494e-01
|
||||
-1.794e+00 6.810e-01
|
||||
-2.123e+00 6.992e-01
|
||||
-2.817e+00 9.699e-01
|
||||
-3.029e+00 9.065e-01
|
||||
-4.139e+00 1.273e+00
|
||||
-7.000e+00 1.331e+00
|
||||
7.840e-02 1.316e-01
|
||||
3.388e-01 1.716e-01
|
||||
6.616e-01 2.532e-01
|
||||
1.349e-01 1.501e-01
|
||||
9.909e-01 3.749e-01
|
||||
8.538e-01 3.311e-01
|
||||
1.397e+00 4.674e-01
|
||||
6.412e-01 2.404e-01
|
||||
***************************************
|
||||
|
||||
==> Swiftv(𝛌=5468)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
-8.802e-01 3.314e-01
|
||||
-1.441e+00 2.854e-01
|
||||
-2.883e+00 4.166e-01
|
||||
-4.483e+00 1.547e-01
|
||||
-1.045e+01 2.071e-01
|
||||
-1.041e+01 3.271e-01
|
||||
-1.042e+01 1.274e-01
|
||||
-1.078e+01 5.879e+03
|
||||
-4.161e-01 9.142e-01
|
||||
-7.660e-01 1.130e+00
|
||||
-2.020e+00 1.179e+00
|
||||
-3.201e+00 1.143e+00
|
||||
-2.699e+00 1.335e+00
|
||||
-2.859e+00 1.254e+00
|
||||
-3.653e+00 1.312e+00
|
||||
-5.501e+00 1.434e+00
|
||||
1.900e-02 2.266e-01
|
||||
3.560e-01 3.285e-01
|
||||
1.007e+00 3.907e-01
|
||||
-1.270e+00 3.724e-01
|
||||
1.858e+00 4.560e-01
|
||||
9.302e-01 4.464e-01
|
||||
-1.502e+00 4.728e-01
|
||||
2.653e+00 4.849e-01
|
||||
***************************************
|
||||
|
||||
==> tables <==
|
||||
|
||||
==> u(𝛌=3471)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
-4.289e-01 2.733e-01
|
||||
-9.017e-01 2.118e-01
|
||||
-2.581e+00 3.058e-01
|
||||
-2.751e+00 2.567e-01
|
||||
-2.764e+00 1.891e-01
|
||||
-3.105e+00 2.139e-01
|
||||
-3.298e+00 1.431e-01
|
||||
-3.870e+00 2.222e-01
|
||||
-1.925e-01 3.059e-01
|
||||
-5.093e-01 3.513e-01
|
||||
-1.923e+00 6.007e-01
|
||||
-2.204e+00 6.982e-01
|
||||
-2.568e+00 8.240e-01
|
||||
-2.768e+00 8.080e-01
|
||||
-3.425e+00 9.196e-01
|
||||
-5.982e+00 1.396e+00
|
||||
1.646e-01 1.104e-01
|
||||
4.376e-01 1.640e-01
|
||||
7.202e-01 1.945e-01
|
||||
9.556e-01 2.502e-01
|
||||
4.856e-01 1.911e-01
|
||||
2.985e-01 1.448e-01
|
||||
-6.801e-01 2.362e-01
|
||||
-2.416e+00 4.887e-01
|
||||
***************************************
|
||||
|
||||
==> UVM2(𝛌=2246)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
-1.003e-01 3.133e-01
|
||||
-5.235e-01 2.446e-01
|
||||
-1.863e+00 3.017e-01
|
||||
-2.159e+00 1.943e-01
|
||||
-2.718e+00 2.380e-01
|
||||
-3.504e+00 4.807e-01
|
||||
-1.139e+01 1.494e+04
|
||||
-4.119e+00 8.683e-01
|
||||
-2.431e-02 1.310e+00
|
||||
-3.121e-01 1.446e+00
|
||||
-1.526e+00 1.472e+00
|
||||
-1.877e+00 1.484e+00
|
||||
-2.455e+00 1.493e+00
|
||||
-2.945e+00 1.508e+00
|
||||
-3.356e+00 1.510e+00
|
||||
-3.602e+00 1.508e+00
|
||||
-1.862e-04 4.058e-01
|
||||
1.380e-01 4.538e-01
|
||||
2.475e-01 4.630e-01
|
||||
3.102e-02 4.416e-01
|
||||
6.040e-01 4.608e-01
|
||||
-1.488e-01 4.709e-01
|
||||
1.171e+00 4.637e-01
|
||||
1.617e+00 4.645e-01
|
||||
***************************************
|
||||
|
||||
==> UVW1(𝛌=2600)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
-2.417e-01 3.190e-01
|
||||
-7.027e-01 2.509e-01
|
||||
-1.998e+00 3.947e-01
|
||||
-2.387e+00 1.542e-01
|
||||
-2.912e+00 1.882e-01
|
||||
-4.026e+00 2.472e-01
|
||||
-1.069e+01 1.123e-01
|
||||
-9.716e+00 2.578e+02
|
||||
-9.867e-02 6.982e-01
|
||||
-4.150e-01 8.299e-01
|
||||
-1.587e+00 8.040e-01
|
||||
-1.963e+00 9.829e-01
|
||||
-2.539e+00 1.024e+00
|
||||
-3.274e+00 1.148e+00
|
||||
-3.460e+00 1.229e+00
|
||||
-6.742e+00 1.428e+00
|
||||
6.734e-02 2.545e-01
|
||||
1.353e-01 2.890e-01
|
||||
3.585e-01 2.644e-01
|
||||
1.587e-01 2.261e-01
|
||||
2.959e-01 2.424e-01
|
||||
1.111e-01 2.882e-01
|
||||
2.364e-01 2.837e-01
|
||||
7.463e-01 3.484e-01
|
||||
***************************************
|
||||
|
||||
==> UVW2(𝛌=1928)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
-1.648e-01 2.888e-01
|
||||
-4.974e-01 2.247e-01
|
||||
-1.686e+00 2.872e-01
|
||||
-1.897e+00 1.430e-01
|
||||
-2.600e+00 1.717e-01
|
||||
-2.801e+00 2.023e-01
|
||||
-3.814e+00 1.056e-01
|
||||
-9.888e+00 3.109e+02
|
||||
-6.101e-02 2.539e-01
|
||||
-3.020e-01 3.388e-01
|
||||
-1.433e+00 4.484e-01
|
||||
-1.756e+00 5.656e-01
|
||||
-2.485e+00 7.064e-01
|
||||
-2.799e+00 8.287e-01
|
||||
-3.317e+00 9.293e-01
|
||||
-7.000e+00 1.393e+00
|
||||
4.622e-02 9.341e-02
|
||||
1.473e-01 1.171e-01
|
||||
2.243e-01 1.191e-01
|
||||
1.606e-01 1.081e-01
|
||||
2.365e-01 1.249e-01
|
||||
1.484e-01 1.182e-01
|
||||
8.110e-01 2.388e-01
|
||||
9.689e-01 2.720e-01
|
||||
***************************************
|
||||
|
||||
==> V(𝛌=5404)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
-1.044e+00 2.743e-01
|
||||
-1.457e+00 2.109e-01
|
||||
-2.848e+00 2.322e-01
|
||||
-2.957e+00 1.713e-01
|
||||
-4.232e+00 2.257e-01
|
||||
-4.052e+00 1.742e-01
|
||||
-4.515e+00 1.593e-01
|
||||
-5.250e+00 3.713e-01
|
||||
-5.026e-01 1.143e+00
|
||||
-7.878e-01 1.224e+00
|
||||
-2.118e+00 1.280e+00
|
||||
-2.276e+00 1.263e+00
|
||||
-3.421e+00 1.357e+00
|
||||
-3.378e+00 1.364e+00
|
||||
-4.224e+00 1.433e+00
|
||||
-4.930e+00 1.494e+00
|
||||
9.142e-02 3.750e-01
|
||||
4.814e-01 3.731e-01
|
||||
8.550e-01 3.884e-01
|
||||
3.929e-01 3.841e-01
|
||||
3.975e-01 3.691e-01
|
||||
5.311e-01 3.951e-01
|
||||
-1.203e-01 3.770e-01
|
||||
-2.905e+00 5.243e-01
|
||||
***************************************
|
||||
|
||||
==> z(𝛌=9157)≻HST(𝛌=1367)_Δt=0.1 <==
|
||||
# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
|
||||
6.533e-02 2.816e-01
|
||||
-9.694e-02 2.103e-01
|
||||
-1.175e+00 2.185e-01
|
||||
-1.525e+00 1.729e-01
|
||||
-2.166e+00 1.536e-01
|
||||
-2.492e+00 1.538e-01
|
||||
-3.258e+00 1.803e-01
|
||||
-9.328e+00 7.607e+01
|
||||
-1.852e+00 2.783e-01
|
||||
-2.156e+00 2.086e-01
|
||||
-3.798e+00 2.699e-01
|
||||
-4.146e+00 2.393e-01
|
||||
-4.577e+00 2.541e-01
|
||||
-4.338e+00 1.832e-01
|
||||
-4.672e+00 1.459e-01
|
||||
-5.230e+00 2.311e-01
|
||||
-9.339e-01 1.264e+00
|
||||
-1.131e+00 1.351e+00
|
||||
-2.536e+00 1.393e+00
|
||||
-3.031e+00 1.413e+00
|
||||
-4.060e+00 1.498e+00
|
||||
-3.801e+00 1.463e+00
|
||||
-4.139e+00 1.444e+00
|
||||
-5.151e+00 1.548e+00
|
||||
2.445e-01 3.942e-01
|
||||
8.131e-01 4.160e-01
|
||||
1.621e+00 4.545e-01
|
||||
1.460e+00 4.611e-01
|
||||
2.909e-01 4.349e-01
|
||||
6.991e-01 4.403e-01
|
||||
-5.831e-01 4.255e-01
|
||||
-2.921e-01 4.034e-01
|
||||
***************************************
|
13
pres/log
13
pres/log
@ -11,18 +11,24 @@ Given calibration data in the form of light curves produced by a known transfer
|
||||
|
||||
Did things, forgot to keep track
|
||||
|
||||
|
||||
Week 3-4ish:
|
||||
A script has been created to read an output file from psdlag and then plot the source and reprocessed power spectral densities, the cross-correlation power spectral density, and the phase difference between the signals.
|
||||
|
||||
scripts/power_lag_plot.pl
|
||||
|
||||
|
||||
2016-
|
||||
|
||||
|
||||
2016-07-11
|
||||
|
||||
created gnuplot templates for plotting the timelags and power spectrum densities
|
||||
- still need to modify this for points to be evenly spaced in log space
|
||||
|
||||
2016-07-14
|
||||
Ran through results with simple error checking and determined they made sense.
|
||||
Project can continue.
|
||||
|
||||
2016-07-15
|
||||
Determined frequency binning for PSD and timelags:
|
||||
|
||||
@ -34,3 +40,8 @@ scripts/power_lag_plot.pl
|
||||
0.16658029 0.25819945
|
||||
0.25819945 0.40020915
|
||||
0.40020915 0.62032418
|
||||
|
||||
|
||||
|
||||
2016-07-18
|
||||
Ran first MCMC run, need to analyze results
|
||||
|
@ -12,19 +12,84 @@ echo Using lightcurves in $lc_dir.
|
||||
|
||||
for echo_curve in $lc_dir/*
|
||||
do
|
||||
# Determine band and inputs for band
|
||||
echo_band=$(basename $echo_curve|sed 's@\(.*\)\.lc@\1@')
|
||||
echo -n $(date "+%R")\: Running psdlag for $echo_band against $ref_band.
|
||||
newfilename="analyses/${echo_band}≻${ref_band}_Δt=${timestep}"
|
||||
case $echo_band in
|
||||
"g(𝛌=4775)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -9.745e-01 -1.384e+00 -2.748e+00 -3.305e+00 -3.314e+00 -3.389e+00 -4.198e+00 -4.465e+00 -4.700e-01 -7.487e-01 -2.046e+00 -2.428e+00 -2.953e+00 -3.086e+00 -3.761e+00 -4.290e+00 9.862e-02 3.899e-01 8.650e-01 5.516e-01 2.228e-01 9.508e-01 -2.872e-01 9.059e-02"
|
||||
;;
|
||||
"HST(𝛌=1158)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 2.402e-01 1.329e-01 -8.990e-01 -1.229e+00 -1.825e+00 -2.236e+00 -2.758e+00 -4.749e+00 1.518e-01 1.794e-02 -1.040e+00 -1.377e+00 -1.989e+00 -2.375e+00 -2.886e+00 -3.737e+00 -3.034e-02 -6.855e-02 -6.829e-02 -4.907e-02 -5.681e-02 -1.644e-01 -3.273e-01 1.333e-01"
|
||||
;;
|
||||
"HST(𝛌=1367)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 6.534e-02 -9.689e-02 -1.175e+00 -1.525e+00 -2.167e+00 -2.492e+00 -3.258e+00 -1.115e+01 6.378e-02 -9.676e-02 -1.181e+00 -1.520e+00 -2.180e+00 -2.532e+00 -3.129e+00 -3.794e+00 -8.535e-05 -5.126e-04 -4.066e-04 1.080e-03 -2.682e-03 -5.462e-03 -2.209e-02 1.564e-01"
|
||||
;;
|
||||
"HST(𝛌=1479)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 1.224e-01 -7.607e-02 -1.174e+00 -1.500e+00 -2.075e+00 -2.554e+00 -3.337e+00 -1.152e+01 9.340e-02 -8.842e-02 -1.173e+00 -1.512e+00 -2.113e+00 -2.547e+00 -3.178e+00 -7.000e+00 -6.035e-02 2.044e-02 1.258e-01 5.764e-02 2.580e-02 -4.737e-02 9.579e-02 6.411e-01"
|
||||
;;
|
||||
"HST(𝛌=1746)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 3.934e-02 -1.638e-01 -1.339e+00 -1.655e+00 -2.232e+00 -2.793e+00 -4.025e+00 -1.072e+01 5.043e-02 -1.304e-01 -1.254e+00 -1.587e+00 -2.207e+00 -2.693e+00 -3.209e+00 -3.820e+00 3.040e-02 2.438e-02 1.450e-01 1.458e-01 1.772e-01 2.076e-01 -1.390e-01 3.929e-01"
|
||||
;;
|
||||
"i(𝛌=7647)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.510e+00 -1.872e+00 -3.441e+00 -3.722e+00 -4.794e+00 -4.336e+00 -4.644e+00 -5.061e+00 -7.402e-01 -9.967e-01 -2.393e+00 -2.723e+00 -3.701e+00 -3.742e+00 -4.317e+00 -4.979e+00 2.638e-01 7.095e-01 1.297e+00 1.035e+00 9.267e-01 8.982e-01 7.318e-01 9.092e-01"
|
||||
;;
|
||||
"I(𝛌=8560)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.497e+00 -1.786e+00 -3.555e+00 -4.203e+00 -3.826e+00 -4.096e+00 -4.304e+00 -1.109e+01 -7.542e-01 -9.532e-01 -2.397e+00 -3.214e+00 -3.008e+00 -3.489e+00 -4.353e+00 -4.061e+00 1.675e-01 7.660e-01 1.781e+00 1.971e+00 6.572e-01 8.520e-01 2.405e+00 -1.573e-01"
|
||||
;;
|
||||
"r(𝛌=6175)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.372e+00 -1.786e+00 -3.214e+00 -3.558e+00 -4.416e+00 -4.366e+00 -4.559e+00 -4.586e+00 -6.693e-01 -5.998e+00 -2.390e+00 -2.927e+00 -3.411e+00 -3.731e+00 -3.930e+00 -7.000e+00 2.905e-01 1.835e+00 8.760e-01 1.476e+00 -3.280e-01 1.439e+00 9.464e-02 1.532e+00"
|
||||
;;
|
||||
"R(𝛌=6439)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.287e+00 -1.580e+00 -2.879e+00 -3.569e+00 -3.661e+00 -4.471e+00 -3.998e+00 -5.001e+00 -6.583e-01 -6.900e+00 -2.060e+00 -2.727e+00 -2.954e+00 -3.654e+00 -4.406e+00 -5.203e+00 1.448e-01 -1.004e-01 9.754e-01 1.140e+00 9.284e-01 1.116e+00 2.092e+00 8.634e-01"
|
||||
;;
|
||||
"Swiftb(𝛌=4392)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -5.521e-01 -1.135e+00 -2.691e+00 -3.626e+00 -4.510e+00 -1.036e+01 -1.076e+01 -1.213e+01 -2.515e-01 -6.156e-01 -1.924e+00 -2.417e+00 -3.053e+00 -6.570e+00 -3.450e+00 -6.333e+00 7.408e-02 3.089e-01 5.057e-01 4.831e-01 5.779e-01 6.277e-01 4.823e-01 6.110e-01"
|
||||
;;
|
||||
"Swiftu(𝛌=3465)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -3.178e-01 -7.537e-01 -2.402e+00 -2.624e+00 -3.413e+00 -3.834e+00 -1.236e+01 -1.071e+01 -1.376e-01 -4.256e-01 -1.794e+00 -2.123e+00 -2.817e+00 -3.029e+00 -4.139e+00 -7.000e+00 7.840e-02 3.388e-01 6.616e-01 1.349e-01 9.909e-01 8.538e-01 1.397e+00 6.412e-01"
|
||||
;;
|
||||
"Swiftv(𝛌=5468)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -8.802e-01 -1.441e+00 -2.883e+00 -4.483e+00 -1.045e+01 -1.041e+01 -1.042e+01 -1.078e+01 -4.161e-01 -7.660e-01 -2.020e+00 -3.201e+00 -2.699e+00 -2.859e+00 -3.653e+00 -5.501e+00 1.900e-02 3.560e-01 1.007e+00 -1.270e+00 1.858e+00 9.302e-01 -1.502e+00 2.653e+00"
|
||||
;;
|
||||
"u(𝛌=3471)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -4.289e-01 -9.017e-01 -2.581e+00 -2.751e+00 -2.764e+00 -3.105e+00 -3.298e+00 -3.870e+00 -1.925e-01 -5.093e-01 -1.923e+00 -2.204e+00 -2.568e+00 -2.768e+00 -3.425e+00 -5.982e+00 1.646e-01 4.376e-01 7.202e-01 9.556e-01 4.856e-01 2.985e-01 -6.801e-01 -2.416e+00"
|
||||
;;
|
||||
"UVM2(𝛌=2246)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.003e-01 -5.235e-01 -1.863e+00 -2.159e+00 -2.718e+00 -3.504e+00 -1.139e+01 -4.119e+00 -2.431e-02 -3.121e-01 -1.526e+00 -1.877e+00 -2.455e+00 -2.945e+00 -3.356e+00 -3.602e+00 -1.862e-04 1.380e-01 2.475e-01 3.102e-02 6.040e-01 -1.488e-01 1.171e+00 1.617e+00"
|
||||
;;
|
||||
"UVW1(𝛌=2600)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -2.417e-01 -7.027e-01 -1.998e+00 -2.387e+00 -2.912e+00 -4.026e+00 -1.069e+01 -9.716e+00 -9.867e-02 -4.150e-01 -1.587e+00 -1.963e+00 -2.539e+00 -3.274e+00 -3.460e+00 -6.742e+00 6.734e-02 1.353e-01 3.585e-01 1.587e-01 2.959e-01 1.111e-01 2.364e-01 7.463e-01"
|
||||
;;
|
||||
"UVW2(𝛌=1928)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.648e-01 -4.974e-01 -1.686e+00 -1.897e+00 -2.600e+00 -2.801e+00 -3.814e+00 -9.888e+00 -6.101e-02 -3.020e-01 -1.433e+00 -1.756e+00 -2.485e+00 -2.799e+00 -3.317e+00 -7.000e+00 4.622e-02 1.473e-01 2.243e-01 1.606e-01 2.365e-01 1.484e-01 8.110e-01 9.689e-01"
|
||||
;;
|
||||
"V(𝛌=5404)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.044e+00 -1.457e+00 -2.848e+00 -2.957e+00 -4.232e+00 -4.052e+00 -4.515e+00 -5.250e+00 -5.026e-01 -7.878e-01 -2.118e+00 -2.276e+00 -3.421e+00 -3.378e+00 -4.224e+00 -4.930e+00 9.142e-02 4.814e-01 8.550e-01 3.929e-01 3.975e-01 5.311e-01 -1.203e-01 -2.905e+00"
|
||||
;;
|
||||
"z(𝛌=9157)")
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.852e+00 -2.156e+00 -3.798e+00 -4.146e+00 -4.577e+00 -4.338e+00 -4.672e+00 -5.230e+00 -9.339e-01 -1.131e+00 -2.536e+00 -3.031e+00 -4.060e+00 -3.801e+00 -4.139e+00 -5.151e+00 2.445e-01 8.131e-01 1.621e+00 1.460e+00 2.909e-01 6.991e-01 -5.831e-01 -2.921e-01"
|
||||
;;
|
||||
*)
|
||||
echo "Unknown band: using default parameters."
|
||||
initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -9.745e-01 -1.384e+00 -2.748e+00 -3.305e+00 -3.314e+00 -3.389e+00 -4.198e+00 -4.465e+00 -4.700e-01 -7.487e-01 -2.046e+00 -2.428e+00 -2.953e+00 -3.086e+00 -3.761e+00 -4.290e+00 9.862e-02 3.899e-01 8.650e-01 5.516e-01 2.228e-01 9.508e-01 -2.872e-01 9.059e-02"
|
||||
;;
|
||||
esac
|
||||
|
||||
# Write temporary input file
|
||||
echo "2" > tmp.psdlagargs
|
||||
echo "$ref_curve 0" >> tmp.psdlagargs
|
||||
echo "$echo_curve 0" >> tmp.psdlagargs
|
||||
echo "8 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603" >> tmp.psdlagargs
|
||||
echo "0" >> tmp.psdlagargs
|
||||
echo ".9999 .9999 .9999 0.49 0.42 0.30 0.36 1.2 -0.60 -0.89 -0.72 2.39 -0.74 -1.5 -1 -3.1 -1.3 -2.3 -2 -1.9 -1.6 -2.6 -2 -1.3 -1.8 -2.6 -1.2 0 -2 -2.7 -1.3 0.63" >> tmp.psdlagargs
|
||||
echo "$initial_params" >> tmp.psdlagargs
|
||||
echo "0:0 0" >> tmp.psdlagargs
|
||||
echo $error_type >> tmp.psdlagargs
|
||||
echo "0" >> tmp.psdlagargs
|
||||
echo "100 50 50 mcmc.dat" >> tmp.psdlagargs
|
||||
|
||||
# Run psdlag with inputs
|
||||
time bin/psdlag tmp.psdlagargs >> $newfilename
|
||||
echo ""
|
||||
done
|
Loading…
Reference in New Issue
Block a user