223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity
DAMASK_EICMD/code/prec_single.f90

83 lines
2.8 KiB
Fortran
Raw Normal View History

changed enconding of all source files to UTF-8 without BOM (signature) Codepage 65001
2011-04-07 12:50:28 +05:30
! Copyright 2011 Max-Planck-Institut für Eisenforschung GmbH
added copyright text to all f90 (free) format files
2011-04-04 19:39:54 +05:30
!
! This file is part of DAMASK,
changed enconding of all source files to UTF-8 without BOM (signature) Codepage 65001
2011-04-07 12:50:28 +05:30
! the Düsseldorf Advanced MAterial Simulation Kit.
added copyright text to all f90 (free) format files
2011-04-04 19:39:54 +05:30
!
! DAMASK is free software: you can redistribute it and/or modify
! it under the terms of the GNU General Public License as published by
! the Free Software Foundation, either version 3 of the License, or
! (at your option) any later version.
!
! DAMASK is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
! GNU General Public License for more details.
!
! You should have received a copy of the GNU General Public License
! along with DAMASK. If not, see <http://www.gnu.org/licenses/>.
!
!##############################################################
moved DAMASK_NaN to prec.f90 respectively prec_single.f90 as it is precision dependent
2011-10-20 22:16:11 +05:30
!* $Id$
added single precision libraries for FFTW First try of implement single precision crystal plasticity, not working yet. polishing text about geometry construction. polishing postResults, still having problems concerning machines without MSC installation
2011-02-25 00:09:57 +05:30
!##############################################################
utf8 encoding output to screen now for all solvers using intel ifort>10 or any gfortran (old versions might not support it, but we don't use them anyway)
2012-04-28 16:16:41 +05:30
#ifdef __INTEL_COMPILER
#if __INTEL_COMPILER<1200
#define LEGACY_COMPILER
#endif
#endif
added error for to IO to substitute stop statement in kdtree2 explicitly defined all functions in as either public or private in the modules to have a quick overview on all functions and parameters that are available
2012-03-06 20:22:48 +05:30
module prec
added single precision libraries for FFTW First try of implement single precision crystal plasticity, not working yet. polishing text about geometry construction. polishing postResults, still having problems concerning machines without MSC installation
2011-02-25 00:09:57 +05:30
!##############################################################
added error for to IO to substitute stop statement in kdtree2 explicitly defined all functions in as either public or private in the modules to have a quick overview on all functions and parameters that are available
2012-03-06 20:22:48 +05:30
implicit none
private
added single precision libraries for FFTW First try of implement single precision crystal plasticity, not working yet. polishing text about geometry construction. polishing postResults, still having problems concerning machines without MSC installation
2011-02-25 00:09:57 +05:30
! *** Precision of real and integer variables ***
added error for to IO to substitute stop statement in kdtree2 explicitly defined all functions in as either public or private in the modules to have a quick overview on all functions and parameters that are available
2012-03-06 20:22:48 +05:30
integer, parameter, public :: pReal = selected_real_kind(6,37) ! 6 significant digits, up to 1e+-37
integer, parameter, public :: pInt = selected_int_kind(9) ! up to +- 1e9
integer, parameter, public :: pLongInt = 4 ! should be 64bit
real(pReal), parameter, public :: tol_math_check = 1.0e-5_pReal
real(pReal), parameter, public :: tol_gravityNodePos = 1.0e-36_pReal
small polishing, exchanged stops in math.f90 with calls to IO_error
2012-02-13 19:38:07 +05:30
added compilation_info.f90, to store compilation date,time and compiler info in each file added #include statement at each init() routine
2012-01-31 20:24:49 +05:30
! NaN is precision dependent
moved DAMASK_NaN to prec.f90 respectively prec_single.f90 as it is precision dependent
2011-10-20 22:16:11 +05:30
! from http://www.hpc.unimelb.edu.au/doc/f90lrm/dfum_035.html
! copy can be found in documentation/Code/Fortran
utf8 encoding output to screen now for all solvers using intel ifort>10 or any gfortran (old versions might not support it, but we don't use them anyway)
2012-04-28 16:16:41 +05:30
#ifdef LEGACY_COMPILER
added error for to IO to substitute stop statement in kdtree2 explicitly defined all functions in as either public or private in the modules to have a quick overview on all functions and parameters that are available
2012-03-06 20:22:48 +05:30
real(pReal), parameter, public :: DAMASK_NaN = Z'Z'7F800001', pReal'
small polishing, exchanged stops in math.f90 with calls to IO_error
2012-02-13 19:38:07 +05:30
#else
added error for to IO to substitute stop statement in kdtree2 explicitly defined all functions in as either public or private in the modules to have a quick overview on all functions and parameters that are available
2012-03-06 20:22:48 +05:30
real(pReal), parameter, public :: DAMASK_NaN = real(Z'7F800001', pReal)
small polishing, exchanged stops in math.f90 with calls to IO_error
2012-02-13 19:38:07 +05:30
#endif
added single precision libraries for FFTW First try of implement single precision crystal plasticity, not working yet. polishing text about geometry construction. polishing postResults, still having problems concerning machines without MSC installation
2011-02-25 00:09:57 +05:30
added error for to IO to substitute stop statement in kdtree2 explicitly defined all functions in as either public or private in the modules to have a quick overview on all functions and parameters that are available
2012-03-06 20:22:48 +05:30
type, public :: p_vec
real(pReal), dimension(:), pointer :: p
end type p_vec
public :: prec_init
contains
added single precision libraries for FFTW First try of implement single precision crystal plasticity, not working yet. polishing text about geometry construction. polishing postResults, still having problems concerning machines without MSC installation
2011-02-25 00:09:57 +05:30
subroutine prec_init
added error for to IO to substitute stop statement in kdtree2 explicitly defined all functions in as either public or private in the modules to have a quick overview on all functions and parameters that are available
2012-03-06 20:22:48 +05:30
use, intrinsic :: iso_fortran_env ! to get compiler_version and compiler_options (at least for gfortran 4.6 at the moment)
implicit none
added compilation_info.f90, to store compilation date,time and compiler info in each file added #include statement at each init() routine
2012-01-31 20:24:49 +05:30
openmp parallelization working again (at least for j2 and nonlocal constitutive model). In order to keep it like that, please follow these simple rules: DON'T use implicit array subscripts: example: real, dimension(3,3) :: A,B A(:,2) = B(:,1) <--- DON'T USE A(1:3,2) = B(1:3,1) <--- BETTER USE In many cases the use of explicit array subscripts is inevitable for parallelization. Additionally, it is an easy means to prevent memory leaks. Enclose all write statements with the following: !$OMP CRITICAL (write2out) <your write statement> !$OMP END CRITICAL (write2out) Whenever you change something in the code and are not sure if it affects parallelization and leads to nonconforming behavior, please ask me and/or Franz to check this.
2011-03-17 16:16:17 +05:30
!$OMP CRITICAL (write2out)
utf8 encoding output to screen now for all solvers using intel ifort>10 or any gfortran (old versions might not support it, but we don't use them anyway)
2012-04-28 16:16:41 +05:30
#ifndef LEGACY_COMPILER
open (6, encoding='UTF-8')
#endif
added error for to IO to substitute stop statement in kdtree2 explicitly defined all functions in as either public or private in the modules to have a quick overview on all functions and parameters that are available
2012-03-06 20:22:48 +05:30
write(6,*)
write(6,*) '<<<+- prec_single init -+>>>'
write(6,*) '$Id$'
made code standard conform to Fortran 2008 (ignoring warning concerning comments beyond character 132). Basically, changing "x" format specifier to "Nx" ("1x") plus removing $ format specifier added compiler switches for gfortran and ifort to check for standard conformity old gnu compilers <4.4 are not longer supported because they don't provide the c binding for fftw
2012-02-01 00:48:55 +05:30
#include "compilation_info.f90"
added error for to IO to substitute stop statement in kdtree2 explicitly defined all functions in as either public or private in the modules to have a quick overview on all functions and parameters that are available
2012-03-06 20:22:48 +05:30
write(6,'(a,i3)') ' Bytes for pReal: ',pReal
write(6,'(a,i3)') ' Bytes for pInt: ',pInt
write(6,'(a,i3)') ' Bytes for pLongInt: ',pLongInt
write(6,'(a,e10.3)') ' NaN: ',DAMASK_NAN
write(6,'(a,l3)') ' NaN /= NaN: ',DAMASK_NaN/=DAMASK_NaN
if (DAMASK_NaN == DAMASK_NaN) call quit(9000)
write(6,*)
openmp parallelization working again (at least for j2 and nonlocal constitutive model). In order to keep it like that, please follow these simple rules: DON'T use implicit array subscripts: example: real, dimension(3,3) :: A,B A(:,2) = B(:,1) <--- DON'T USE A(1:3,2) = B(1:3,1) <--- BETTER USE In many cases the use of explicit array subscripts is inevitable for parallelization. Additionally, it is an easy means to prevent memory leaks. Enclose all write statements with the following: !$OMP CRITICAL (write2out) <your write statement> !$OMP END CRITICAL (write2out) Whenever you change something in the code and are not sure if it affects parallelization and leads to nonconforming behavior, please ask me and/or Franz to check this.
2011-03-17 16:16:17 +05:30
!$OMP END CRITICAL (write2out)
added single precision libraries for FFTW First try of implement single precision crystal plasticity, not working yet. polishing text about geometry construction. polishing postResults, still having problems concerning machines without MSC installation
2011-02-25 00:09:57 +05:30
added error for to IO to substitute stop statement in kdtree2 explicitly defined all functions in as either public or private in the modules to have a quick overview on all functions and parameters that are available
2012-03-06 20:22:48 +05:30
end subroutine prec_init
added single precision libraries for FFTW First try of implement single precision crystal plasticity, not working yet. polishing text about geometry construction. polishing postResults, still having problems concerning machines without MSC installation
2011-02-25 00:09:57 +05:30
added error for to IO to substitute stop statement in kdtree2 explicitly defined all functions in as either public or private in the modules to have a quick overview on all functions and parameters that are available
2012-03-06 20:22:48 +05:30
end module prec