introduced pLongInt for Long Integers

trunk/CPFEM.f90

@@ -482,14 +482,14 @@
    i,&              ! integration point number
    e&               ! element number
  )
- use prec, only: pReal,pInt,reltol_Outer
+ use prec, only: pReal,pInt,pLongInt,reltol_Outer
  use constitutive, only: constitutive_dotState,constitutive_sizeDotState,&
                          constitutive_state_old,constitutive_state_new
  use debug
  
  logical integrateState
 
- integer(8) tick,tock,tickrate,maxticks
+ integer(pLongInt) tick,tock,tickrate,maxticks
  integer(pInt) g,i,e,mySize
  real(pReal), dimension(6) :: Tstar_v
  real(pReal) dt
@@ -659,7 +659,7 @@
  logical failed
  integer(pInt) cp_en, ip, grain
  integer(pInt) iInner,dummy, i,j,k,l,m,n
- integer(8) tick,tock,tickrate,maxticks
+ integer(pLongInt) tick,tock,tickrate,maxticks
  real(pReal) dt, Temperature, det, p_hydro, leapfrog,maxleap
  real(pReal), dimension(6) :: Tstar_v
  real(pReal), dimension(9,9) :: dLp,dTdLp,dRdLp,invdRdLp,eye2

trunk/debug.f90

@@ -9,8 +9,8 @@
  integer(pInt), dimension(nCutback+1) :: debug_cutbackDistribution = 0_pInt
  integer(pInt), dimension(nInner) :: debug_InnerLoopDistribution = 0_pInt
  integer(pInt), dimension(nOuter) :: debug_OuterLoopDistribution = 0_pInt
- integer(8) :: debug_cumLpTicks = 0_pInt
+ integer(pLongInt) :: debug_cumLpTicks = 0_pInt
- integer(8) :: debug_cumDotStateTicks = 0_pInt
+ integer(pLongInt) :: debug_cumDotStateTicks = 0_pInt
  integer(pInt) :: debug_cumLpCalls = 0_pInt
  integer(pInt) :: debug_cumDotStateCalls = 0_pInt
  logical :: debugger = .false.
@@ -28,7 +28,7 @@
  implicit none
 
  integer(pInt) i,integral
- integer(8) tickrate
+ integer(pLongInt) tickrate
 
  write(6,*)
  write(6,*) 'DEBUG Info'
@@ -56,7 +56,7 @@
  write(6,*)
  
  integral = 0_pInt
- write(6,*)	'distribution_InnerLoop :'
+ write(6,*)	'distribution_InnerLoop :'
  do i=1,nInner
    if (debug_InnerLoopDistribution(i) /= 0) then
      integral = integral + i*debug_InnerLoopDistribution(i)

trunk/math.f90

@@ -694,7 +694,9 @@ end subroutine
  real(pReal),dimension(3,3),intent(in)  :: A
  real(pReal) DetA
 
- real(pReal),dimension(3,3) :: math_inv3x3 = 0.0_pReal
+ real(pReal),dimension(3,3) :: math_inv3x3
+ 
+ math_inv3x3 = 0.0_pReal
 
  DetA =   A(1,1) * ( A(2,2) * A(3,3) - A(2,3) * A(3,2) )&
         - A(1,2) * ( A(2,1) * A(3,3) - A(2,3) * A(3,1) )&

trunk/prec.f90

@@ -8,6 +8,8 @@
 !    *** Precision of real and integer variables ***
  integer, parameter :: pReal = selected_real_kind(15,300)  ! 15 significant digits, up to 1e+-300
  integer, parameter :: pInt  = selected_int_kind(9)        ! up to +- 1e9
+ integer, parameter :: pLongInt  = 4                       ! should be 64bit
+
 
  type :: p_vec
      real(pReal), dimension(:), pointer :: p