Initial commit

1 files changed, 114 insertions, 0 deletions

diff --git a/math/ieee/chap1/const.f b/math/ieee/chap1/const.f
new file mode 100644
index 00000000..54eb4e41
--- /dev/null
+++ b/math/ieee/chap1/const.f
@@ -0,0 +1,114 @@
+c
+c-----------------------------------------------------------------------
+c subroutine: const
+c   computes the multipliers for the various modules
+c-----------------------------------------------------------------------
+c
+      subroutine const(co3,co8,co16,co9,cdc,cdd)
+      double precision dtheta,dtwopi,dsq32,dsq2
+      double precision dcos1,dcos2,dcos3,dcos4
+      double precision dsin1,dsin2,dsin3,dsin4
+      dimension co3(3),co8(8),co16(18),co9(11),cdc(9),cdd(6)
+      dtwopi=8.0d0*datan(1.0d0)
+      dsq32=dsqrt(0.75d0)
+      dsq2=dsqrt(0.5d0)
+c
+c   multipliers for the three point module
+c
+      co3(1)=1.0
+      co3(2)=-1.5
+      co3(3)=-dsq32
+c
+c   multipliers for the five point module
+c
+      dtheta=dtwopi/5.0d0
+      dcos1=dcos(dtheta)
+      dcos2=dcos(2.0d0*dtheta)
+      dsin1=dsin(dtheta)
+      dsin2=dsin(2.0d0*dtheta)
+      cdd(1)=1.0
+      cdd(2)=-1.25
+      cdd(3)=-dsin1-dsin2
+      cdd(4)=0.5*(dcos1-dcos2)
+      cdd(5)=dsin1-dsin2
+      cdd(6)=dsin2
+c
+c
+c   multipliers for the seven point module
+c
+      dtheta=dtwopi/7.0d0
+      dcos1=dcos(dtheta)
+      dcos2=dcos(2.0d0*dtheta)
+      dcos3=dcos(3.0d0*dtheta)
+      dsin1=dsin(dtheta)
+      dsin2=dsin(2.0d0*dtheta)
+      dsin3=dsin(3.0d0*dtheta)
+      cdc(1)=1.0
+      cdc(2)=-7.0d0/6.0d0
+      cdc(3)=-(dsin1+dsin2-dsin3)/3.0d0
+      cdc(4)=(dcos1+dcos2-2.0d0*dcos3)/3.0d0
+      cdc(5)=(2.0d0*dcos1-dcos2-dcos3)/3.0d0
+      cdc(6)=-(2.0d0*dsin1-dsin2+dsin3)/3.0d0
+      cdc(7)=-(dsin1+dsin2+2.0d0*dsin3)/3.0d0
+      cdc(8)=(dcos1-2.0d0*dcos2+dcos3)/3.0d0
+      cdc(9)=-(dsin1-2.0d0*dsin2-dsin3)/3.0d0
+c
+c   multipliers for the eight point module
+c
+      co8(1)=1.0
+      co8(2)=1.0
+      co8(3)=1.0
+      co8(4)=-1.0
+      co8(5)=1.0
+      co8(6)=-dsq2
+      co8(7)=-1.0
+      co8(8)=dsq2
+c
+c   multipliers for the nine point module
+c
+      dtheta=dtwopi/9.0d0
+      dcos1=dcos(dtheta)
+      dcos2=dcos(2.0d0*dtheta)
+      dcos4=dcos(4.0d0*dtheta)
+      dsin1=dsin(dtheta)
+      dsin2=dsin(2.0d0*dtheta)
+      dsin4=dsin(4.0d0*dtheta)
+      co9(1)=1.0
+      co9(2)=-1.5
+      co9(3)=-dsq32
+      co9(4)=0.5
+      co9(5)=(2.0d0*dcos1-dcos2-dcos4)/3.0d0
+      co9(6)=(dcos1-2.0d0*dcos2+dcos4)/3.0d0
+      co9(7)=(dcos1+dcos2-2.0d0*dcos4)/3.0d0
+      co9(8)=-(2.0d0*dsin1+dsin2-dsin4)/3.0d0
+      co9(9)=-(dsin1+2.0d0*dsin2+dsin4)/3.0d0
+      co9(10)=-(dsin1-dsin2-2.0d0*dsin4)/3.0d0
+      co9(11)=-dsq32
+c
+c   multipliers for the sixteen point module
+c
+      dtheta=dtwopi/16.0d0
+      dcos1=dcos(dtheta)
+      dcos3=dcos(3.0d0*dtheta)
+      dsin1=dsin(dtheta)
+      dsin3=dsin(3.0d0*dtheta)
+      co16(1)=1.0
+      co16(2)=1.0
+      co16(3)=1.0
+      co16(4)=-1.0
+      co16(5)=1.0
+      co16(6)=-dsq2
+      co16(7)=-1.0
+      co16(8)=dsq2
+      co16(9)=1.0
+      co16(10)=-(dsin1-dsin3)
+      co16(11)=-dsq2
+      co16(12)=-co16(10)
+      co16(13)=-1.0
+      co16(14)=-(dsin1+dsin3)
+      co16(15)=dsq2
+      co16(16)=-co16(14)
+      co16(17)=-dsin3
+      co16(18)=dcos3
+      return
+      end