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.f90

104 lines
5.1 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,
merged precision info into the file prec.f90 and removed prec_single.f90 precision information is now parsed as a makro, defined by the compiler, eg. -DFLOAT=8 for 8 byte long floats or -- in case no values is given -- for the FEM wrapper by the interface files. Please don't forget to run setup_code.py !!!!!!!!!!! abaqus_v6.env is adjusted to specify the standard precision (float=8, int=4) and got some extra compilation options, some deprecated warning supressors are removed Optional Marc modification will follow that allow automatic selection of suitable precision
2012-08-28 21:38:17 +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/>.
!
changed comments in prec and IO to be read by doxygen, fine tuned output of doxygen and added precompilation for documentation
2012-08-09 16:31:53 +05:30
!--------------------------------------------------------------------------------------------------
! $Id$
!--------------------------------------------------------------------------------------------------
!> @author Franz Roters, Max-Planck-Institut für Eisenforschung GmbH
!> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH
!> @author Christoph Kords, Max-Planck-Institut für Eisenforschung GmbH
!> @author Martin Diehl, Max-Planck-Institut für Eisenforschung GmbH
merged precision info into the file prec.f90 and removed prec_single.f90 precision information is now parsed as a makro, defined by the compiler, eg. -DFLOAT=8 for 8 byte long floats or -- in case no values is given -- for the FEM wrapper by the interface files. Please don't forget to run setup_code.py !!!!!!!!!!! abaqus_v6.env is adjusted to specify the standard precision (float=8, int=4) and got some extra compilation options, some deprecated warning supressors are removed Optional Marc modification will follow that allow automatic selection of suitable precision
2012-08-28 21:38:17 +05:30
!> @brief setting precision for real and int type depending on makros "FLOAT" and "INT"
changed comments in prec and IO to be read by doxygen, fine tuned output of doxygen and added precompilation for documentation
2012-08-09 16:31:53 +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
fixed bug causing trouble with ifort 11 (abaqus)
2012-08-29 00:40:54 +05:30
#if __INTEL_COMPILER<1200
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
#define LEGACY_COMPILER
#endif
#endif
changed comments in prec and IO to be read by doxygen, fine tuned output of doxygen and added precompilation for documentation
2012-08-09 16:31:53 +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
module prec
changed comments in prec and IO to be read by doxygen, fine tuned output of doxygen and added precompilation for documentation
2012-08-09 16:31:53 +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
changed comments in prec and IO to be read by doxygen, fine tuned output of doxygen and added precompilation for documentation
2012-08-09 16:31:53 +05:30
merged precision info into the file prec.f90 and removed prec_single.f90 precision information is now parsed as a makro, defined by the compiler, eg. -DFLOAT=8 for 8 byte long floats or -- in case no values is given -- for the FEM wrapper by the interface files. Please don't forget to run setup_code.py !!!!!!!!!!! abaqus_v6.env is adjusted to specify the standard precision (float=8, int=4) and got some extra compilation options, some deprecated warning supressors are removed Optional Marc modification will follow that allow automatic selection of suitable precision
2012-08-28 21:38:17 +05:30
#if (FLOAT==4)
integer, parameter, public :: pReal = 4 !< floating point single precition (was selected_real_kind(6,37), number with 6 significant digits, up to 1e+-37)
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
merged precision info into the file prec.f90 and removed prec_single.f90 precision information is now parsed as a makro, defined by the compiler, eg. -DFLOAT=8 for 8 byte long floats or -- in case no values is given -- for the FEM wrapper by the interface files. Please don't forget to run setup_code.py !!!!!!!!!!! abaqus_v6.env is adjusted to specify the standard precision (float=8, int=4) and got some extra compilation options, some deprecated warning supressors are removed Optional Marc modification will follow that allow automatic selection of suitable precision
2012-08-28 21:38:17 +05:30
real(pReal), parameter, public :: DAMASK_NaN = Z'7F800001' !< quiet NaN for single precision (from http://www.hpc.unimelb.edu.au/doc/f90lrm/dfum_035.html, copy can be found in documentation/Code/Fortran)
small polishing, exchanged stops in math.f90 with calls to IO_error
2012-02-13 19:38:07 +05:30
#else
merged precision info into the file prec.f90 and removed prec_single.f90 precision information is now parsed as a makro, defined by the compiler, eg. -DFLOAT=8 for 8 byte long floats or -- in case no values is given -- for the FEM wrapper by the interface files. Please don't forget to run setup_code.py !!!!!!!!!!! abaqus_v6.env is adjusted to specify the standard precision (float=8, int=4) and got some extra compilation options, some deprecated warning supressors are removed Optional Marc modification will follow that allow automatic selection of suitable precision
2012-08-28 21:38:17 +05:30
real(pReal), parameter, public :: DAMASK_NaN = real(Z'7F800001', pReal) !< quiet NaN for single precision (from http://www.hpc.unimelb.edu.au/doc/f90lrm/dfum_035.html, copy can be found in documentation/Code/Fortran)
#endif
#elif (FLOAT==8)
fixed bug causing trouble with ifort 11 (abaqus)
2012-08-29 00:40:54 +05:30
integer, parameter, public :: pReal = 8 !< floating point double precision (was selected_real_kind(15,300), number with 15 significant digits, up to 1e+-300)
merged precision info into the file prec.f90 and removed prec_single.f90 precision information is now parsed as a makro, defined by the compiler, eg. -DFLOAT=8 for 8 byte long floats or -- in case no values is given -- for the FEM wrapper by the interface files. Please don't forget to run setup_code.py !!!!!!!!!!! abaqus_v6.env is adjusted to specify the standard precision (float=8, int=4) and got some extra compilation options, some deprecated warning supressors are removed Optional Marc modification will follow that allow automatic selection of suitable precision
2012-08-28 21:38:17 +05:30
#ifdef LEGACY_COMPILER
real(pReal), parameter, public :: DAMASK_NaN = Z'7FF8000000000000' !< quiet NaN for double precision (from http://www.hpc.unimelb.edu.au/doc/f90lrm/dfum_035.html, copy can be found in documentation/Code/Fortran)
#else
real(pReal), parameter, public :: DAMASK_NaN = real(Z'7FF8000000000000', pReal) !< quiet NaN for double precision (from http://www.hpc.unimelb.edu.au/doc/f90lrm/dfum_035.html, copy can be found in documentation/Code/Fortran)
#endif
small polishing, exchanged stops in math.f90 with calls to IO_error
2012-02-13 19:38:07 +05:30
#endif
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
merged precision info into the file prec.f90 and removed prec_single.f90 precision information is now parsed as a makro, defined by the compiler, eg. -DFLOAT=8 for 8 byte long floats or -- in case no values is given -- for the FEM wrapper by the interface files. Please don't forget to run setup_code.py !!!!!!!!!!! abaqus_v6.env is adjusted to specify the standard precision (float=8, int=4) and got some extra compilation options, some deprecated warning supressors are removed Optional Marc modification will follow that allow automatic selection of suitable precision
2012-08-28 21:38:17 +05:30
#if (INT==4)
integer, parameter, public :: pInt = 4 !< integer representation 32 bit (was selected_int_kind(9), number with at least up to +- 1e9)
#elif (INT==8)
integer, parameter, public :: pInt = 8 !< integer representation 64 bit (was selected_int_kind(12), number with at least up to +- 1e12)
#endif
integer, parameter, public :: pLongInt = 8 !< integer representation 64 bit (was selected_int_kind(12), number with at least up to +- 1e12)
real(pReal), parameter, public :: tol_math_check = 1.0e-8_pReal
real(pReal), parameter, public :: tol_gravityNodePos = 1.0e-100_pReal
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
added version information to all files do NOT edit text like this: $Id: constitutive_phenopowerlaw.f90 406 2009-08-31 14:13:10Z MPIE\f.roters $
2009-08-31 20:39:15 +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
public :: prec_init
contains
merged precision info into the file prec.f90 and removed prec_single.f90 precision information is now parsed as a makro, defined by the compiler, eg. -DFLOAT=8 for 8 byte long floats or -- in case no values is given -- for the FEM wrapper by the interface files. Please don't forget to run setup_code.py !!!!!!!!!!! abaqus_v6.env is adjusted to specify the standard precision (float=8, int=4) and got some extra compilation options, some deprecated warning supressors are removed Optional Marc modification will follow that allow automatic selection of suitable precision
2012-08-28 21:38:17 +05:30
changed comments in prec and IO to be read by doxygen, fine tuned output of doxygen and added precompilation for documentation
2012-08-09 16:31:53 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief reporting precision and checking if DAMASK_NaN is set correctly
!--------------------------------------------------------------------------------------------------
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
subroutine prec_init
changed comments in prec and IO to be read by doxygen, fine tuned output of doxygen and added precompilation for documentation
2012-08-09 16:31:53 +05:30
use, intrinsic :: iso_fortran_env ! to get compiler_version and compiler_options (at least for gfortran 4.6 at the moment)
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
added "-assume byterecl" in Makefile for ifortran to specify with "recl=" the size of the file in bytes instead of 4byte chunk when opening files. (same as in gfortran per default). no important changes in other files
2012-03-20 17:56:21 +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 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
corrected NaN in prec_single.f90 (was already silent but with typo), polished output and added switch in order to check range on all files except of prec.f90 (in which NaN is defined)
2012-05-11 18:16:17 +05:30
write(6,'(a,e10.3)') ' NaN: ', DAMASK_NaN
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,l3)') ' NaN /= NaN: ',DAMASK_NaN/=DAMASK_NaN
write(6,*)
merged precision info into the file prec.f90 and removed prec_single.f90 precision information is now parsed as a makro, defined by the compiler, eg. -DFLOAT=8 for 8 byte long floats or -- in case no values is given -- for the FEM wrapper by the interface files. Please don't forget to run setup_code.py !!!!!!!!!!! abaqus_v6.env is adjusted to specify the standard precision (float=8, int=4) and got some extra compilation options, some deprecated warning supressors are removed Optional Marc modification will follow that allow automatic selection of suitable precision
2012-08-28 21:38:17 +05:30
if (DAMASK_NaN == DAMASK_NaN) call quit(9000)
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 version information to all files do NOT edit text like this: $Id: constitutive_phenopowerlaw.f90 406 2009-08-31 14:13:10Z MPIE\f.roters $
2009-08-31 20:39:15 +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 version information to all files do NOT edit text like this: $Id: constitutive_phenopowerlaw.f90 406 2009-08-31 14:13:10Z MPIE\f.roters $
2009-08-31 20:39:15 +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