Commit 14346a90 authored by Matthias Redies's avatar Matthias Redies

redo stuff for .f-files

parent 0ba38f36
......@@ -83,7 +83,7 @@ C ---> save and data statements
c
save one,zero,half,eps,eps1
data zero/0.0/,one/1.0/,half/0.5/,
+ eps/1.0*10.0**-8/,eps1/1.0*10.0**-5/
+ eps/1e-8/,eps1/1e-5/
c
c-----------------------------------------------------------------------
if (kpri .ge. 3) then
......
......@@ -96,7 +96,7 @@
IF( all( matmul(rot(:,:,i),rot(:,:,j))
& .eq.reshape((/1,0,0,0,1,0,0,0,1/),(/3,3/)))
& .and.all(modulo(matmul(rot(:,:,i),rtau(:,j))+rtau(:,i),1.0)
& .lt.10.0**-10) )THEN
& .lt.1e-10) )THEN
IF(invtab(i).ne.0) STOP 'kptgen: inverse operation
& & already defined.'
invtab(i) = j
......
......@@ -109,7 +109,7 @@ C ---> save and data statements
c
save one,zero,half,eps,eps1,iplus,iminus
data zero/0.0/,one/1.0/,half/0.5/,
+ eps/1.0*10.0**-8/,eps1/1.0*10.0**-5/,iplus/1/,iminus/-1/
+ eps/1e-8/,eps1/1e-5/,iplus/1/,iminus/-1/
c
c-----------------------------------------------------------------------
c
......
......@@ -43,8 +43,8 @@ C Fred Hutson
C
C*************************************************************************
C
data eps/10.0**-10/
data eps1/10.0**-5/
data eps/1e-10/
data eps1/1e-5/
pi = 4.0*atan(1.0)
C
C CONSTRUCT THE FIRST TETRAHEDRON
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment