diff options
Diffstat (limited to 'src/qfunc.x')
| -rw-r--r-- | src/qfunc.x | 29 |
1 files changed, 29 insertions, 0 deletions
diff --git a/src/qfunc.x b/src/qfunc.x new file mode 100644 index 0000000..d26919d --- /dev/null +++ b/src/qfunc.x @@ -0,0 +1,29 @@ + qf(i,j,k) = amp*(1.00000000000000d0+c*cos(phigrd(k))**2)* + $ (exp(-((etagrd(i)-eta0)**2/sigma**2))+exp(-((etagrd(i)+eta0)**2 + $ /sigma**2)))*sin(qgrd(j))**n + + qfetaeta(i,j,k) = amp*(1.00000000000000d0+c*cos(phigrd(k)) + $ **2)*((4.0000000000000d0*(etagrd(i)-eta0)**2*exp(-((etagrd(i)-eta0) + $ **2/sigma**2)))/sigma**4-(2.00000000000000d0*exp(-((etagrd(i) + $ -eta0)**2/sigma**2)))/sigma**2+(4.0000000000000d0*(etagrd(i)+ + $ eta0)**2*exp(-((etagrd(i)+eta0)**2/sigma**2)))/sigma**4-(2.00 + $ 000000000000d0*exp(-((etagrd(i)+eta0)**2/sigma**2)))/sigma**2 + $ )*sin(qgrd(j))**n + + qfqq(i,j,k) = amp*(-1.00000000000000d0+n)*n*(1.00000000 + $ 000000d0+c*cos(phigrd(k))**2)*cos(qgrd(j))**2*(exp(-((etagrd(i)-eta0)**2/ + $ sigma**2))+exp(-((etagrd(i)+eta0)**2/sigma**2)))*sin(qgrd(j))**(-2. + $ 00000000000000d0+n)-amp*n*(1.00000000000000d0+c*cos(phigrd(k)) + $ **2)*(exp(-((etagrd(i)-eta0)**2/sigma**2))+exp(-((etagrd(i)+eta0)**2 + $ /sigma**2)))*sin(qgrd(j))**n + + qfphi(i,j,k) = -2.00000000000000d0*amp*c*cos(phigrd(k))*(exp(- + $ ((etagrd(i)-eta0)**2/sigma**2))+exp(-((etagrd(i)+eta0)**2/sigma**2)) + $ )*sin(phigrd(k))*sin(qgrd(j))**n + + qfphiphi(i,j,k) = -2.00000000000000d0*amp*c*cos(phigrd(k))**2* + $ (exp(-((etagrd(i)-eta0)**2/sigma**2))+exp(-((etagrd(i)+eta0)**2/ + $ sigma**2)))*sin(qgrd(j))**n+2.00000000000000d0*amp*c*(exp(-(( + $ etagrd(i)-eta0)**2/sigma**2))+exp(-((etagrd(i)+eta0)**2/sigma**2)))* + $ sin(phigrd(k))**2*sin(qgrd(j))**n + |