spherical_deriv.c checking numerical derivatives rval[0]=0.250000 rval[1]=0.500000 rval[2]=0.750000 rval[3]=1.000000 rval[4]=1.250000 rval[5]=1.500000 rval[6]=1.750000 tval[0]=-3.141593 tval[1]=-2.748894 tval[2]=-2.356194 tval[3]=-1.963495 tval[4]=-1.570796 tval[5]=-1.178097 tval[6]=-0.785398 tval[7]=-0.392699 tval[8]=0.000000 tval[9]=0.392699 tval[10]=0.785398 tval[11]=1.178097 tval[12]=1.570796 tval[13]=1.963495 tval[14]=2.356194 tval[15]=2.748894 pval[0]=-1.570796 pval[1]=-1.308997 pval[2]=-1.047198 pval[3]=-0.785398 pval[4]=-0.523599 pval[5]=-0.261799 pval[6]=-0.000000 pval[7]=0.261799 pval[8]=0.523599 pval[9]=0.785398 pval[10]=1.047198 pval[11]=1.308997 rtp[0][0][0]=0.094804 rtp[0][0][11]=0.345683 rtp[0][15][0]=0.365189 rtp[0][15][11]=1.331585 rtp[6][0][0]=0.004871 rtp[6][0][11]=0.017762 rtp[6][15][0]=0.018765 rtp[6][15][11]=0.068421 at 0,0,0 r=0.250000, t=-3.141593, p=-1.570796 analytic dudr=-0.291705, dudt=0.047665, dudp=0.089595 computed dudr=-0.307615, dudt=-895.069686, dudp=6.462242 at 0,0,11 r=0.250000, t=-3.141593, p=1.308997 analytic dudr=-1.063639, dudt=0.173801, dudp=-0.555893 computed dudr=-1.121652, dudt=-3263.680274, dudp=8.564273 at 0,15,0 r=0.250000, t=2.748894, p=-1.570796 analytic dudr=-1.123660, dudt=-0.353634, dudp=0.345123 computed dudr=-1.184947, dudt=-2310.652724, dudp=24.892853 at 0,15,11 r=0.250000, t=2.748894, p=1.308997 analytic dudr=-4.097186, dudt=-1.289450, dudp=-2.141325 computed dudr=-4.320656, dudt=-8425.301212, dudp=32.989972 at 6,0,0 r=1.750000, t=-3.141593, p=-1.570796 analytic dudr=-0.005391, dudt=0.002449, dudp=0.004604 computed dudr=-0.004336, dudt=-45.991723, dudp=0.332052 at 6,0,11 r=1.750000, t=-3.141593, p=1.308997 analytic dudr=-0.019658, dudt=0.008930, dudp=-0.028564 computed dudr=-0.015811, dudt=-167.698986, dudp=0.440061 at 6,15,0 r=1.750000, t=2.748894, p=-1.570796 analytic dudr=-0.020767, dudt=-0.018171, dudp=0.017734 computed dudr=-0.016704, dudt=-118.729191, dudp=1.279079 at 6,15,11 r=1.750000, t=2.748894, p=1.308997 analytic dudr=-0.075723, dudt=-0.066256, dudp=-0.110029 computed dudr=-0.060906, dudt=-432.920616, dudp=1.695137 rcoef[0]=0.666667 rcoef[1]=-4.800000 rcoef[2]=15.000000 rcoef[3]=-26.666667 rcoef[4]=30.000000 rcoef[5]=-24.000000 rcoef[6]=9.800000 tcoef[0]=-0.169765 tcoef[1]=2.728370 tcoef[2]=-20.567716 tcoef[3]=96.553999 tcoef[4]=-315.994905 tcoef[5]=764.707671 tcoef[6]=-1416.125316 tcoef[7]=2048.324118 tcoef[8]=-2340.941849 tcoef[9]=2124.187974 tcoef[10]=-1529.415341 tcoef[11]=868.985989 tcoef[12]=-386.215995 tcoef[13]=133.690152 tcoef[14]=-38.197186 tcoef[15]=8.449801 pcoef[0]=-0.347247 pcoef[1]=4.201690 pcoef[2]=-23.342725 pcoef[3]=78.781697 pcoef[4]=-180.072450 pcoef[5]=294.118335 pcoef[6]=-352.942002 pcoef[7]=315.126787 pcoef[8]=-210.084525 pcoef[9]=105.042262 pcoef[10]=-42.016905 pcoef[11]=11.535082 /nmm=5, mb=2 at 3,4,3 r=1.000000, t=-1.570796, p=-0.785398 analytic dudr=-0.170256, dudt=0.076721, dudp=0.130715 computed dudr=-0.165901, dudt=0.076087, dudp=0.130263 rcoef[0]=0.333333 rcoef[1]=-2.666667 rcoef[2]=-0.000000 rcoef[3]=2.666667 rcoef[4]=-0.333333 tcoef[0]=0.212207 tcoef[1]=-1.697653 tcoef[2]=-0.000000 tcoef[3]=1.697653 tcoef[4]=-0.212207 pcoef[0]=0.318310 pcoef[1]=-2.546479 pcoef[2]=0.000000 pcoef[3]=2.546479 pcoef[4]=-0.318310 /nmm=7, mb=3 at 3,4,3 r=1.000000, t=-1.570796, p=-0.785398 analytic dudr=-0.170256, dudt=0.076721, dudp=0.130715 computed dudr=-0.170370, dudt=0.076932, dudp=0.130279 rcoef[0]=-0.066667 rcoef[1]=0.600000 rcoef[2]=-3.000000 rcoef[3]=-0.000000 rcoef[4]=3.000000 rcoef[5]=-0.600000 rcoef[6]=0.066667 tcoef[0]=-0.042441 tcoef[1]=0.381972 tcoef[2]=-1.909859 tcoef[3]=0.000000 tcoef[4]=1.909859 tcoef[5]=-0.381972 tcoef[6]=0.042441 pcoef[0]=-0.063662 pcoef[1]=0.572958 pcoef[2]=-2.864789 pcoef[3]=0.000000 pcoef[4]=2.864789 pcoef[5]=-0.572958 pcoef[6]=0.063662 end spherical_deriv check