Commit 5e85f6a6 authored by lvs's avatar lvs
Browse files

added grad correction for axisangle

parent d251789a
......@@ -206,11 +206,11 @@ PDB := 3ce3
NRUN := 100
POPSIZE := 150
TESTNAME := test
TESTLS := fire
TESTLS := sd
test: odock
$(BIN_DIR)/$(TARGET) -ffile ./input/$(PDB)/derived/$(PDB)_protein.maps.fld -lfile ./input/$(PDB)/derived/$(PDB)_ligand.pdbqt -nrun $(NRUN) -psize $(POPSIZE) -resnam $(TESTNAME) -gfpop 0 -lsmet $(TESTLS) -lsit 10
$(BIN_DIR)/$(TARGET) -ffile ./input/$(PDB)/derived/$(PDB)_protein.maps.fld -lfile ./input/$(PDB)/derived/$(PDB)_ligand.pdbqt -nrun $(NRUN) -psize $(POPSIZE) -resnam $(TESTNAME) -gfpop 0 -lsmet $(TESTLS) -lsit 30
ASTEX_PDB := 2bsm
ASTEX_NRUN:= 10
......
......@@ -1086,11 +1086,11 @@ void gpu_calc_gradient(
// Using interpolation on out-of-bounds elements results in hang
if (index_theta <= 0) {
//printf("WARNING: index_theta: %u\n", index_theta);
dependence_on_theta = dependence_on_theta_const[0];
dependence_on_theta = dependence_on_theta_const[0]; //printf("%f\n",dependence_on_theta_const[0]);
}
else if (index_theta >= 999){
//printf("WARNING: index_theta: %u\n", index_theta);
dependence_on_theta = dependence_on_theta_const[999];
dependence_on_theta = dependence_on_theta_const[999]; //printf("%f\n",dependence_on_theta_const[999]);
}
else {
X0_theta = angle_const[index_theta];
......@@ -1110,11 +1110,11 @@ void gpu_calc_gradient(
// Using interpolation on previous and/or next elements results in hang
if (index_rotangle <= 0) {
//printf("WARNING: index_rotangle: %u\n", index_rotangle);
dependence_on_rotangle = dependence_on_rotangle_const[0];
dependence_on_rotangle = dependence_on_rotangle_const[0]; //printf("%f\n",dependence_on_rotangle_const[0]);
}
else if (index_rotangle >= 999){
//printf("WARNING: index_rotangle: %u\n", index_rotangle);
dependence_on_rotangle = dependence_on_rotangle_const[999];
dependence_on_rotangle = dependence_on_rotangle_const[999]; //printf("%f\n",dependence_on_rotangle_const[999]);
}
else {
X0_rotangle = angle_const[index_rotangle];
......
/*
Gradient corrections
The corrections for dE/d_phi depend on two variables: rotangle and theta.
The correction for dE/d_theta depends on a single variable: rotangle.
These arrays are passed to constant device on device.
*/
/*const*/ float angle [1000] = {
3.14159265e-03, 9.42477796e-03, 1.57079633e-02,
2.19911486e-02, 2.82743339e-02, 3.45575192e-02,
4.08407045e-02, 4.71238898e-02, 5.34070751e-02,
5.96902604e-02, 6.59734457e-02, 7.22566310e-02,
7.85398163e-02, 8.48230016e-02, 9.11061870e-02,
9.73893723e-02, 1.03672558e-01, 1.09955743e-01,
1.16238928e-01, 1.22522113e-01, 1.28805299e-01,
1.35088484e-01, 1.41371669e-01, 1.47654855e-01,
1.53938040e-01, 1.60221225e-01, 1.66504411e-01,
1.72787596e-01, 1.79070781e-01, 1.85353967e-01,
1.91637152e-01, 1.97920337e-01, 2.04203522e-01,
2.10486708e-01, 2.16769893e-01, 2.23053078e-01,
2.29336264e-01, 2.35619449e-01, 2.41902634e-01,
2.48185820e-01, 2.54469005e-01, 2.60752190e-01,
2.67035376e-01, 2.73318561e-01, 2.79601746e-01,
2.85884931e-01, 2.92168117e-01, 2.98451302e-01,
3.04734487e-01, 3.11017673e-01, 3.17300858e-01,
3.23584043e-01, 3.29867229e-01, 3.36150414e-01,
3.42433599e-01, 3.48716785e-01, 3.54999970e-01,
3.61283155e-01, 3.67566340e-01, 3.73849526e-01,
3.80132711e-01, 3.86415896e-01, 3.92699082e-01,
3.98982267e-01, 4.05265452e-01, 4.11548638e-01,
4.17831823e-01, 4.24115008e-01, 4.30398194e-01,
4.36681379e-01, 4.42964564e-01, 4.49247749e-01,
4.55530935e-01, 4.61814120e-01, 4.68097305e-01,
4.74380491e-01, 4.80663676e-01, 4.86946861e-01,
4.93230047e-01, 4.99513232e-01, 5.05796417e-01,
5.12079603e-01, 5.18362788e-01, 5.24645973e-01,
5.30929158e-01, 5.37212344e-01, 5.43495529e-01,
5.49778714e-01, 5.56061900e-01, 5.62345085e-01,
5.68628270e-01, 5.74911456e-01, 5.81194641e-01,
5.87477826e-01, 5.93761012e-01, 6.00044197e-01,
6.06327382e-01, 6.12610567e-01, 6.18893753e-01,
6.25176938e-01, 6.31460123e-01, 6.37743309e-01,
6.44026494e-01, 6.50309679e-01, 6.56592865e-01,
6.62876050e-01, 6.69159235e-01, 6.75442421e-01,
6.81725606e-01, 6.88008791e-01, 6.94291976e-01,
7.00575162e-01, 7.06858347e-01, 7.13141532e-01,
7.19424718e-01, 7.25707903e-01, 7.31991088e-01,
7.38274274e-01, 7.44557459e-01, 7.50840644e-01,
7.57123830e-01, 7.63407015e-01, 7.69690200e-01,
7.75973385e-01, 7.82256571e-01, 7.88539756e-01,
7.94822941e-01, 8.01106127e-01, 8.07389312e-01,
8.13672497e-01, 8.19955683e-01, 8.26238868e-01,
8.32522053e-01, 8.38805239e-01, 8.45088424e-01,
8.51371609e-01, 8.57654794e-01, 8.63937980e-01,
8.70221165e-01, 8.76504350e-01, 8.82787536e-01,
8.89070721e-01, 8.95353906e-01, 9.01637092e-01,
9.07920277e-01, 9.14203462e-01, 9.20486648e-01,
9.26769833e-01, 9.33053018e-01, 9.39336203e-01,
9.45619389e-01, 9.51902574e-01, 9.58185759e-01,
9.64468945e-01, 9.70752130e-01, 9.77035315e-01,
9.83318501e-01, 9.89601686e-01, 9.95884871e-01,
1.00216806e+00, 1.00845124e+00, 1.01473443e+00,
1.02101761e+00, 1.02730080e+00, 1.03358398e+00,
1.03986717e+00, 1.04615035e+00, 1.05243354e+00,
1.05871672e+00, 1.06499991e+00, 1.07128309e+00,
1.07756628e+00, 1.08384947e+00, 1.09013265e+00,
1.09641584e+00, 1.10269902e+00, 1.10898221e+00,
1.11526539e+00, 1.12154858e+00, 1.12783176e+00,
1.13411495e+00, 1.14039813e+00, 1.14668132e+00,
1.15296450e+00, 1.15924769e+00, 1.16553087e+00,
1.17181406e+00, 1.17809725e+00, 1.18438043e+00,
1.19066362e+00, 1.19694680e+00, 1.20322999e+00,
1.20951317e+00, 1.21579636e+00, 1.22207954e+00,
1.22836273e+00, 1.23464591e+00, 1.24092910e+00,
1.24721228e+00, 1.25349547e+00, 1.25977865e+00,
1.26606184e+00, 1.27234502e+00, 1.27862821e+00,
1.28491140e+00, 1.29119458e+00, 1.29747777e+00,
1.30376095e+00, 1.31004414e+00, 1.31632732e+00,
1.32261051e+00, 1.32889369e+00, 1.33517688e+00,
1.34146006e+00, 1.34774325e+00, 1.35402643e+00,
1.36030962e+00, 1.36659280e+00, 1.37287599e+00,
1.37915917e+00, 1.38544236e+00, 1.39172555e+00,
1.39800873e+00, 1.40429192e+00, 1.41057510e+00,
1.41685829e+00, 1.42314147e+00, 1.42942466e+00,
1.43570784e+00, 1.44199103e+00, 1.44827421e+00,
1.45455740e+00, 1.46084058e+00, 1.46712377e+00,
1.47340695e+00, 1.47969014e+00, 1.48597333e+00,
1.49225651e+00, 1.49853970e+00, 1.50482288e+00,
1.51110607e+00, 1.51738925e+00, 1.52367244e+00,
1.52995562e+00, 1.53623881e+00, 1.54252199e+00,
1.54880518e+00, 1.55508836e+00, 1.56137155e+00,
1.56765473e+00, 1.57393792e+00, 1.58022110e+00,
1.58650429e+00, 1.59278748e+00, 1.59907066e+00,
1.60535385e+00, 1.61163703e+00, 1.61792022e+00,
1.62420340e+00, 1.63048659e+00, 1.63676977e+00,
1.64305296e+00, 1.64933614e+00, 1.65561933e+00,
1.66190251e+00, 1.66818570e+00, 1.67446888e+00,
1.68075207e+00, 1.68703525e+00, 1.69331844e+00,
1.69960163e+00, 1.70588481e+00, 1.71216800e+00,
1.71845118e+00, 1.72473437e+00, 1.73101755e+00,
1.73730074e+00, 1.74358392e+00, 1.74986711e+00,
1.75615029e+00, 1.76243348e+00, 1.76871666e+00,
1.77499985e+00, 1.78128303e+00, 1.78756622e+00,
1.79384941e+00, 1.80013259e+00, 1.80641578e+00,
1.81269896e+00, 1.81898215e+00, 1.82526533e+00,
1.83154852e+00, 1.83783170e+00, 1.84411489e+00,
1.85039807e+00, 1.85668126e+00, 1.86296444e+00,
1.86924763e+00, 1.87553081e+00, 1.88181400e+00,
1.88809718e+00, 1.89438037e+00, 1.90066356e+00,
1.90694674e+00, 1.91322993e+00, 1.91951311e+00,
1.92579630e+00, 1.93207948e+00, 1.93836267e+00,
1.94464585e+00, 1.95092904e+00, 1.95721222e+00,
1.96349541e+00, 1.96977859e+00, 1.97606178e+00,
1.98234496e+00, 1.98862815e+00, 1.99491134e+00,
2.00119452e+00, 2.00747771e+00, 2.01376089e+00,
2.02004408e+00, 2.02632726e+00, 2.03261045e+00,
2.03889363e+00, 2.04517682e+00, 2.05146000e+00,
2.05774319e+00, 2.06402637e+00, 2.07030956e+00,
2.07659274e+00, 2.08287593e+00, 2.08915911e+00,
2.09544230e+00, 2.10172549e+00, 2.10800867e+00,
2.11429186e+00, 2.12057504e+00, 2.12685823e+00,
2.13314141e+00, 2.13942460e+00, 2.14570778e+00,
2.15199097e+00, 2.15827415e+00, 2.16455734e+00,
2.17084052e+00, 2.17712371e+00, 2.18340689e+00,
2.18969008e+00, 2.19597326e+00, 2.20225645e+00,
2.20853964e+00, 2.21482282e+00, 2.22110601e+00,
2.22738919e+00, 2.23367238e+00, 2.23995556e+00,
2.24623875e+00, 2.25252193e+00, 2.25880512e+00,
2.26508830e+00, 2.27137149e+00, 2.27765467e+00,
2.28393786e+00, 2.29022104e+00, 2.29650423e+00,
2.30278742e+00, 2.30907060e+00, 2.31535379e+00,
2.32163697e+00, 2.32792016e+00, 2.33420334e+00,
2.34048653e+00, 2.34676971e+00, 2.35305290e+00,
2.35933608e+00, 2.36561927e+00, 2.37190245e+00,
2.37818564e+00, 2.38446882e+00, 2.39075201e+00,
2.39703519e+00, 2.40331838e+00, 2.40960157e+00,
2.41588475e+00, 2.42216794e+00, 2.42845112e+00,
2.43473431e+00, 2.44101749e+00, 2.44730068e+00,
2.45358386e+00, 2.45986705e+00, 2.46615023e+00,
2.47243342e+00, 2.47871660e+00, 2.48499979e+00,
2.49128297e+00, 2.49756616e+00, 2.50384934e+00,
2.51013253e+00, 2.51641572e+00, 2.52269890e+00,
2.52898209e+00, 2.53526527e+00, 2.54154846e+00,
2.54783164e+00, 2.55411483e+00, 2.56039801e+00,
2.56668120e+00, 2.57296438e+00, 2.57924757e+00,
2.58553075e+00, 2.59181394e+00, 2.59809712e+00,
2.60438031e+00, 2.61066350e+00, 2.61694668e+00,
2.62322987e+00, 2.62951305e+00, 2.63579624e+00,
2.64207942e+00, 2.64836261e+00, 2.65464579e+00,
2.66092898e+00, 2.66721216e+00, 2.67349535e+00,
2.67977853e+00, 2.68606172e+00, 2.69234490e+00,
2.69862809e+00, 2.70491127e+00, 2.71119446e+00,
2.71747765e+00, 2.72376083e+00, 2.73004402e+00,
2.73632720e+00, 2.74261039e+00, 2.74889357e+00,
2.75517676e+00, 2.76145994e+00, 2.76774313e+00,
2.77402631e+00, 2.78030950e+00, 2.78659268e+00,
2.79287587e+00, 2.79915905e+00, 2.80544224e+00,
2.81172542e+00, 2.81800861e+00, 2.82429180e+00,
2.83057498e+00, 2.83685817e+00, 2.84314135e+00,
2.84942454e+00, 2.85570772e+00, 2.86199091e+00,
2.86827409e+00, 2.87455728e+00, 2.88084046e+00,
2.88712365e+00, 2.89340683e+00, 2.89969002e+00,
2.90597320e+00, 2.91225639e+00, 2.91853958e+00,
2.92482276e+00, 2.93110595e+00, 2.93738913e+00,
2.94367232e+00, 2.94995550e+00, 2.95623869e+00,
2.96252187e+00, 2.96880506e+00, 2.97508824e+00,
2.98137143e+00, 2.98765461e+00, 2.99393780e+00,
3.00022098e+00, 3.00650417e+00, 3.01278735e+00,
3.01907054e+00, 3.02535373e+00, 3.03163691e+00,
3.03792010e+00, 3.04420328e+00, 3.05048647e+00,
3.05676965e+00, 3.06305284e+00, 3.06933602e+00,
3.07561921e+00, 3.08190239e+00, 3.08818558e+00,
3.09446876e+00, 3.10075195e+00, 3.10703513e+00,
3.11331832e+00, 3.11960151e+00, 3.12588469e+00,
3.13216788e+00, 3.13845106e+00, 3.14473425e+00,
3.15101743e+00, 3.15730062e+00, 3.16358380e+00,
3.16986699e+00, 3.17615017e+00, 3.18243336e+00,
3.18871654e+00, 3.19499973e+00, 3.20128291e+00,
3.20756610e+00, 3.21384928e+00, 3.22013247e+00,
3.22641566e+00, 3.23269884e+00, 3.23898203e+00,
3.24526521e+00, 3.25154840e+00, 3.25783158e+00,
3.26411477e+00, 3.27039795e+00, 3.27668114e+00,
3.28296432e+00, 3.28924751e+00, 3.29553069e+00,
3.30181388e+00, 3.30809706e+00, 3.31438025e+00,
3.32066343e+00, 3.32694662e+00, 3.33322981e+00,
3.33951299e+00, 3.34579618e+00, 3.35207936e+00,
3.35836255e+00, 3.36464573e+00, 3.37092892e+00,
3.37721210e+00, 3.38349529e+00, 3.38977847e+00,
3.39606166e+00, 3.40234484e+00, 3.40862803e+00,
3.41491121e+00, 3.42119440e+00, 3.42747759e+00,
3.43376077e+00, 3.44004396e+00, 3.44632714e+00,
3.45261033e+00, 3.45889351e+00, 3.46517670e+00,
3.47145988e+00, 3.47774307e+00, 3.48402625e+00,
3.49030944e+00, 3.49659262e+00, 3.50287581e+00,
3.50915899e+00, 3.51544218e+00, 3.52172536e+00,
3.52800855e+00, 3.53429174e+00, 3.54057492e+00,
3.54685811e+00, 3.55314129e+00, 3.55942448e+00,
3.56570766e+00, 3.57199085e+00, 3.57827403e+00,
3.58455722e+00, 3.59084040e+00, 3.59712359e+00,
3.60340677e+00, 3.60968996e+00, 3.61597314e+00,
3.62225633e+00, 3.62853951e+00, 3.63482270e+00,
3.64110589e+00, 3.64738907e+00, 3.65367226e+00,
3.65995544e+00, 3.66623863e+00, 3.67252181e+00,
3.67880500e+00, 3.68508818e+00, 3.69137137e+00,
3.69765455e+00, 3.70393774e+00, 3.71022092e+00,
3.71650411e+00, 3.72278729e+00, 3.72907048e+00,
3.73535367e+00, 3.74163685e+00, 3.74792004e+00,
3.75420322e+00, 3.76048641e+00, 3.76676959e+00,
3.77305278e+00, 3.77933596e+00, 3.78561915e+00,
3.79190233e+00, 3.79818552e+00, 3.80446870e+00,
3.81075189e+00, 3.81703507e+00, 3.82331826e+00,
3.82960144e+00, 3.83588463e+00, 3.84216782e+00,
3.84845100e+00, 3.85473419e+00, 3.86101737e+00,
3.86730056e+00, 3.87358374e+00, 3.87986693e+00,
3.88615011e+00, 3.89243330e+00, 3.89871648e+00,
3.90499967e+00, 3.91128285e+00, 3.91756604e+00,
3.92384922e+00, 3.93013241e+00, 3.93641559e+00,
3.94269878e+00, 3.94898197e+00, 3.95526515e+00,
3.96154834e+00, 3.96783152e+00, 3.97411471e+00,
3.98039789e+00, 3.98668108e+00, 3.99296426e+00,
3.99924745e+00, 4.00553063e+00, 4.01181382e+00,
4.01809700e+00, 4.02438019e+00, 4.03066337e+00,
4.03694656e+00, 4.04322975e+00, 4.04951293e+00,
4.05579612e+00, 4.06207930e+00, 4.06836249e+00,
4.07464567e+00, 4.08092886e+00, 4.08721204e+00,
4.09349523e+00, 4.09977841e+00, 4.10606160e+00,
4.11234478e+00, 4.11862797e+00, 4.12491115e+00,
4.13119434e+00, 4.13747752e+00, 4.14376071e+00,
4.15004390e+00, 4.15632708e+00, 4.16261027e+00,
4.16889345e+00, 4.17517664e+00, 4.18145982e+00,
4.18774301e+00, 4.19402619e+00, 4.20030938e+00,
4.20659256e+00, 4.21287575e+00, 4.21915893e+00,
4.22544212e+00, 4.23172530e+00, 4.23800849e+00,
4.24429167e+00, 4.25057486e+00, 4.25685805e+00,
4.26314123e+00, 4.26942442e+00, 4.27570760e+00,
4.28199079e+00, 4.28827397e+00, 4.29455716e+00,
4.30084034e+00, 4.30712353e+00, 4.31340671e+00,
4.31968990e+00, 4.32597308e+00, 4.33225627e+00,
4.33853945e+00, 4.34482264e+00, 4.35110583e+00,
4.35738901e+00, 4.36367220e+00, 4.36995538e+00,
4.37623857e+00, 4.38252175e+00, 4.38880494e+00,
4.39508812e+00, 4.40137131e+00, 4.40765449e+00,
4.41393768e+00, 4.42022086e+00, 4.42650405e+00,
4.43278723e+00, 4.43907042e+00, 4.44535360e+00,
4.45163679e+00, 4.45791998e+00, 4.46420316e+00,
4.47048635e+00, 4.47676953e+00, 4.48305272e+00,
4.48933590e+00, 4.49561909e+00, 4.50190227e+00,
4.50818546e+00, 4.51446864e+00, 4.52075183e+00,
4.52703501e+00, 4.53331820e+00, 4.53960138e+00,
4.54588457e+00, 4.55216776e+00, 4.55845094e+00,
4.56473413e+00, 4.57101731e+00, 4.57730050e+00,
4.58358368e+00, 4.58986687e+00, 4.59615005e+00,
4.60243324e+00, 4.60871642e+00, 4.61499961e+00,
4.62128279e+00, 4.62756598e+00, 4.63384916e+00,
4.64013235e+00, 4.64641553e+00, 4.65269872e+00,
4.65898191e+00, 4.66526509e+00, 4.67154828e+00,
4.67783146e+00, 4.68411465e+00, 4.69039783e+00,
4.69668102e+00, 4.70296420e+00, 4.70924739e+00,
4.71553057e+00, 4.72181376e+00, 4.72809694e+00,
4.73438013e+00, 4.74066331e+00, 4.74694650e+00,
4.75322968e+00, 4.75951287e+00, 4.76579606e+00,
4.77207924e+00, 4.77836243e+00, 4.78464561e+00,
4.79092880e+00, 4.79721198e+00, 4.80349517e+00,
4.80977835e+00, 4.81606154e+00, 4.82234472e+00,
4.82862791e+00, 4.83491109e+00, 4.84119428e+00,
4.84747746e+00, 4.85376065e+00, 4.86004384e+00,
4.86632702e+00, 4.87261021e+00, 4.87889339e+00,
4.88517658e+00, 4.89145976e+00, 4.89774295e+00,
4.90402613e+00, 4.91030932e+00, 4.91659250e+00,
4.92287569e+00, 4.92915887e+00, 4.93544206e+00,
4.94172524e+00, 4.94800843e+00, 4.95429161e+00,
4.96057480e+00, 4.96685799e+00, 4.97314117e+00,
4.97942436e+00, 4.98570754e+00, 4.99199073e+00,
4.99827391e+00, 5.00455710e+00, 5.01084028e+00,
5.01712347e+00, 5.02340665e+00, 5.02968984e+00,
5.03597302e+00, 5.04225621e+00, 5.04853939e+00,
5.05482258e+00, 5.06110576e+00, 5.06738895e+00,
5.07367214e+00, 5.07995532e+00, 5.08623851e+00,
5.09252169e+00, 5.09880488e+00, 5.10508806e+00,
5.11137125e+00, 5.11765443e+00, 5.12393762e+00,
5.13022080e+00, 5.13650399e+00, 5.14278717e+00,
5.14907036e+00, 5.15535354e+00, 5.16163673e+00,
5.16791992e+00, 5.17420310e+00, 5.18048629e+00,
5.18676947e+00, 5.19305266e+00, 5.19933584e+00,
5.20561903e+00, 5.21190221e+00, 5.21818540e+00,
5.22446858e+00, 5.23075177e+00, 5.23703495e+00,
5.24331814e+00, 5.24960132e+00, 5.25588451e+00,
5.26216769e+00, 5.26845088e+00, 5.27473407e+00,
5.28101725e+00, 5.28730044e+00, 5.29358362e+00,
5.29986681e+00, 5.30614999e+00, 5.31243318e+00,
5.31871636e+00, 5.32499955e+00, 5.33128273e+00,
5.33756592e+00, 5.34384910e+00, 5.35013229e+00,
5.35641547e+00, 5.36269866e+00, 5.36898184e+00,
5.37526503e+00, 5.38154822e+00, 5.38783140e+00,
5.39411459e+00, 5.40039777e+00, 5.40668096e+00,
5.41296414e+00, 5.41924733e+00, 5.42553051e+00,
5.43181370e+00, 5.43809688e+00, 5.44438007e+00,
5.45066325e+00, 5.45694644e+00, 5.46322962e+00,
5.46951281e+00, 5.47579600e+00, 5.48207918e+00,
5.48836237e+00, 5.49464555e+00, 5.50092874e+00,
5.50721192e+00, 5.51349511e+00, 5.51977829e+00,
5.52606148e+00, 5.53234466e+00, 5.53862785e+00,
5.54491103e+00, 5.55119422e+00, 5.55747740e+00,
5.56376059e+00, 5.57004377e+00, 5.57632696e+00,
5.58261015e+00, 5.58889333e+00, 5.59517652e+00,
5.60145970e+00, 5.60774289e+00, 5.61402607e+00,
5.62030926e+00, 5.62659244e+00, 5.63287563e+00,
5.63915881e+00, 5.64544200e+00, 5.65172518e+00,
5.65800837e+00, 5.66429155e+00, 5.67057474e+00,
5.67685793e+00, 5.68314111e+00, 5.68942430e+00,
5.69570748e+00, 5.70199067e+00, 5.70827385e+00,
5.71455704e+00, 5.72084022e+00, 5.72712341e+00,
5.73340659e+00, 5.73968978e+00, 5.74597296e+00,
5.75225615e+00, 5.75853933e+00, 5.76482252e+00,
5.77110570e+00, 5.77738889e+00, 5.78367208e+00,
5.78995526e+00, 5.79623845e+00, 5.80252163e+00,
5.80880482e+00, 5.81508800e+00, 5.82137119e+00,
5.82765437e+00, 5.83393756e+00, 5.84022074e+00,
5.84650393e+00, 5.85278711e+00, 5.85907030e+00,
5.86535348e+00, 5.87163667e+00, 5.87791985e+00,
5.88420304e+00, 5.89048623e+00, 5.89676941e+00,
5.90305260e+00, 5.90933578e+00, 5.91561897e+00,
5.92190215e+00, 5.92818534e+00, 5.93446852e+00,
5.94075171e+00, 5.94703489e+00, 5.95331808e+00,
5.95960126e+00, 5.96588445e+00, 5.97216763e+00,
5.97845082e+00, 5.98473401e+00, 5.99101719e+00,
5.99730038e+00, 6.00358356e+00, 6.00986675e+00,
6.01614993e+00, 6.02243312e+00, 6.02871630e+00,
6.03499949e+00, 6.04128267e+00, 6.04756586e+00,
6.05384904e+00, 6.06013223e+00, 6.06641541e+00,
6.07269860e+00, 6.07898178e+00, 6.08526497e+00,
6.09154816e+00, 6.09783134e+00, 6.10411453e+00,
6.11039771e+00, 6.11668090e+00, 6.12296408e+00,
6.12924727e+00, 6.13553045e+00, 6.14181364e+00,
6.14809682e+00, 6.15438001e+00, 6.16066319e+00,
6.16694638e+00, 6.17322956e+00, 6.17951275e+00,
6.18579593e+00, 6.19207912e+00, 6.19836231e+00,
6.20464549e+00, 6.21092868e+00, 6.21721186e+00,
6.22349505e+00, 6.22977823e+00, 6.23606142e+00,
6.24234460e+00, 6.24862779e+00, 6.25491097e+00,
6.26119416e+00, 6.26747734e+00, 6.27376053e+00,
6.28004371e+00};
/*const*/ float dependence_on_theta [1000] = {
1.01328679e+05, 1.12547638e+04, 4.05232581e+03,
2.06779540e+03, 1.25114454e+03, 8.37717298e+02,
5.99884079e+02, 4.50639506e+02, 3.50909153e+02,
2.80987869e+02, 2.30079555e+02, 1.91865997e+02,
1.62447378e+02, 1.39323295e+02, 1.20810744e+02,
1.05768556e+02, 9.33790456e+01, 8.30488939e+01,
7.43478122e+01, 6.69503661e+01, 6.06086671e+01,
5.51328881e+01, 5.03690865e+01, 4.62009357e+01,
4.25334622e+01, 3.92888692e+01, 3.64046398e+01,
3.38293598e+01, 3.15201220e+01, 2.94418547e+01,
2.75646975e+01, 2.58634663e+01, 2.43169068e+01,
2.29068181e+01, 2.16175522e+01, 2.04357596e+01,
1.93497322e+01, 1.83495933e+01, 1.74262607e+01,
1.65721732e+01, 1.57804920e+01, 1.50454362e+01,
1.43619076e+01, 1.37249896e+01, 1.31302931e+01,
1.25744454e+01, 1.20540522e+01, 1.15662457e+01,
1.11083211e+01, 1.06778297e+01, 1.02725643e+01,
9.89095962e+00, 9.53098292e+00, 9.19091021e+00,
8.86941981e+00, 8.56516420e+00, 8.27683055e+00,
8.00351016e+00, 7.74417395e+00, 7.49783768e+00,
7.26358224e+00, 7.04068327e+00, 6.82843071e+00,
6.62617100e+00, 6.43323563e+00, 6.24910558e+00,
6.07322479e+00, 5.90511845e+00, 5.74435998e+00,
5.59054492e+00, 5.44322874e+00, 5.30207084e+00,
5.16671513e+00, 5.03688962e+00, 4.91229002e+00,
4.79262152e+00, 4.67763433e+00, 4.56708629e+00,
4.46074531e+00, 4.35843392e+00, 4.25994365e+00,
4.16507787e+00, 4.07366131e+00, 3.98555075e+00,
3.90055547e+00, 3.81855434e+00, 3.73939582e+00,
3.66296101e+00, 3.58912508e+00, 3.51775796e+00,
3.44877483e+00, 3.38207528e+00, 3.31755731e+00,
3.25512820e+00, 3.19469147e+00, 3.13617023e+00,
3.07947700e+00, 3.02453649e+00, 2.97128591e+00,
2.91966199e+00, 2.86959943e+00, 2.82103199e+00,
2.77390801e+00, 2.72816702e+00, 2.68374666e+00,
2.64059954e+00, 2.59870176e+00, 2.55800095e+00,
2.51844312e+00, 2.47999224e+00, 2.44260567e+00,
2.40624915e+00, 2.37088649e+00, 2.33648495e+00,
2.30300222e+00, 2.27040763e+00, 2.23868041e+00,
2.20778423e+00, 2.17768903e+00, 2.14837859e+00,
2.11981847e+00, 2.09198790e+00, 2.06486401e+00,
2.03842543e+00, 2.01264338e+00, 1.98751345e+00,
1.96300108e+00, 1.93908987e+00, 1.91576471e+00,
1.89300774e+00, 1.87080063e+00, 1.84912465e+00,
1.82796755e+00, 1.80731332e+00, 1.78714983e+00,
1.76745921e+00, 1.74822437e+00, 1.72944845e+00,
1.71110419e+00, 1.69318232e+00, 1.67567214e+00,
1.65856508e+00, 1.64185374e+00, 1.62550144e+00,
1.60952917e+00, 1.59391742e+00, 1.57865499e+00,
1.56373223e+00, 1.54914103e+00, 1.53487630e+00,
1.52092919e+00, 1.50728597e+00, 1.49394282e+00,
1.48088649e+00, 1.46813113e+00, 1.45564755e+00,
1.44343447e+00, 1.43148805e+00, 1.41979994e+00,
1.40836624e+00, 1.39718117e+00, 1.38623693e+00,
1.37552794e+00, 1.36504952e+00, 1.35479612e+00,
1.34476443e+00, 1.33495009e+00, 1.32534595e+00,
1.31595004e+00, 1.30675442e+00, 1.29775655e+00,
1.28895261e+00, 1.28034123e+00, 1.27191563e+00,
1.26367372e+00, 1.25560738e+00, 1.24771840e+00,
1.23999876e+00, 1.23244970e+00, 1.22506494e+00,
1.21784158e+00, 1.21077672e+00, 1.20386854e+00,
1.19711371e+00, 1.19050870e+00, 1.18405310e+00,
1.17774293e+00, 1.17157442e+00, 1.16554523e+00,
1.15965446e+00, 1.15389859e+00, 1.14827612e+00,
1.14278441e+00, 1.13742036e+00, 1.13218695e+00,
1.12707674e+00, 1.12208979e+00, 1.11722327e+00,
1.11247469e+00, 1.10784488e+00, 1.10333019e+00,
1.09893053e+00, 1.09464302e+00, 1.09046596e+00,
1.08639874e+00, 1.08243974e+00, 1.07858726e+00,
1.07484016e+00, 1.07119683e+00, 1.06765606e+00,
1.06421637e+00, 1.06087710e+00, 1.05763739e+00,
1.05449544e+00, 1.05145051e+00, 1.04850133e+00,
1.04564706e+00, 1.04288669e+00, 1.04021982e+00,
1.03764428e+00, 1.03515979e+00, 1.03276563e+00,
1.03046110e+00, 1.02824533e+00, 1.02611737e+00,
1.02407718e+00, 1.02212363e+00, 1.02025601e+00,
1.01847410e+00, 1.01677680e+00, 1.01516379e+00,
1.01363472e+00, 1.01218899e+00, 1.01082595e+00,
1.00954541e+00, 1.00834692e+00, 1.00723016e+00,
1.00619332e+00, 1.00524261e+00, 1.00436638e+00,
1.00357242e+00, 1.00285864e+00, 1.00222494e+00,
1.00167091e+00, 1.00119567e+00, 1.00080040e+00,
1.00048446e+00, 1.00024750e+00, 1.00008955e+00,
1.00001048e+00, 1.00001075e+00, 1.00008956e+00,
1.00024739e+00, 1.00048461e+00, 1.00080045e+00,
1.00119631e+00, 1.00167001e+00, 1.00222414e+00,
1.00285843e+00, 1.00357191e+00, 1.00436647e+00,
1.00524044e+00, 1.00619454e+00, 1.00722884e+00,
1.00834671e+00, 1.00954521e+00, 1.01082532e+00,
1.01218827e+00, 1.01363405e+00, 1.01516288e+00,
1.01677631e+00, 1.01847395e+00, 1.02025403e+00,
1.02212285e+00, 1.02407617e+00, 1.02611652e+00,
1.02824420e+00, 1.03046001e+00, 1.03276496e+00,
1.03515831e+00, 1.03764436e+00, 1.04022059e+00,
1.04288835e+00, 1.04564966e+00, 1.04850418e+00,
1.05145455e+00, 1.05449808e+00, 1.05764151e+00,
1.06087832e+00, 1.06421675e+00, 1.06765622e+00,
1.07119678e+00, 1.07484016e+00, 1.07858747e+00,
1.08243978e+00, 1.08639907e+00, 1.09046707e+00,
1.09464385e+00, 1.09893094e+00, 1.10333166e+00,
1.10784580e+00, 1.11247538e+00, 1.11722302e+00,
1.12208988e+00, 1.12707828e+00, 1.13218948e+00,
1.13742508e+00, 1.14278683e+00, 1.14827796e+00,
1.15389975e+00, 1.15965467e+00, 1.16554407e+00,
1.17157268e+00, 1.17774064e+00, 1.18405153e+00,
1.19050829e+00, 1.19711157e+00, 1.20386739e+00,
1.21078626e+00, 1.21784364e+00, 1.22506490e+00,
1.23244925e+00, 1.23999768e+00, 1.24771390e+00,
1.25560613e+00, 1.26367176e+00, 1.27191600e+00,
1.28034110e+00, 1.28895156e+00, 1.29775501e+00,
1.30675229e+00, 1.31594690e+00, 1.32534484e+00,
1.33495197e+00, 1.34476199e+00, 1.35479270e+00,
1.36504503e+00, 1.37552221e+00, 1.38623094e+00,
1.39717884e+00, 1.40837167e+00, 1.41980294e+00,
1.43148933e+00, 1.44343851e+00, 1.45565145e+00,
1.46813609e+00, 1.48089881e+00, 1.49394864e+00,
1.50729112e+00, 1.52093155e+00, 1.53487604e+00,
1.54914196e+00, 1.56373686e+00, 1.57865894e+00,
1.59392144e+00, 1.60953337e+00, 1.62550528e+00,
1.64184790e+00, 1.65856462e+00, 1.67567248e+00,
1.69318220e+00, 1.71110379e+00, 1.72944932e+00,
1.74823137e+00, 1.76746170e+00, 1.78715285e+00,
1.80731754e+00, 1.82797388e+00, 1.84912923e+00,
1.87080253e+00, 1.89301215e+00, 1.91577034e+00,
1.93909449e+00, 1.96300533e+00, 1.98751647e+00,
2.01265063e+00, 2.03842484e+00, 2.06486280e+00,
2.09198023e+00, 2.11980214e+00, 2.14836366e+00,
2.17767324e+00, 2.20777724e+00, 2.23866942e+00,
2.27040453e+00, 2.30299720e+00, 2.33648709e+00,
2.37089275e+00, 2.40624943e+00, 2.44260837e+00,
2.47999577e+00, 2.51844673e+00, 2.55800224e+00,
2.59870653e+00, 2.64059610e+00, 2.68372572e+00,
2.72815253e+00, 2.77388772e+00, 2.82101920e+00,
2.86958647e+00, 2.91965192e+00, 2.97127496e+00,
3.02450906e+00, 3.07944172e+00, 3.13614110e+00,
3.19466985e+00, 3.25511032e+00, 3.31754917e+00,
3.38208071e+00, 3.44879433e+00, 3.51775774e+00,
3.58910367e+00, 3.66291523e+00, 3.73934686e+00,
3.81850354e+00, 3.90050956e+00, 3.98550877e+00,
4.07365553e+00, 4.16507153e+00, 4.25994051e+00,
4.35843505e+00, 4.46073629e+00, 4.56705238e+00,
4.67759894e+00, 4.79257854e+00, 4.91224585e+00,
5.03685373e+00, 5.16668317e+00, 5.30201735e+00,
5.44318666e+00, 5.59051119e+00, 5.74440639e+00,
5.90517017e+00, 6.07325325e+00, 6.24912158e+00,
6.43325815e+00, 6.62613577e+00, 6.82846496e+00,
7.04065759e+00, 7.26355195e+00, 7.49781116e+00,
7.74415234e+00, 8.00348897e+00, 8.27681715e+00,
8.56496165e+00, 8.86921713e+00, 9.19069767e+00,
9.53087662e+00, 9.89104353e+00, 1.02727244e+01,
1.06778672e+01, 1.11082077e+01, 1.15661555e+01,
1.20538878e+01, 1.25743641e+01, 1.31300086e+01,
1.37241540e+01, 1.43618859e+01, 1.50459005e+01,
1.57809517e+01, 1.65726144e+01, 1.74266981e+01,
1.83500881e+01, 1.93504067e+01, 2.04364838e+01,
2.16180945e+01, 2.29073294e+01, 2.43175531e+01,
2.58646704e+01, 2.75658519e+01, 2.94431019e+01,
3.15206588e+01, 3.38305537e+01, 3.64063882e+01,
3.92903385e+01, 4.25349736e+01, 4.62024177e+01,
5.03692303e+01, 5.51315261e+01, 6.06079334e+01,
6.69473866e+01, 7.43440281e+01, 8.30479655e+01,
9.33771614e+01, 1.05768695e+02, 1.20807161e+02,
1.39316160e+02, 1.62453759e+02, 1.91872994e+02,
2.30095280e+02, 2.80990973e+02, 3.50916358e+02,
4.50628770e+02, 5.99838542e+02, 8.37606689e+02,
1.25105588e+03, 2.06769124e+03, 4.05233265e+03,
1.12559219e+04, 1.01273379e+05, 1.01339517e+05,
1.12547674e+04, 4.05231614e+03, 2.06779417e+03,
1.25114857e+03, 8.37718494e+02, 5.99885295e+02,
4.50639589e+02, 3.50907722e+02, 2.80985545e+02,
2.30076856e+02, 1.91863865e+02, 1.62446368e+02,
1.39323680e+02, 1.20810259e+02, 1.05769149e+02,
9.33790183e+01, 8.30488328e+01, 7.43482299e+01,
6.69506610e+01, 6.06088520e+01, 5.51323789e+01,
5.03691519e+01, 4.62009519e+01, 4.25334553e+01,
3.92888954e+01, 3.64046734e+01, 3.38293614e+01,
3.15201425e+01, 2.94418838e+01, 2.75647804e+01,
2.58634902e+01, 2.43169254e+01, 2.29068548e+01,
2.16174820e+01, 2.04357631e+01, 1.93497241e+01,
1.83496026e+01, 1.74262199e+01, 1.65721321e+01,
1.57804689e+01, 1.50454574e+01, 1.43618634e+01,
1.37249152e+01, 1.31303039e+01, 1.25744447e+01,
1.20540863e+01, 1.15662599e+01, 1.11083237e+01,
1.06777609e+01, 1.02725654e+01, 9.89096570e+00,
9.53098356e+00, 9.19089160e+00, 8.86940718e+00,
8.56516496e+00, 8.27683013e+00, 8.00352002e+00,
7.74416014e+00, 7.49782391e+00, 7.26358505e+00,
7.04068369e+00, 6.82844334e+00, 6.62616388e+00,
6.43323687e+00, 6.24910641e+00, 6.07323018e+00,
5.90512609e+00, 5.74436305e+00, 5.59054792e+00,
5.44322889e+00, 5.30207331e+00, 5.16671529e+00,
5.03689044e+00, 4.91229000e+00, 4.79260205e+00,
4.67762631e+00, 4.56708817e+00, 4.46074456e+00,
4.35843395e+00, 4.25994400e+00, 4.16508670e+00,
4.07366933e+00, 3.98555050e+00, 3.90055633e+00,
3.81855592e+00, 3.73939623e+00, 3.66296119e+00,
3.58912525e+00, 3.51775818e+00, 3.44877174e+00,
3.38207469e+00, 3.31755643e+00, 3.25512893e+00,
3.19469126e+00, 3.13617318e+00, 3.07947611e+00,
3.02453605e+00, 2.97129126e+00, 2.91966190e+