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

1472 lines
45 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/>.
!
!##############################################################
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
!* $Id$
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
!##############################################################
MODULE IO
!##############################################################
removed implicit type castings
2012-02-10 16:54:53 +05:30
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
CONTAINS
!---------------------------
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
! function IO_abaqus_assembleInputFile
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
! subroutine IO_open_file(unit,relPath)
! subroutine IO_open_inputFile(unit, model)
! subroutine IO_open_logFile(unit)
change of expon in fixedFloatNoE to type int, general polishing
2007-03-21 22:34:10 +05:30
! function IO_hybridIA(Nast,ODFfileName)
! private function hybridIA_reps(dV_V,steps,C)
small polishing
2007-04-10 16:51:34 +05:30
! function IO_stringPos(line,maxN)
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
! function IO_stringValue(line,positions,pos)
! function IO_floatValue(line,positions,pos)
! function IO_intValue(line,positions,pos)
small polishing
2007-04-10 16:51:34 +05:30
! function IO_fixedStringValue(line,ends,pos)
! function IO_fixedFloatValue(line,ends,pos)
! function IO_fixedFloatNoEValue(line,ends,pos)
! function IO_fixedIntValue(line,ends,pos)
corrected misspelling of function 'IO_continousIntValues'
2009-04-02 18:32:25 +05:30
! function IO_continousIntValues(unit,maxN)
small polishing
2007-04-10 16:51:34 +05:30
! function IO_lc(line)
! subroutine IO_lcInplace(line)
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
! subroutine IO_error(ID)
added subroutine to throw warnings (instead of terminal errors)
2009-03-31 14:51:57 +05:30
! subroutine IO_warning(ID)
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
!---------------------------
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
!********************************************************************
! output version number
!********************************************************************
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
subroutine IO_init ()
removed implicit type castings
2012-02-10 16:54:53 +05:30
use, intrinsic :: iso_fortran_env
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
!$OMP CRITICAL (write2out)
write(6,*)
write(6,*) '<<<+- IO 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"
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
call flush(6)
!$OMP END CRITICAL (write2out)
convergence of state in crystallite is now tested as follows: (state < relevant state) or (residuum < relative tolerance * state) since the relevant value for the state variables depend on their nature and can vary by large scales (e.g. volume fraction: 1e-10, dislocation density: 1e5) it is not possible to set a unique value. instead the constitutive law has to decide what is relevant. therefore, all constitutive laws now read in parameters from the material.config that determine the values for relevantState [@luc: in dislobased law relevant State is for the moment generally set to 1e-200, so no additional parameters necessary in material.config. if you also want this feature, we can still implement it, no big deal] - added sanity checks in constitutive_nonlocal.f90 - corrected coordinate transformation for backstress calculation in constitutive_nonlocal.f90 - corrected equations for evolution of dipole dislocation densities (athermal annihilation and formation by glide)
2009-09-18 21:07:14 +05:30
endsubroutine
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
!********************************************************************
! AP: 12.07.10
! create a new input file for abaqus simulations
! by removing all comment lines and including "include"s
!********************************************************************
recursive function IO_abaqus_assembleInputFile(unit1,unit2) result(createSuccess)
use prec
renamed interface from "mpie_interface" to "DAMASK_interface"
2011-05-11 22:31:03 +05:30
use DAMASK_interface
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
implicit none
character(len=300) line,fname
integer(pInt), intent(in) :: unit1, unit2
logical createSuccess,fexist
restarting now works with Abaqus (standard, cannot test explicit due to lack of license) you have to specify the job you are restarting from in the job description (cae), if you prepare your input file by hand this is the first line after *Heading example: if the first job was using Oldjob.inp the first entry in the job description needs to be Oldjob (without the .inp) as for Marc restart works only from last converged increment, i.e. ther restart writing should be specified like this: *retsart, write, frequency=1, overlay Overlay is not essential but saves a lot of disk space and as stated before you can only restart from the last converged increment anyway
2011-07-18 14:45:20 +05:30
integer(pInt), parameter :: maxNchunks = 6
integer(pInt), dimension(1+2*maxNchunks) :: positions
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
do
read(unit2,'(A300)',END=220) line
restarting now works with Abaqus (standard, cannot test explicit due to lack of license) you have to specify the job you are restarting from in the job description (cae), if you prepare your input file by hand this is the first line after *Heading example: if the first job was using Oldjob.inp the first entry in the job description needs to be Oldjob (without the .inp) as for Marc restart works only from last converged increment, i.e. ther restart writing should be specified like this: *retsart, write, frequency=1, overlay Overlay is not essential but saves a lot of disk space and as stated before you can only restart from the last converged increment anyway
2011-07-18 14:45:20 +05:30
! line = IO_lc(trim(line))
! do not change the whole line to lower case, file names in Linux are case sensitive!
positions = IO_stringPos(line,maxNchunks)
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
! call IO_lcInPlace(line)
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
if (IO_lc(IO_StringValue(line,positions,1_pInt))=='*include') then
fname = trim(getSolverWorkingDirectoryName())//trim(line(9_pInt+scan(line(9_pInt:),'='):))
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
inquire(file=fname, exist=fexist)
if (.not.(fexist)) then
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)
write(6,*)'ERROR: file does not exist error in IO_abaqus_assembleInputFile'
write(6,*)'filename: ', trim(fname)
!$OMP END CRITICAL (write2out)
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
createSuccess = .false.
return
endif
open(unit2+1,err=200,status='old',file=fname)
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
if (IO_abaqus_assembleInputFile(unit1,unit2+1_pInt)) then
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
createSuccess=.true.
close(unit2+1)
else
createSuccess=.false.
return
endif
else if (line(1:2) /= '**') then
write(unit1,'(A)') trim(line)
endif
enddo
220 createSuccess = .true.
return
200 createSuccess =.false.
did a lot of polishing: - removed unnecessary "return" before end of subroutine or function: - changed undetermined array length (:) to (1:3) To prevent problems with some code analysing tools: - "3D oneliner loops" (with ";) only for "do" and "enddo" at the same time - removed line continuation in OMP statements made the makefile more flexible, removed heap-arrays switch
2011-09-13 21:24:06 +05:30
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
end function
!***********************************************************
! check if the input file for Abaqus contains part info
!***********************************************************
function IO_abaqus_hasNoPart(unit)
use prec, only: pInt
implicit none
integer(pInt) unit
integer(pInt), parameter :: maxNchunks = 1
integer(pInt), dimension(1+2*maxNchunks) :: pos
logical IO_abaqus_hasNoPart
character(len=300) line
IO_abaqus_hasNoPart = .true.
610 FORMAT(A300)
rewind(unit)
do
read(unit,610,END=620) line
pos = IO_stringPos(line,maxNchunks)
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
if (IO_lc(IO_stringValue(line,pos,1_pInt)) == '*part' ) then
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
IO_abaqus_hasNoPart = .false.
exit
endif
enddo
polishing !incremental update, wait for commit of damask_spectral.f90 before checking out
2011-10-18 14:52:33 +05:30
620 endfunction
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
!********************************************************************
! open existing file to given unit
! path to file is relative to working directory
!********************************************************************
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
logical function IO_open_file_stat(unit,relPath)
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
use prec, only: pInt
renamed interface from "mpie_interface" to "DAMASK_interface"
2011-05-11 22:31:03 +05:30
use DAMASK_interface
removed implicit type castings
2012-02-10 16:54:53 +05:30
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
implicit none
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
integer(pInt), intent(in) :: unit
character(len=*), intent(in) :: relPath
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
character(len=1024) path
integer(pInt) stat
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
path = trim(getSolverWorkingDirectoryName())//relPath
open(unit,status='old',iostat=stat,file=path)
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
IO_open_file_stat = (stat == 0_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
endfunction
!********************************************************************
! open existing file to given unit
! path to file is relative to working directory
!********************************************************************
subroutine IO_open_file(unit,relPath)
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
use prec, only: pInt
use DAMASK_interface
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
implicit none
integer(pInt), intent(in) :: unit
character(len=*), intent(in) :: relPath
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
character(len=1024) path
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
integer(pInt) stat
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
path = trim(getSolverWorkingDirectoryName())//relPath
open(unit,status='old',iostat=stat,file=path)
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
if (stat /= 0_pInt) call IO_error(100_pInt,ext_msg=path)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
endsubroutine
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
!********************************************************************
! open FEM inputfile to given unit
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
! AP: 12.07.10
! : changed the function to open *.inp_assembly, which is basically
! the input file without comment lines and possibly assembled includes
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
!********************************************************************
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
subroutine IO_open_inputFile(unit,model)
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
use prec, only: pReal, pInt
renamed interface from "mpie_interface" to "DAMASK_interface"
2011-05-11 22:31:03 +05:30
use DAMASK_interface
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
implicit none
quite some changes: # non-greedy memory allocation # generation of outputConstitutive to allow for script-based T16 extraction # exchange of phenomenological by more general phenopowerlaw # lattice is based on slip and twin families which can be treated as individual entities (switched on/off, separate hardening, etc.) # nicer debugging output # changed some error/warning codes # plus potentially some minor additional brushes here and there
2009-07-22 21:37:19 +05:30
integer(pInt), intent(in) :: unit
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
character(len=*), intent(in) :: model
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
character(len=1024) path
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
integer(pInt) stat
changed fileopen routines in IO.f90 to use solver dependant functions defined in the interface routines DID NOT test ABAQUS versions yet, not shure whether (V)GETOUTDIR returns path with or without terminating slash
2010-05-10 20:32:59 +05:30
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
if (FEsolver == 'Abaqus') then
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
path = trim(getSolverWorkingDirectoryName())//trim(model)//InputFileExtension
open(unit+1,status='old',iostat=stat,file=path)
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
if (stat /= 0_pInt) call IO_error(100_pInt,ext_msg=path)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
path = trim(getSolverWorkingDirectoryName())//trim(model)//InputFileExtension//'_assembly'
open(unit,iostat=stat,file=path)
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
if (stat /= 0_pInt) call IO_error(100_pInt,ext_msg=path)
if (IO_abaqus_assembleInputFile(unit,unit+1_pInt)) call IO_error(103_pInt) ! strip comments and concatenate any "include"s
removed implicit type castings
2012-02-10 16:54:53 +05:30
close(unit+1_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
else
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
path = trim(getSolverWorkingDirectoryName())//trim(model)//InputFileExtension
open(unit,status='old',iostat=stat,file=path)
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
if (stat /= 0_pInt) call IO_error(100_pInt,ext_msg=path)
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
endif
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
endsubroutine
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
!********************************************************************
! open FEM logfile to given unit
!********************************************************************
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
subroutine IO_open_logFile(unit)
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
use prec, only: pReal, pInt
renamed interface from "mpie_interface" to "DAMASK_interface"
2011-05-11 22:31:03 +05:30
use DAMASK_interface
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
implicit none
integer(pInt), intent(in) :: unit
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
character(len=1024) path
integer(pInt) stat
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
path = trim(getSolverWorkingDirectoryName())//trim(getSolverJobName())//LogFileExtension
open(unit,status='old',iostat=stat,file=path)
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
if (stat /= 0) call IO_error(100_pInt,ext_msg=path)
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
endsubroutine
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
quite some changes: # non-greedy memory allocation # generation of outputConstitutive to allow for script-based T16 extraction # exchange of phenomenological by more general phenopowerlaw # lattice is based on slip and twin families which can be treated as individual entities (switched on/off, separate hardening, etc.) # nicer debugging output # changed some error/warning codes # plus potentially some minor additional brushes here and there
2009-07-22 21:37:19 +05:30
!********************************************************************
! open (write) file related to current job
! but with different extension to given unit
!********************************************************************
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
logical function IO_open_jobFile_stat(unit,newExt)
quite some changes: # non-greedy memory allocation # generation of outputConstitutive to allow for script-based T16 extraction # exchange of phenomenological by more general phenopowerlaw # lattice is based on slip and twin families which can be treated as individual entities (switched on/off, separate hardening, etc.) # nicer debugging output # changed some error/warning codes # plus potentially some minor additional brushes here and there
2009-07-22 21:37:19 +05:30
use prec, only: pReal, pInt
renamed interface from "mpie_interface" to "DAMASK_interface"
2011-05-11 22:31:03 +05:30
use DAMASK_interface
quite some changes: # non-greedy memory allocation # generation of outputConstitutive to allow for script-based T16 extraction # exchange of phenomenological by more general phenopowerlaw # lattice is based on slip and twin families which can be treated as individual entities (switched on/off, separate hardening, etc.) # nicer debugging output # changed some error/warning codes # plus potentially some minor additional brushes here and there
2009-07-22 21:37:19 +05:30
implicit none
integer(pInt), intent(in) :: unit
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
character(len=*), intent(in) :: newExt
character(len=1024) path
integer(pInt) stat
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
path = trim(getSolverWorkingDirectoryName())//trim(getSolverJobName())//'.'//newExt
open(unit,status='old',iostat=stat,file=path)
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
IO_open_jobFile_stat = (stat == 0_pInt)
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
endfunction
quite some changes: # non-greedy memory allocation # generation of outputConstitutive to allow for script-based T16 extraction # exchange of phenomenological by more general phenopowerlaw # lattice is based on slip and twin families which can be treated as individual entities (switched on/off, separate hardening, etc.) # nicer debugging output # changed some error/warning codes # plus potentially some minor additional brushes here and there
2009-07-22 21:37:19 +05:30
introduced alternative location for material configuration. If present, the >>>JOB.materialConfig<<< file takes precedence over the standard >>>material.config<<<
2011-08-02 15:44:16 +05:30
!********************************************************************
! open (write) file related to current job
! but with different extension to given unit
!********************************************************************
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
subroutine IO_open_jobFile(unit,newExt)
introduced alternative location for material configuration. If present, the >>>JOB.materialConfig<<< file takes precedence over the standard >>>material.config<<<
2011-08-02 15:44:16 +05:30
use prec, only: pReal, pInt
use DAMASK_interface
implicit none
integer(pInt), intent(in) :: unit
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
character(len=*), intent(in) :: newExt
character(len=1024) path
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
integer(pInt) stat
introduced alternative location for material configuration. If present, the >>>JOB.materialConfig<<< file takes precedence over the standard >>>material.config<<<
2011-08-02 15:44:16 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
path = trim(getSolverWorkingDirectoryName())//trim(getSolverJobName())//'.'//newExt
open(unit,status='old',iostat=stat,file=path)
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
if (stat /= 0_pInt) call IO_error(100_pInt,ext_msg=path)
introduced alternative location for material configuration. If present, the >>>JOB.materialConfig<<< file takes precedence over the standard >>>material.config<<<
2011-08-02 15:44:16 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
endsubroutine
!********************************************************************
! open (write) file related to current job
! but with different extension to given unit
!********************************************************************
subroutine IO_write_jobFile(unit,newExt)
use prec, only: pReal, pInt
use DAMASK_interface
implicit none
integer(pInt), intent(in) :: unit
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
character(len=*), intent(in) :: newExt
character(len=1024) path
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
integer(pInt) stat
path = trim(getSolverWorkingDirectoryName())//trim(getSolverJobName())//'.'//newExt
open(unit,status='replace',iostat=stat,file=path)
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
if (stat /= 0_pInt) call IO_error(100_pInt,ext_msg=path)
introduced alternative location for material configuration. If present, the >>>JOB.materialConfig<<< file takes precedence over the standard >>>material.config<<<
2011-08-02 15:44:16 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
endsubroutine
introduced alternative location for material configuration. If present, the >>>JOB.materialConfig<<< file takes precedence over the standard >>>material.config<<<
2011-08-02 15:44:16 +05:30
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
!********************************************************************
! open (write) binary file related to current job
! but with different extension to given unit
!********************************************************************
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
subroutine IO_write_jobBinaryFile(unit,newExt,recMultiplier)
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
use prec, only: pReal, pInt
renamed interface from "mpie_interface" to "DAMASK_interface"
2011-05-11 22:31:03 +05:30
use DAMASK_interface
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
implicit none
integer(pInt), intent(in) :: unit
integer(pInt), intent(in), optional :: recMultiplier
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
character(len=*), intent(in) :: newExt
character(len=1024) path
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
integer(pInt) stat
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
path = trim(getSolverWorkingDirectoryName())//trim(getSolverJobName())//'.'//newExt
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
if (present(recMultiplier)) then
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
open(unit,status='replace',form='unformatted',access='direct',recl=pReal*recMultiplier,iostat=stat,file=path)
else
open(unit,status='replace',form='unformatted',access='direct',recl=pReal,iostat=stat,file=path)
endif
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
if (stat /= 0_pInt) call IO_error(100_pInt,ext_msg=path)
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
endsubroutine
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
!********************************************************************
! open (read) binary file related to restored job
! and with different extension to given unit
!********************************************************************
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
subroutine IO_read_jobBinaryFile(unit,newExt,jobName,recMultiplier)
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
use prec, only: pReal, pInt
renamed interface from "mpie_interface" to "DAMASK_interface"
2011-05-11 22:31:03 +05:30
use DAMASK_interface
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
implicit none
integer(pInt), intent(in) :: unit
integer(pInt), intent(in), optional :: recMultiplier
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
character(len=*), intent(in) :: newExt, jobName
character(len=1024) path
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
integer(pInt) stat
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
path = trim(getSolverWorkingDirectoryName())//trim(jobName)//'.'//newExt
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
if (present(recMultiplier)) then
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
open(unit,status='old',form='unformatted',access='direct',recl=pReal*recMultiplier,iostat=stat,file=path)
else
open(unit,status='old',form='unformatted',access='direct',recl=pReal,iostat=stat,file=path)
endif
if (stat /= 0) then
call IO_error(100_pInt,ext_msg=path)
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
endif
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
endsubroutine
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
hybridIA sampling included
2007-03-21 18:02:15 +05:30
!********************************************************************
! hybrid IA repetition counter
!********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
function hybridIA_reps(dV_V,steps,C)
hybridIA sampling included
2007-03-21 18:02:15 +05:30
use prec, only: pReal, pInt
implicit none
integer(pInt), intent(in), dimension(3) :: steps
integer(pInt) hybridIA_reps, phi1,Phi,phi2
real(pReal), intent(in), dimension(steps(3),steps(2),steps(1)) :: dV_V
real(pReal), intent(in) :: C
hybridIA_reps = 0_pInt
do phi1=1,steps(1)
do Phi =1,steps(2)
do phi2=1,steps(3)
hybridIA_reps = hybridIA_reps+nint(C*dV_V(phi2,Phi,phi1), pInt)
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
enddo
enddo
enddo
hybridIA sampling included
2007-03-21 18:02:15 +05:30
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endfunction
hybridIA sampling included
2007-03-21 18:02:15 +05:30
!********************************************************************
! hybrid IA sampling of ODFfile
!********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
function IO_hybridIA(Nast,ODFfileName)
hybridIA sampling included
2007-03-21 18:02:15 +05:30
use prec, only: pReal, pInt
implicit none
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
real(pReal), parameter :: pi = 3.14159265358979323846264338327950288419716939937510_pReal
real(pReal), parameter :: inRad = pi/180.0_pReal
hybridIA sampling included
2007-03-21 18:02:15 +05:30
character(len=*) ODFfileName
character(len=80) line
character(len=*), parameter :: fileFormat = '(A80)'
fixed allocatable hybridIA return value, finalized IO_lc
2007-03-26 20:33:21 +05:30
integer(pInt) i,j,bin,Nast,NnonZero,Nset,Nreps,reps,phi1,Phi,phi2
integer(pInt), dimension(7) :: pos
hybridIA sampling included
2007-03-21 18:02:15 +05:30
integer(pInt), dimension(3) :: steps
integer(pInt), dimension(:), allocatable :: binSet
real(pReal) center,sum_dV_V,prob,dg_0,C,lowerC,upperC,rnd
real(pReal), dimension(3) :: limits,deltas
real(pReal), dimension(:,:,:), allocatable :: dV_V
fixed allocatable hybridIA return value, finalized IO_lc
2007-03-26 20:33:21 +05:30
real(pReal), dimension(3,Nast) :: IO_hybridIA
hybridIA sampling included
2007-03-21 18:02:15 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
call IO_open_file(999_pInt,ODFfileName)
hybridIA sampling included
2007-03-21 18:02:15 +05:30
!--- parse header of ODF file ---
!--- limits in phi1, Phi, phi2 ---
read(999,fmt=fileFormat,end=100) line
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
pos = IO_stringPos(line,3_pInt)
hybridIA sampling included
2007-03-21 18:02:15 +05:30
if (pos(1).ne.3) goto 100
do i=1,3
removed implicit type castings
2012-02-10 16:54:53 +05:30
limits(i) = IO_floatValue(line,pos,i)*inRad
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
enddo
hybridIA sampling included
2007-03-21 18:02:15 +05:30
!--- deltas in phi1, Phi, phi2 ---
read(999,fmt=fileFormat,end=100) line
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
pos = IO_stringPos(line,3_pInt)
hybridIA sampling included
2007-03-21 18:02:15 +05:30
if (pos(1).ne.3) goto 100
do i=1,3
removed implicit type castings
2012-02-10 16:54:53 +05:30
deltas(i) = IO_floatValue(line,pos,i)*inRad
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
enddo
hybridIA sampling included
2007-03-21 18:02:15 +05:30
steps = nint(limits/deltas,pInt)
allocate(dV_V(steps(3),steps(2),steps(1)))
!--- box boundary/center at origin? ---
read(999,fmt=fileFormat,end=100) line
if (index(IO_lc(line),'bound')>0) then
center = 0.5_pReal
else
center = 0.0_pReal
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endif
hybridIA sampling included
2007-03-21 18:02:15 +05:30
!--- skip blank line ---
read(999,fmt=fileFormat,end=100) line
sum_dV_V = 0.0_pReal
dV_V = 0.0_pReal
dg_0 = deltas(1)*deltas(3)*2.0_pReal*sin(deltas(2)/2.0_pReal)
NnonZero = 0_pInt
do phi1=1,steps(1)
do Phi=1,steps(2)
do phi2=1,steps(3)
changed file read formats from explicit types to more general (and not harmful) '*'
2008-03-12 19:23:00 +05:30
read(999,fmt=*,end=100) prob
hybridIA sampling included
2007-03-21 18:02:15 +05:30
if (prob > 0.0_pReal) then
NnonZero = NnonZero+1
sum_dV_V = sum_dV_V+prob
else
prob = 0.0_pReal
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endif
hybridIA sampling included
2007-03-21 18:02:15 +05:30
dV_V(phi2,Phi,phi1) = prob*dg_0*sin((Phi-1.0_pReal+center)*deltas(2))
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
enddo
enddo
enddo
replaces Nast with Nreps in hybridIA formerly, if the number of orientations requested was (far) less than the number of nonzero bins in the texture representation, you only got a subset of low phi1 values..! (that is the top corner in Euler space -- Hi Bing ;-) )
2009-08-13 18:51:22 +05:30
hybridIA sampling included
2007-03-21 18:02:15 +05:30
dV_V = dV_V/sum_dV_V ! normalize to 1
!--- now fix bounds ---
put some comments to hybridIA
2009-10-14 18:51:03 +05:30
Nset = max(Nast,NnonZero) ! if less than non-zero voxel count requested, sample at least that much
hybridIA sampling included
2007-03-21 18:02:15 +05:30
lowerC = 0.0_pReal
upperC = real(Nset, pReal)
do while (hybridIA_reps(dV_V,steps,upperC) < Nset)
lowerC = upperC
upperC = upperC*2.0_pReal
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
enddo
hybridIA sampling included
2007-03-21 18:02:15 +05:30
!--- binary search for best C ---
do
C = (upperC+lowerC)/2.0_pReal
Nreps = hybridIA_reps(dV_V,steps,C)
if (abs(upperC-lowerC) < upperC*1.0e-14_pReal) then
C = upperC
Nreps = hybridIA_reps(dV_V,steps,C)
exit
elseif (Nreps < Nset) then
lowerC = C
elseif (Nreps > Nset) then
upperC = C
else
exit
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endif
enddo
replaces Nast with Nreps in hybridIA formerly, if the number of orientations requested was (far) less than the number of nonzero bins in the texture representation, you only got a subset of low phi1 values..! (that is the top corner in Euler space -- Hi Bing ;-) )
2009-08-13 18:51:22 +05:30
hybridIA sampling included
2007-03-21 18:02:15 +05:30
allocate(binSet(Nreps))
removed implicit type castings
2012-02-10 16:54:53 +05:30
bin = 0_pInt ! bin counter
hybridIA sampling included
2007-03-21 18:02:15 +05:30
i = 1 ! set counter
do phi1=1,steps(1)
do Phi=1,steps(2)
do phi2=1,steps(3)
reps = nint(C*dV_V(phi2,Phi,phi1), pInt)
binSet(i:i+reps-1) = bin
bin = bin+1 ! advance bin
i = i+reps ! advance set
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
enddo
enddo
enddo
hybridIA sampling included
2007-03-21 18:02:15 +05:30
do i=1,Nast
if (i < Nast) then
call random_number(rnd)
replaces Nast with Nreps in hybridIA formerly, if the number of orientations requested was (far) less than the number of nonzero bins in the texture representation, you only got a subset of low phi1 values..! (that is the top corner in Euler space -- Hi Bing ;-) )
2009-08-13 18:51:22 +05:30
j = nint(rnd*(Nreps-i)+i+0.5_pReal,pInt)
hybridIA sampling included
2007-03-21 18:02:15 +05:30
else
j = i
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endif
hybridIA sampling included
2007-03-21 18:02:15 +05:30
bin = binSet(j)
removed implicit type castings
2012-02-10 16:54:53 +05:30
IO_hybridIA(1,i) = deltas(1)*(real(mod(bin/(steps(3)*steps(2)),steps(1)),pReal)+center) ! phi1
IO_hybridIA(2,i) = deltas(2)*(real(mod(bin/ steps(3) ,steps(2)),pReal)+center) ! Phi
IO_hybridIA(3,i) = deltas(3)*(real(mod(bin ,steps(3)),pReal)+center) ! phi2
hybridIA sampling included
2007-03-21 18:02:15 +05:30
binSet(j) = binSet(i)
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
enddo
hybridIA sampling included
2007-03-21 18:02:15 +05:30
close(999)
return
! on error
removed implicit type castings
2012-02-10 16:54:53 +05:30
100 IO_hybridIA = -1.0_pReal
hybridIA sampling included
2007-03-21 18:02:15 +05:30
close(999)
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endfunction
hybridIA sampling included
2007-03-21 18:02:15 +05:30
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
!********************************************************************
! identifies lines without content
!********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
pure function IO_isBlank (line)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
use prec, only: pInt
implicit none
character(len=*), intent(in) :: line
character(len=*), parameter :: blank = achar(32)//achar(9)//achar(10)//achar(13) ! whitespaces
character(len=*), parameter :: comment = achar(35) ! comment id '#'
integer(pInt) posNonBlank, posComment
logical IO_isBlank
posNonBlank = verify(line,blank)
posComment = scan(line,comment)
IO_isBlank = posNonBlank == 0 .or. posNonBlank == posComment
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endfunction
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
!********************************************************************
! get tagged content of line
!********************************************************************
did a lot of polishing: - removed unnecessary "return" before end of subroutine or function: - changed undetermined array length (:) to (1:3) To prevent problems with some code analysing tools: - "3D oneliner loops" (with ";) only for "do" and "enddo" at the same time - removed line continuation in OMP statements made the makefile more flexible, removed heap-arrays switch
2011-09-13 21:24:06 +05:30
pure function IO_getTag (line,openChar,closeChar)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
use prec, only: pInt
implicit none
character(len=*), intent(in) :: line,openChar,closeChar
character(len=*), parameter :: sep=achar(32)//achar(9)//achar(10)//achar(13) ! whitespaces
character(len=len_trim(line)) IO_getTag
integer(pInt) left,right
IO_getTag = ''
left = scan(line,openChar)
right = scan(line,closeChar)
if (left == verify(line,sep) .and. right > left) & ! openChar is first and closeChar occurs
IO_getTag = line(left+1:right-1)
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endfunction
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
!*********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
function IO_countSections(file,part)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
!*********************************************************************
use prec, only: pInt
implicit none
!* Definition of variables
integer(pInt), intent(in) :: file
character(len=*), intent(in) :: part
integer(pInt) IO_countSections
character(len=1024) line
IO_countSections = 0
line = ''
rewind(file)
do while (IO_getTag(line,'<','>') /= part) ! search for part
read(file,'(a1024)',END=100) line
enddo
do
read(file,'(a1024)',END=100) line
if (IO_isBlank(line)) cycle ! skip empty lines
if (IO_getTag(line,'<','>') /= '') exit ! stop at next part
if (IO_getTag(line,'[',']') /= '') & ! found [section] identifier
IO_countSections = IO_countSections + 1
enddo
polishing !incremental update, wait for commit of damask_spectral.f90 before checking out
2011-10-18 14:52:33 +05:30
100 endfunction
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
!*********************************************************************
! return array of myTag counts within <part> for at most N[sections]
!*********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
function IO_countTagInPart(file,part,myTag,Nsections)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
use prec, only: pInt
implicit none
!* Definition of variables
integer(pInt), intent(in) :: file, Nsections
character(len=*), intent(in) :: part, myTag
integer(pInt), dimension(Nsections) :: IO_countTagInPart, counter
integer(pInt), parameter :: maxNchunks = 1
integer(pInt), dimension(1+2*maxNchunks) :: positions
integer(pInt) section
character(len=1024) line,tag
counter = 0_pInt
section = 0_pInt
line = ''
rewind(file)
do while (IO_getTag(line,'<','>') /= part) ! search for part
read(file,'(a1024)',END=100) line
enddo
do
read(file,'(a1024)',END=100) line
if (IO_isBlank(line)) cycle ! skip empty lines
if (IO_getTag(line,'<','>') /= '') exit ! stop at next part
if (IO_getTag(line,'[',']') /= '') & ! found [section] identifier
section = section + 1
if (section > 0) then
positions = IO_stringPos(line,maxNchunks)
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
tag = IO_lc(IO_stringValue(line,positions,1_pInt)) ! extract key
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
if (tag == myTag) & ! match
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
counter(section) = counter(section) + 1_pInt
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
endif
enddo
100 IO_countTagInPart = counter
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endfunction
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
new function IO_spotTagInPart
2009-04-03 16:00:18 +05:30
!*********************************************************************
! return array of myTag presence within <part> for at most N[sections]
!*********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
function IO_spotTagInPart(file,part,myTag,Nsections)
new function IO_spotTagInPart
2009-04-03 16:00:18 +05:30
use prec, only: pInt
implicit none
!* Definition of variables
integer(pInt), intent(in) :: file, Nsections
character(len=*), intent(in) :: part, myTag
logical, dimension(Nsections) :: IO_spotTagInPart
integer(pInt), parameter :: maxNchunks = 1
integer(pInt), dimension(1+2*maxNchunks) :: positions
integer(pInt) section
character(len=1024) line,tag
IO_spotTagInPart = .false. ! assume to nowhere spot tag
section = 0_pInt
line = ''
rewind(file)
do while (IO_getTag(line,'<','>') /= part) ! search for part
read(file,'(a1024)',END=100) line
enddo
do
read(file,'(a1024)',END=100) line
if (IO_isBlank(line)) cycle ! skip empty lines
if (IO_getTag(line,'<','>') /= '') exit ! stop at next part
if (IO_getTag(line,'[',']') /= '') & ! found [section] identifier
section = section + 1
if (section > 0) then
positions = IO_stringPos(line,maxNchunks)
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
tag = IO_lc(IO_stringValue(line,positions,1_pInt)) ! extract key
new function IO_spotTagInPart
2009-04-03 16:00:18 +05:30
if (tag == myTag) & ! match
IO_spotTagInPart(section) = .true.
endif
enddo
polishing !incremental update, wait for commit of damask_spectral.f90 before checking out
2011-10-18 14:52:33 +05:30
100 endfunction
new function IO_spotTagInPart
2009-04-03 16:00:18 +05:30
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
!********************************************************************
! locate at most N space-separated parts in line
corrected line parsing now counts all parts
2009-12-15 21:33:53 +05:30
! return array containing number of parts in line and
! the left/right positions of at most N to be used by IO_xxxVal
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
!********************************************************************
corrected line parsing now counts all parts
2009-12-15 21:33:53 +05:30
! pure function IO_stringPos (line,N)
function IO_stringPos (line,N)
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
use prec, only: pReal,pInt
implicit none
rendered most of the functions to be "pure"
2009-01-16 20:55:37 +05:30
character(len=*), intent(in) :: line
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
character(len=*), parameter :: sep=achar(44)//achar(32)//achar(9)//achar(10)//achar(13) ! comma and whitespaces
changed Line Ending to Unix style (LF)
2009-01-20 00:40:58 +05:30
integer(pInt), intent(in) :: N
corrected line parsing now counts all parts
2009-12-15 21:33:53 +05:30
integer(pInt) left,right
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
integer(pInt) IO_stringPos(1+N*2)
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 20:08:22 +05:30
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
IO_stringPos = -1
IO_stringPos(1) = 0
corrected line parsing now counts all parts
2009-12-15 21:33:53 +05:30
right = 0
do while (verify(line(right+1:),sep)>0)
left = right + verify(line(right+1:),sep)
right = left + scan(line(left:),sep) - 2
* comments are now recognized by IO even when they do not appear at the beginning of the line * no need for the zero entries in the material.config anymore
2011-05-30 14:39:19 +05:30
if ( line(left:left) == '#' ) then
exit
endif
corrected line parsing now counts all parts
2009-12-15 21:33:53 +05:30
if ( IO_stringPos(1)<N ) then
IO_stringPos(1+IO_stringPos(1)*2+1) = left
IO_stringPos(1+IO_stringPos(1)*2+2) = right
endif
IO_stringPos(1) = IO_stringPos(1)+1
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
enddo
corrected line parsing now counts all parts
2009-12-15 21:33:53 +05:30
polishing !incremental update, wait for commit of damask_spectral.f90 before checking out
2011-10-18 14:52:33 +05:30
endfunction
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
!********************************************************************
! read string value at pos from line
!********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
pure function IO_stringValue (line,positions,pos)
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
use prec, only: pReal,pInt
implicit none
rendered most of the functions to be "pure"
2009-01-16 20:55:37 +05:30
character(len=*), intent(in) :: line
integer(pInt), intent(in) :: positions(*),pos
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
character(len=1+positions(pos*2+1)-positions(pos*2)) IO_stringValue
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 20:08:22 +05:30
if (positions(1) < pos) then
IO_stringValue = ''
fixed flaws in fixedFormat number extraction
2007-03-28 15:30:49 +05:30
else
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 20:08:22 +05:30
IO_stringValue = line(positions(pos*2):positions(pos*2+1))
fixed flaws in fixedFormat number extraction
2007-03-28 15:30:49 +05:30
endif
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endfunction
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
added fixed format line reading
2007-03-21 20:15:03 +05:30
!********************************************************************
polished some syntax problems
2007-03-21 20:19:21 +05:30
! read string value at pos from fixed format line
added fixed format line reading
2007-03-21 20:15:03 +05:30
!********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
pure function IO_fixedStringValue (line,ends,pos)
added fixed format line reading
2007-03-21 20:15:03 +05:30
use prec, only: pReal,pInt
implicit none
rendered most of the functions to be "pure"
2009-01-16 20:55:37 +05:30
character(len=*), intent(in) :: line
integer(pInt), intent(in) :: ends(*),pos
fixed flaws in fixedFormat number extraction
2007-03-28 15:30:49 +05:30
character(len=ends(pos+1)-ends(pos)) IO_fixedStringValue
added fixed format line reading
2007-03-21 20:15:03 +05:30
fixed flaws in fixedFormat number extraction
2007-03-28 15:30:49 +05:30
IO_fixedStringValue = line(ends(pos)+1:ends(pos+1))
added fixed format line reading
2007-03-21 20:15:03 +05:30
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endfunction
added fixed format line reading
2007-03-21 20:15:03 +05:30
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
!********************************************************************
! read float value at pos from line
!********************************************************************
polishing !incremental update, wait for commit of damask_spectral.f90 before checking out
2011-10-18 14:52:33 +05:30
pure function IO_floatValue (line,positions,myPos)
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
use prec, only: pReal,pInt
implicit none
rendered most of the functions to be "pure"
2009-01-16 20:55:37 +05:30
character(len=*), intent(in) :: line
polishing !incremental update, wait for commit of damask_spectral.f90 before checking out
2011-10-18 14:52:33 +05:30
integer(pInt), intent(in) :: positions(*),myPos
changed Line Ending to Unix style (LF)
2009-01-20 00:40:58 +05:30
real(pReal) IO_floatValue
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
polishing !incremental update, wait for commit of damask_spectral.f90 before checking out
2011-10-18 14:52:33 +05:30
if (positions(1) < myPos) then
now safe to request values from config-parsing which are not present. 3rd and up value of key "xyz" would return "0" or "0.0" in this example: xyz 1.234 4.321
2009-06-19 12:39:39 +05:30
IO_floatValue = 0.0_pReal
else
polishing !incremental update, wait for commit of damask_spectral.f90 before checking out
2011-10-18 14:52:33 +05:30
read(UNIT=line(positions(myPos*2):positions(myPos*2+1)),ERR=100,FMT=*) IO_floatValue
fixed flaws in fixedFormat number extraction
2007-03-28 15:30:49 +05:30
endif
now safe to request values from config-parsing which are not present. 3rd and up value of key "xyz" would return "0" or "0.0" in this example: xyz 1.234 4.321
2009-06-19 12:39:39 +05:30
return
changed _Xvalue() return codes on error to huge()
2007-04-26 18:10:06 +05:30
100 IO_floatValue = huge(1.0_pReal)
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endfunction
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
added fixed format line reading
2007-03-21 20:15:03 +05:30
!********************************************************************
polished some syntax problems
2007-03-21 20:19:21 +05:30
! read float value at pos from fixed format line
added fixed format line reading
2007-03-21 20:15:03 +05:30
!********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
pure function IO_fixedFloatValue (line,ends,pos)
added fixed format line reading
2007-03-21 20:15:03 +05:30
use prec, only: pReal,pInt
implicit none
rendered most of the functions to be "pure"
2009-01-16 20:55:37 +05:30
character(len=*), intent(in) :: line
integer(pInt), intent(in) :: ends(*),pos
changed Line Ending to Unix style (LF)
2009-01-20 00:40:58 +05:30
real(pReal) IO_fixedFloatValue
added fixed format line reading
2007-03-21 20:15:03 +05:30
changed _Xvalue() return codes on error to huge()
2007-04-26 18:10:06 +05:30
read(UNIT=line(ends(pos-1)+1:ends(pos)),ERR=100,FMT=*) IO_fixedFloatValue
added fixed format line reading
2007-03-21 20:15:03 +05:30
return
changed _Xvalue() return codes on error to huge()
2007-04-26 18:10:06 +05:30
100 IO_fixedFloatValue = huge(1.0_pReal)
added fixed format line reading
2007-03-21 20:15:03 +05:30
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endfunction
added fixed format line reading
2007-03-21 20:15:03 +05:30
!********************************************************************
polished some syntax problems
2007-03-21 20:19:21 +05:30
! read float x.y+z value at pos from format line line
added fixed format line reading
2007-03-21 20:15:03 +05:30
!********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
pure function IO_fixedNoEFloatValue (line,ends,pos)
added fixed format line reading
2007-03-21 20:15:03 +05:30
use prec, only: pReal,pInt
implicit none
rendered most of the functions to be "pure"
2009-01-16 20:55:37 +05:30
character(len=*), intent(in) :: line
changed Line Ending to Unix style (LF)
2009-01-20 00:40:58 +05:30
integer(pInt), intent(in) :: ends(*),pos
integer(pInt) pos_exp,expon
real(pReal) IO_fixedNoEFloatValue,base
added fixed format line reading
2007-03-21 20:15:03 +05:30
fixed flaws in fixedFormat number extraction
2007-03-28 15:30:49 +05:30
pos_exp = scan(line(ends(pos)+1:ends(pos+1)),'+-',back=.true.)
added fixed format line reading
2007-03-21 20:15:03 +05:30
if (pos_exp > 1) then
changed file read formats from explicit types to more general (and not harmful) '*'
2008-03-12 19:23:00 +05:30
read(UNIT=line(ends(pos)+1:ends(pos)+pos_exp-1),ERR=100,FMT=*) base
read(UNIT=line(ends(pos)+pos_exp:ends(pos+1)),ERR=100,FMT=*) expon
added fixed format line reading
2007-03-21 20:15:03 +05:30
else
changed _Xvalue() return codes on error to huge()
2007-04-26 18:10:06 +05:30
read(UNIT=line(ends(pos)+1:ends(pos+1)),ERR=100,FMT=*) base
change of expon in fixedFloatNoE to type int, general polishing
2007-03-21 22:34:10 +05:30
expon = 0_pInt
added fixed format line reading
2007-03-21 20:15:03 +05:30
endif
IO_fixedNoEFloatValue = base*10.0_pReal**expon
return
changed _Xvalue() return codes on error to huge()
2007-04-26 18:10:06 +05:30
100 IO_fixedNoEFloatValue = huge(1.0_pReal)
added fixed format line reading
2007-03-21 20:15:03 +05:30
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endfunction
added fixed format line reading
2007-03-21 20:15:03 +05:30
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
!********************************************************************
! read int value at pos from line
!********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
pure function IO_intValue (line,positions,pos)
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
use prec, only: pReal,pInt
implicit none
rendered most of the functions to be "pure"
2009-01-16 20:55:37 +05:30
character(len=*), intent(in) :: line
integer(pInt), intent(in) :: positions(*),pos
changed Line Ending to Unix style (LF)
2009-01-20 00:40:58 +05:30
integer(pInt) IO_intValue
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
now safe to request values from config-parsing which are not present. 3rd and up value of key "xyz" would return "0" or "0.0" in this example: xyz 1.234 4.321
2009-06-19 12:39:39 +05:30
if (positions(1) < pos) then
IO_intValue = 0_pInt
else
changed file read formats from explicit types to more general (and not harmful) '*'
2008-03-12 19:23:00 +05:30
read(UNIT=line(positions(pos*2):positions(pos*2+1)),ERR=100,FMT=*) IO_intValue
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 20:08:22 +05:30
endif
now safe to request values from config-parsing which are not present. 3rd and up value of key "xyz" would return "0" or "0.0" in this example: xyz 1.234 4.321
2009-06-19 12:39:39 +05:30
return
changed _Xvalue() return codes on error to huge()
2007-04-26 18:10:06 +05:30
100 IO_intValue = huge(1_pInt)
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endfunction
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
added fixed format line reading
2007-03-21 20:15:03 +05:30
!********************************************************************
! read int value at pos from fixed format line
!********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
pure function IO_fixedIntValue (line,ends,pos)
added fixed format line reading
2007-03-21 20:15:03 +05:30
use prec, only: pReal,pInt
implicit none
rendered most of the functions to be "pure"
2009-01-16 20:55:37 +05:30
character(len=*), intent(in) :: line
integer(pInt), intent(in) :: ends(*),pos
changed Line Ending to Unix style (LF)
2009-01-20 00:40:58 +05:30
integer(pInt) IO_fixedIntValue
added fixed format line reading
2007-03-21 20:15:03 +05:30
changed file read formats from explicit types to more general (and not harmful) '*'
2008-03-12 19:23:00 +05:30
read(UNIT=line(ends(pos)+1:ends(pos+1)),ERR=100,FMT=*) IO_fixedIntValue
added fixed format line reading
2007-03-21 20:15:03 +05:30
return
changed _Xvalue() return codes on error to huge()
2007-04-26 18:10:06 +05:30
100 IO_fixedIntValue = huge(1_pInt)
added fixed format line reading
2007-03-21 20:15:03 +05:30
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endfunction
added fixed format line reading
2007-03-21 20:15:03 +05:30
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 20:08:22 +05:30
!********************************************************************
! change character in line to lower case
!********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
pure function IO_lc (line)
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
use prec, only: pInt
implicit none
rendered most of the functions to be "pure"
2009-01-16 20:55:37 +05:30
character (len=*), intent(in) :: line
hybridIA sampling included
2007-03-21 18:02:15 +05:30
character (len=len(line)) IO_lc
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
integer(pInt) i
hybridIA sampling included
2007-03-21 18:02:15 +05:30
IO_lc = line
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 20:08:22 +05:30
do i=1,len(line)
if(64<iachar(line(i:i)) .and. iachar(line(i:i))<91) IO_lc(i:i)=achar(iachar(line(i:i))+32)
enddo
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endfunction
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
!********************************************************************
hybridIA sampling included
2007-03-21 18:02:15 +05:30
! in place change of character in line to lower case
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
!********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
subroutine IO_lcInplace (line)
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
use prec, only: pInt
implicit none
character (len=*) line
fixed compiler error in IO_lcInplace
2007-03-28 14:01:12 +05:30
character (len=len(line)) IO_lc
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
integer(pInt) i
fixed compiler error in IO_lcInplace
2007-03-28 14:01:12 +05:30
IO_lc = line
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 20:08:22 +05:30
do i=1,len(line)
if(64<iachar(line(i:i)) .and. iachar(line(i:i))<91) IO_lc(i:i)=achar(iachar(line(i:i))+32)
enddo
fixed compiler error in IO_lcInplace
2007-03-28 14:01:12 +05:30
line = IO_lc
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endsubroutine
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
added subroutine 'IO_skipChunks' that is needed by mesh.f90
2009-04-03 12:34:31 +05:30
!********************************************************************
! read on in file to skip (at least) N chunks (may be over multiple lines)
!********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
subroutine IO_skipChunks (unit,N)
added subroutine 'IO_skipChunks' that is needed by mesh.f90
2009-04-03 12:34:31 +05:30
use prec, only: pReal,pInt
implicit none
integer(pInt) remainingChunks,unit,N
integer(pInt), parameter :: maxNchunks = 64
integer(pInt), dimension(1+2*maxNchunks) :: pos
character(len=300) line
remainingChunks = N
do while (remainingChunks > 0)
read(unit,'(A300)',end=100) line
pos = IO_stringPos(line,maxNchunks)
remainingChunks = remainingChunks - pos(1)
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
enddo
polishing !incremental update, wait for commit of damask_spectral.f90 before checking out
2011-10-18 14:52:33 +05:30
100 endsubroutine
added subroutine 'IO_skipChunks' that is needed by mesh.f90
2009-04-03 12:34:31 +05:30
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
!********************************************************************
! extract value from key=value pair and check whether key matches
!********************************************************************
pure function IO_extractValue (line,key)
added subroutine 'IO_skipChunks' that is needed by mesh.f90
2009-04-03 12:34:31 +05:30
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
use prec, only: pReal,pInt
implicit none
character(len=*), intent(in) :: line,key
character(len=*), parameter :: sep = achar(61) ! '='
integer(pInt) pos
character(len=300) IO_extractValue
IO_extractValue = ''
pos = scan(line,sep)
if (pos > 0 .and. line(:pos-1) == key(:pos-1)) & ! key matches expected key
IO_extractValue = line(pos+1:) ! extract value
endfunction
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 20:08:22 +05:30
!********************************************************************
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
! count lines containig data up to next *keyword
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
! AP: changed the function to neglect comment lines between keyword definitions.
! : is not changed back to the original version since *.inp_assembly does not
! : contain any comment lines (12.07.2010)
added function IO_countContinousIntValues to read number of items in set definitions: 123 to 423 or something like 1 3 5 6 7 9 10 c 12 14
2007-10-15 19:25:52 +05:30
!********************************************************************
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
function IO_countDataLines (unit)
added function IO_countContinousIntValues to read number of items in set definitions: 123 to 423 or something like 1 3 5 6 7 9 10 c 12 14
2007-10-15 19:25:52 +05:30
use prec, only: pReal,pInt
implicit none
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
integer(pInt) IO_countDataLines,unit
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
integer(pInt), parameter :: maxNchunks = 1
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
integer(pInt), dimension(1+2*maxNchunks) :: pos
character(len=300) line,tmp
added function IO_countContinousIntValues to read number of items in set definitions: 123 to 423 or something like 1 3 5 6 7 9 10 c 12 14
2007-10-15 19:25:52 +05:30
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
IO_countDataLines = 0
added function IO_countContinousIntValues to read number of items in set definitions: 123 to 423 or something like 1 3 5 6 7 9 10 c 12 14
2007-10-15 19:25:52 +05:30
do
read(unit,'(A300)',end=100) line
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
pos = IO_stringPos(line,maxNchunks)
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
tmp = IO_lc(IO_stringValue(line,pos,1_pInt))
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
if (tmp(1:1) == '*' .and. tmp(2:2) /= '*') then ! found keyword
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
exit
added function IO_countContinousIntValues to read number of items in set definitions: 123 to 423 or something like 1 3 5 6 7 9 10 c 12 14
2007-10-15 19:25:52 +05:30
else
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
if (tmp(2:2) /= '*') IO_countDataLines = IO_countDataLines + 1_pInt
added function IO_countContinousIntValues to read number of items in set definitions: 123 to 423 or something like 1 3 5 6 7 9 10 c 12 14
2007-10-15 19:25:52 +05:30
endif
enddo
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
100 backspace(unit)
endfunction
!********************************************************************
! count items in consecutive lines
! Marc: ints concatenated by "c" as last char or range of values a "to" b
! Abaqus: triplet of start,stop,inc
!********************************************************************
function IO_countContinousIntValues (unit)
renamed interface from "mpie_interface" to "DAMASK_interface"
2011-05-11 22:31:03 +05:30
use DAMASK_interface
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
use prec, only: pReal,pInt
implicit none
allow for longer pathes when openening files corrected error messages concerning file openening
2010-02-18 13:59:57 +05:30
integer(pInt) unit,l,count
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
integer(pInt) IO_countContinousIntValues
parsing of spectral geom file now accepts two-dimensional microstructure definition (plus range indications by 'x to y') For instance: 1 2 3 4 5 6 7 8 9 10 11 12 ...
2012-01-12 22:31:24 +05:30
integer(pInt), parameter :: maxNchunks = 8192
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
integer(pInt), dimension(1+2*maxNchunks) :: pos
parsing of spectral geom file now accepts two-dimensional microstructure definition (plus range indications by 'x to y') For instance: 1 2 3 4 5 6 7 8 9 10 11 12 ...
2012-01-12 22:31:24 +05:30
character(len=65536) line
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
IO_countContinousIntValues = 0_pInt
select case (FEsolver)
parsing of spectral geom file now accepts two-dimensional microstructure definition (plus range indications by 'x to y') For instance: 1 2 3 4 5 6 7 8 9 10 11 12 ...
2012-01-12 22:31:24 +05:30
case ('Marc','Spectral')
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
do
read(unit,'(A300)',end=100) line
pos = IO_stringPos(line,maxNchunks)
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
if (IO_lc(IO_stringValue(line,pos,2_pInt)) == 'to' ) then ! found range indicator
IO_countContinousIntValues = 1_pInt + IO_intValue(line,pos,3_pInt) - IO_intValue(line,pos,1_pInt)
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
exit ! only one single range indicator allowed
else
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
IO_countContinousIntValues = IO_countContinousIntValues+pos(1)-1_pInt ! add line's count when assuming 'c'
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
if ( IO_lc(IO_stringValue(line,pos,pos(1))) /= 'c' ) then ! line finished, read last value
IO_countContinousIntValues = IO_countContinousIntValues+1
exit ! data ended
endif
endif
enddo
case('Abaqus')
count = IO_countDataLines(unit)
do l = 1,count
backspace(unit)
enddo
do l = 1,count
read(unit,'(A300)',end=100) line
pos = IO_stringPos(line,maxNchunks)
IO_countContinousIntValues = IO_countContinousIntValues + 1 + & ! assuming range generation
removed implicit type castings
2012-02-10 16:54:53 +05:30
(IO_intValue(line,pos,2_pInt)-IO_intValue(line,pos,1_pInt))/&
max(1_pInt,IO_intValue(line,pos,3_pInt))
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
enddo
endselect
polishing !incremental update, wait for commit of damask_spectral.f90 before checking out
2011-10-18 14:52:33 +05:30
100 endfunction
added function IO_countContinousIntValues to read number of items in set definitions: 123 to 423 or something like 1 3 5 6 7 9 10 c 12 14
2007-10-15 19:25:52 +05:30
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
!********************************************************************
! return integer list corrsponding to items in consecutive lines
! Marc: ints concatenated by "c" as last char, range of a "to" b, or named set
! Abaqus: triplet of start,stop,inc or named set
!********************************************************************
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
function IO_continousIntValues (unit,maxN,lookupName,lookupMap,lookupMaxN)
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 20:08:22 +05:30
renamed interface from "mpie_interface" to "DAMASK_interface"
2011-05-11 22:31:03 +05:30
use DAMASK_interface
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 20:08:22 +05:30
use prec, only: pReal,pInt
implicit none
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
integer(pInt) unit,maxN,i,j,l,count,first,last
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 20:08:22 +05:30
integer(pInt), dimension(1+maxN) :: IO_continousIntValues
parsing of spectral geom file now accepts two-dimensional microstructure definition (plus range indications by 'x to y') For instance: 1 2 3 4 5 6 7 8 9 10 11 12 ...
2012-01-12 22:31:24 +05:30
integer(pInt), parameter :: maxNchunks = 8192_pInt
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
integer(pInt), dimension(1+2*maxNchunks) :: pos
fixed broken IO_continousIntValues to be aware of element (or other) sets -This line, and those below, will be ignored-- M IO.f90
2007-10-23 18:38:27 +05:30
character(len=64), dimension(:) :: lookupName
integer(pInt) :: lookupMaxN
integer(pInt), dimension(:,:) :: lookupMap
parsing of spectral geom file now accepts two-dimensional microstructure definition (plus range indications by 'x to y') For instance: 1 2 3 4 5 6 7 8 9 10 11 12 ...
2012-01-12 22:31:24 +05:30
character(len=65536) line
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
logical rangeGeneration
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 20:08:22 +05:30
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
IO_continousIntValues = 0
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
rangeGeneration = .false.
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
select case (FEsolver)
parsing of spectral geom file now accepts two-dimensional microstructure definition (plus range indications by 'x to y') For instance: 1 2 3 4 5 6 7 8 9 10 11 12 ...
2012-01-12 22:31:24 +05:30
case ('Marc','Spectral')
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
do
parsing of spectral geom file now accepts two-dimensional microstructure definition (plus range indications by 'x to y') For instance: 1 2 3 4 5 6 7 8 9 10 11 12 ...
2012-01-12 22:31:24 +05:30
read(unit,'(A65536)',end=100) line
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
pos = IO_stringPos(line,maxNchunks)
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
if (verify(IO_stringValue(line,pos,1_pInt),'0123456789') > 0) then ! a non-int, i.e. set name
do i = 1_pInt,lookupMaxN ! loop over known set names
if (IO_stringValue(line,pos,1_pInt) == lookupName(i)) then ! found matching name
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
IO_continousIntValues = lookupMap(:,i) ! return resp. entity list
exit
endif
enddo
fixed broken IO_continousIntValues to be aware of element (or other) sets -This line, and those below, will be ignored-- M IO.f90
2007-10-23 18:38:27 +05:30
exit
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
else if (pos(1) > 2_pInt .and. IO_lc(IO_stringValue(line,pos,2_pInt)) == 'to' ) then ! found range indicator
do i = IO_intValue(line,pos,1_pInt),IO_intValue(line,pos,3_pInt)
IO_continousIntValues(1) = IO_continousIntValues(1) + 1_pInt
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
IO_continousIntValues(1+IO_continousIntValues(1)) = i
enddo
exit
else
do i = 1,pos(1)-1 ! interpret up to second to last value
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
IO_continousIntValues(1) = IO_continousIntValues(1) + 1_pInt
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
IO_continousIntValues(1+IO_continousIntValues(1)) = IO_intValue(line,pos,i)
enddo
if ( IO_lc(IO_stringValue(line,pos,pos(1))) /= 'c' ) then ! line finished, read last value
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
IO_continousIntValues(1) = IO_continousIntValues(1) + 1_pInt
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
IO_continousIntValues(1+IO_continousIntValues(1)) = IO_intValue(line,pos,pos(1))
exit
endif
fixed broken IO_continousIntValues to be aware of element (or other) sets -This line, and those below, will be ignored-- M IO.f90
2007-10-23 18:38:27 +05:30
endif
enddo
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
case('Abaqus')
count = IO_countDataLines(unit)
do l = 1,count
backspace(unit)
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 20:08:22 +05:30
enddo
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
! check if the element values in the elset are auto generated
backspace(unit)
parsing of spectral geom file now accepts two-dimensional microstructure definition (plus range indications by 'x to y') For instance: 1 2 3 4 5 6 7 8 9 10 11 12 ...
2012-01-12 22:31:24 +05:30
read(unit,'(A65536)',end=100) line
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
pos = IO_stringPos(line,maxNchunks)
do i = 1,pos(1)
if (IO_lc(IO_stringValue(line,pos,i)) == 'generate') rangeGeneration = .true.
enddo
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
do l = 1,count
parsing of spectral geom file now accepts two-dimensional microstructure definition (plus range indications by 'x to y') For instance: 1 2 3 4 5 6 7 8 9 10 11 12 ...
2012-01-12 22:31:24 +05:30
read(unit,'(A65536)',end=100) line
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
pos = IO_stringPos(line,maxNchunks)
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
if (verify(IO_stringValue(line,pos,1_pInt),'0123456789') > 0_pInt) then ! a non-int, i.e. set names follow on this line
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
do i = 1,pos(1) ! loop over set names in line
do j = 1,lookupMaxN ! look thru known set names
if (IO_stringValue(line,pos,i) == lookupName(j)) then ! found matching name
first = 2 + IO_continousIntValues(1) ! where to start appending data
last = first + lookupMap(1,j) - 1 ! up to where to append data
IO_continousIntValues(first:last) = lookupMap(2:1+lookupMap(1,j),j) ! add resp. entity list
IO_continousIntValues(1) = IO_continousIntValues(1) + lookupMap(1,j) ! count them
endif
enddo
enddo
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
else if (rangeGeneration) then ! range generation
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
do i = IO_intValue(line,pos,1_pInt),IO_intValue(line,pos,2_pInt),max(1_pInt,IO_intValue(line,pos,3_pInt))
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
IO_continousIntValues(1) = IO_continousIntValues(1) + 1
IO_continousIntValues(1+IO_continousIntValues(1)) = i
enddo
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
else ! read individual elem nums
do i = 1,pos(1)
! write(*,*)'IO_CIV-int',IO_intValue(line,pos,i)
IO_continousIntValues(1) = IO_continousIntValues(1) + 1
IO_continousIntValues(1+IO_continousIntValues(1)) = IO_intValue(line,pos,i)
enddo
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
endif
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 20:08:22 +05:30
enddo
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
endselect
polishing !incremental update, wait for commit of damask_spectral.f90 before checking out
2011-10-18 14:52:33 +05:30
100 endfunction
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 20:08:22 +05:30
# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
2009-10-12 21:31:49 +05:30
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
!********************************************************************
! write error statements to standard out
! and terminate the Marc run with exit #9xxx
! in ABAQUS either time step is reduced or execution terminated
!********************************************************************
removed deadlock when calling IO_warning(33) polished IO_error and IO_warning checked OpenMP critical statements in DAMASK_spectral.f90 corrected writing of headed in DAMASK_spectral.f90
2011-11-02 20:08:42 +05:30
subroutine IO_error(error_ID,e,i,g,ext_msg)
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
CPFEM: time integration scheme changed to start from guessing of Lp. Newton-Raphson algorithm now accelerates as long as residuum improves. IO: added error 650 for polar decomposition problems prec: some cleanup. Introduced "relevantStrain" for convergence checking in CPFEM
2008-02-15 15:34:49 +05:30
use prec, only: pInt
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
implicit none
removed deadlock when calling IO_warning(33) polished IO_error and IO_warning checked OpenMP critical statements in DAMASK_spectral.f90 corrected writing of headed in DAMASK_spectral.f90
2011-11-02 20:08:42 +05:30
integer(pInt), intent(in) :: error_ID
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
integer(pInt), optional, intent(in) :: e,i,g
character(len=*), optional, intent(in) :: ext_msg
removed deadlock when calling IO_warning(33) polished IO_error and IO_warning checked OpenMP critical statements in DAMASK_spectral.f90 corrected writing of headed in DAMASK_spectral.f90
2011-11-02 20:08:42 +05:30
character(len=1024) msg
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
removed deadlock when calling IO_warning(33) polished IO_error and IO_warning checked OpenMP critical statements in DAMASK_spectral.f90 corrected writing of headed in DAMASK_spectral.f90
2011-11-02 20:08:42 +05:30
select case (error_ID)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
!* file handling errors
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (100_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
msg = 'could not open file:'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (101_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
msg = 'write error for file:'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (102_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
msg = 'could not read file:'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (103_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
msg = 'could not assemble input files'
!* material error messages and related messages in mesh
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (150_pInt)
introduced alternative location for material configuration. If present, the >>>JOB.materialConfig<<< file takes precedence over the standard >>>material.config<<<
2011-08-02 15:44:16 +05:30
msg = 'crystallite index out of bounds'
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
case (151_pInt)
introduced alternative location for material configuration. If present, the >>>JOB.materialConfig<<< file takes precedence over the standard >>>material.config<<<
2011-08-02 15:44:16 +05:30
msg = 'phase index out of bounds'
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
case (152_pInt)
introduced alternative location for material configuration. If present, the >>>JOB.materialConfig<<< file takes precedence over the standard >>>material.config<<<
2011-08-02 15:44:16 +05:30
msg = 'texture index out of bounds'
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
case (153_pInt)
introduced alternative location for material configuration. If present, the >>>JOB.materialConfig<<< file takes precedence over the standard >>>material.config<<<
2011-08-02 15:44:16 +05:30
msg = 'sum of phase fractions differs from 1'
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
case (154_pInt)
msg = 'homogenization index out of bounds'
case (155_pInt)
msg = 'microstructure index out of bounds'
case (156_pInt)
msg = 'reading from ODF file'
case (160_pInt)
msg = 'no entries in config part'
case (170_pInt)
msg = 'no homogenization specified via State Variable 2'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (180_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
msg = 'no microstructure specified via State Variable 3'
!* constitutive error messages
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (200_pInt)
introduced alternative location for material configuration. If present, the >>>JOB.materialConfig<<< file takes precedence over the standard >>>material.config<<<
2011-08-02 15:44:16 +05:30
msg = 'unknown constitution specified'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (205_pInt)
introduced alternative location for material configuration. If present, the >>>JOB.materialConfig<<< file takes precedence over the standard >>>material.config<<<
2011-08-02 15:44:16 +05:30
msg = 'unknown lattice structure encountered'
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
added check for unknown keywords again (in agreement with Philip ;-) indeed 'constitution' and '/nonlocal/' are the only keywords which have to be ignored, since they are parsed elsewhere
2012-02-14 14:52:37 +05:30
case (210_pInt)
msg = 'unknown material parameter for j2 constitutive phase:'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (211_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
msg = 'material parameter for j2 constitutive phase out of bounds:'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (212_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
msg = 'unknown constitutive output for j2 constitution:'
added check for unknown keywords again (in agreement with Philip ;-) indeed 'constitution' and '/nonlocal/' are the only keywords which have to be ignored, since they are parsed elsewhere
2012-02-14 14:52:37 +05:30
case (220_pInt)
msg = 'unknown material parameter for phenopowerlaw constitutive phase:'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (221_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
msg = 'material parameter for phenopowerlaw constitutive phase out of bounds:'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (222_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
msg = 'unknown constitutive output for phenopowerlaw constitution:'
added check for unknown keywords again (in agreement with Philip ;-) indeed 'constitution' and '/nonlocal/' are the only keywords which have to be ignored, since they are parsed elsewhere
2012-02-14 14:52:37 +05:30
case (230_pInt)
msg = 'unknown material parameter for titanmod constitutive phase:'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (231_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
msg = 'material parameter for titanmod constitutive phase out of bounds:'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (232_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
msg = 'unknown constitutive output for titanmod constitution:'
added check for unknown keywords again (in agreement with Philip ;-) indeed 'constitution' and '/nonlocal/' are the only keywords which have to be ignored, since they are parsed elsewhere
2012-02-14 14:52:37 +05:30
case (240_pInt)
msg = 'unknown material parameter for dislotwin constitutive phase:'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (241_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
msg = 'material parameter for dislotwin constitutive phase out of bounds:'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (242_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
msg = 'unknown constitutive output for dislotwin constitution:'
case (243_pInt)
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.) I suggest to kill similar logic from other files if present somewhere...
2012-02-14 05:00:59 +05:30
msg = 'zero stacking fault energy'
introduced a flexibility in cut-backing scheme in homogenization.f90 and in crystallite.f90: (1) subStepSizeHomog and subStepSizeCryst := size of substep when cut-back is applied (initially was hard-coded). (2) stepIncreaseHomog and stepIncreaseCryst := step increase when calculation for substep converge (was also hardcoded). introduced a possibility to choose different finite difference scheme, i.e., forward-, backward- and central-difference, for computing grain numerical tangent. note that central-difference scheme will slow down the computation significantly. please use it only if necessary. parameters to set these new features have been included in numerics.f90 and numerics.config, whereas corresponding error messages have been introduced in the IO.f90
2009-11-10 19:06:27 +05:30
added check for unknown keywords again (in agreement with Philip ;-) indeed 'constitution' and '/nonlocal/' are the only keywords which have to be ignored, since they are parsed elsewhere
2012-02-14 14:52:37 +05:30
case (250_pInt)
msg = 'unknown material parameter for nonlocal constitutive phase:'
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
case (251_pInt)
msg = 'material parameter for nonlocal constitutive phase out of bounds:'
case (252_pInt)
msg = 'unknown constitutive output for nonlocal constitution:'
case (253_pInt)
msg = 'element type not supported for nonlocal constitution'
fixed a potential memory leak for hexagonal structures. added some status output to constitutive_xx
2009-10-16 01:32:52 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
!* numerics error messages
fixed a potential memory leak for hexagonal structures. added some status output to constitutive_xx
2009-10-16 01:32:52 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
case (300_pInt)
msg = 'unknown numerics parameter:'
case (301_pInt)
msg = 'numerics parameter out of bounds:'
!* math errors
case (400_pInt)
msg = 'matrix inversion error'
case (401_pInt)
debugged zoo of rotation operations and transformations all angles now in radians introduced a rudimentary check in math_init to complain (IO_error) about broken transformations (e.g. quat --> R --> quat)
2010-05-06 19:37:21 +05:30
msg = 'math_check: quat -> axisAngle -> quat failed'
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
case (402_pInt)
debugged zoo of rotation operations and transformations all angles now in radians introduced a rudimentary check in math_init to complain (IO_error) about broken transformations (e.g. quat --> R --> quat)
2010-05-06 19:37:21 +05:30
msg = 'math_check: quat -> R -> quat failed'
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
case (403_pInt)
debugged zoo of rotation operations and transformations all angles now in radians introduced a rudimentary check in math_init to complain (IO_error) about broken transformations (e.g. quat --> R --> quat)
2010-05-06 19:37:21 +05:30
msg = 'math_check: quat -> euler -> quat failed'
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
case (404_pInt)
debugged zoo of rotation operations and transformations all angles now in radians introduced a rudimentary check in math_init to complain (IO_error) about broken transformations (e.g. quat --> R --> quat)
2010-05-06 19:37:21 +05:30
msg = 'math_check: R -> euler -> R failed'
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
case (405_pInt)
msg = 'I_TO_HALTON-error: An input base BASE is <= 1'
case (406_pInt)
msg = 'Prime-error: N must be between 0 and PRIME_MAX'
case (450_pInt)
msg = 'unknown symmetry type specified'
!* homogenization errors
case (500_pInt)
msg = 'unknown homogenization specified'
debugged zoo of rotation operations and transformations all angles now in radians introduced a rudimentary check in math_init to complain (IO_error) about broken transformations (e.g. quat --> R --> quat)
2010-05-06 19:37:21 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
!* DAMASK_marc errors
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (700_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
msg = 'invalid materialpoint result requested'
!* errors related to spectral solver
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
small polishing, exchanged stops in math.f90 with calls to IO_error
2012-02-13 19:38:07 +05:30
case (802_pInt)
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
msg = 'non-positive dimension'
case (803_pInt)
msg = 'odd resolution given'
case (808_pInt)
msg = 'precision not suitable for FFTW'
case (809_pInt)
msg = 'initializing FFTW'
case (831_pInt)
msg = 'mask consistency violated in spectral loadcase'
case (832_pInt)
msg = 'ill-defined L (each line should be either fully or not at all defined) in spectral loadcase'
case (834_pInt)
msg = 'negative time increment in spectral loadcase'
case (835_pInt)
msg = 'non-positive increments in spectral loadcase'
case (836_pInt)
msg = 'non-positive result frequency in spectral loadcase'
case (837_pInt)
msg = 'incomplete loadcase'
case (838_pInt)
msg = 'mixed boundary conditions allow rotation'
case (842_pInt)
msg = 'missing header info in spectral mesh'
case (843_pInt)
msg = 'resolution in spectral mesh'
case (844_pInt)
msg = 'dimension in spectral mesh'
case (845_pInt)
msg = 'incomplete information in spectral mesh header'
case (846_pInt)
msg = 'not a rotation defined for loadcase rotation'
case (847_pInt)
msg = 'updating of gamma operator not possible if it is pre calculated'
case (880_pInt)
msg = 'mismatch of microstructure count and a*b*c in geom file'
!* Error messages related to parsing of Abaqus input file
extended IO to cope with different name for solverJob and Model polishing, added error codes added FFTW library files
2011-02-21 20:07:38 +05:30
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (900_pInt)
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
msg = 'PARSE ERROR: Improper definition of nodes in input file (Nnodes < 2)'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (901_pInt)
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
msg = 'PARSE ERROR: No Elements defined in input file (Nelems = 0)'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (902_pInt)
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
msg = 'PARSE ERROR: No Element sets defined in input file (Atleast one *Elset must exist)'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (903_pInt)
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
msg = 'PARSE ERROR: No Materials defined in input file (Look into section assigments)'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (904_pInt)
introduced alternative location for material configuration. If present, the >>>JOB.materialConfig<<< file takes precedence over the standard >>>material.config<<<
2011-08-02 15:44:16 +05:30
msg = 'PARSE ERROR: No elements could be assigned for Elset: '
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (905_pInt)
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
msg = 'PARSE ERROR: Error in mesh_abaqus_map_materials'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (906_pInt)
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
msg = 'PARSE ERROR: Error in mesh_abaqus_count_cpElements'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (907_pInt)
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
msg = 'PARSE ERROR: Incorrect size of mesh_mapFEtoCPelem in mesh_abaqus_map_elements; Size cannot be zero'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (908_pInt)
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
msg = 'PARSE ERROR: Incorrect size of mesh_mapFEtoCPnode in mesh_abaqus_map_nodes; Size cannot be zero'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (909_pInt)
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
msg = 'PARSE ERROR: Incorrect size of mesh_node in mesh_abaqus_build_nodes; must be equal to mesh_Nnodes'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (910_pInt)
introduced alternative location for material configuration. If present, the >>>JOB.materialConfig<<< file takes precedence over the standard >>>material.config<<<
2011-08-02 15:44:16 +05:30
msg = 'PARSE ERROR: Incorrect element type mapping in '
adopted improvements done by Arun Prakash. mesh: elemType identification based on lower case Abaqus now reports more errors IO: new function to inquire whether inputfile contains "parts" new function to assemble multiply included inputfile into a flat one awareness of range generation in element numbers error reporting
2010-07-13 15:56:07 +05:30
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
!* general error messages
case (666_pInt)
msg = 'memory leak detected'
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
case default
added subroutine to throw warnings (instead of terminal errors)
2009-03-31 14:51:57 +05:30
msg = 'Unknown error number...'
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
end select
changed Line Ending to Unix style (LF)
2009-01-20 00:40:58 +05:30
added code lines for OpenMP, however it does not work yet!
2008-05-26 18:41:25 +05:30
!$OMP CRITICAL (write2out)
added subroutine to throw warnings (instead of terminal errors)
2009-03-31 14:51:57 +05:30
write(6,*)
removed deadlock when calling IO_warning(33) polished IO_error and IO_warning checked OpenMP critical statements in DAMASK_spectral.f90 corrected writing of headed in DAMASK_spectral.f90
2011-11-02 20:08:42 +05:30
write(6,'(a38)') '+------------------------------------+'
write(6,'(a38)') '+ error +'
write(6,'(a17,i3,a18)') '+ ',error_ID,' +'
write(6,'(a38)') '+ +'
write(6,'(a2,a)') '+ ', trim(msg)
if (present(ext_msg)) write(6,'(a2,a)') '+ ', trim(ext_msg)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
if (present(e)) then
if (present(i) .and. present(g)) then
+ error box + now tells error number (useful for Abaqus)
2011-05-13 22:25:13 +05:30
write(6,'(a13,i6,a4,i2,a7,i4,a2)') '+ at element ',e,' IP ',i,' grain ',g,' +'
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
else
+ error box + now tells error number (useful for Abaqus)
2011-05-13 22:25:13 +05:30
write(6,'(a18,i6,a14)') '+ at ',e,' +'
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
endif
endif
restructured stress integration. Double loop for stress integration (Lp and stress) with cutback step in IntegrateStress is replaced by single loop for Lp without cutback. loop counter and limits are accordingly renamed and implemented
2009-05-28 22:08:40 +05:30
write(6,'(a38)') '+------------------------------------+'
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
call flush(6)
removed implicit type castings
2012-02-10 16:54:53 +05:30
call quit(9000_pInt+error_ID)
small polishing
2010-02-18 15:42:45 +05:30
!$OMP END CRITICAL (write2out)
changed Line Ending to Unix style (LF)
2009-01-20 00:40:58 +05:30
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
! ABAQUS returns in some cases
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endsubroutine
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
added subroutine to throw warnings (instead of terminal errors)
2009-03-31 14:51:57 +05:30
!********************************************************************
! write warning statements to standard out
!********************************************************************
removed deadlock when calling IO_warning(33) polished IO_error and IO_warning checked OpenMP critical statements in DAMASK_spectral.f90 corrected writing of headed in DAMASK_spectral.f90
2011-11-02 20:08:42 +05:30
subroutine IO_warning(warning_ID,e,i,g,ext_msg)
added subroutine to throw warnings (instead of terminal errors)
2009-03-31 14:51:57 +05:30
use prec, only: pInt
implicit none
removed deadlock when calling IO_warning(33) polished IO_error and IO_warning checked OpenMP critical statements in DAMASK_spectral.f90 corrected writing of headed in DAMASK_spectral.f90
2011-11-02 20:08:42 +05:30
integer(pInt), intent(in) :: warning_ID
added subroutine to throw warnings (instead of terminal errors)
2009-03-31 14:51:57 +05:30
integer(pInt), optional, intent(in) :: e,i,g
character(len=*), optional, intent(in) :: ext_msg
removed deadlock when calling IO_warning(33) polished IO_error and IO_warning checked OpenMP critical statements in DAMASK_spectral.f90 corrected writing of headed in DAMASK_spectral.f90
2011-11-02 20:08:42 +05:30
character(len=1024) msg
added subroutine to throw warnings (instead of terminal errors)
2009-03-31 14:51:57 +05:30
removed deadlock when calling IO_warning(33) polished IO_error and IO_warning checked OpenMP critical statements in DAMASK_spectral.f90 corrected writing of headed in DAMASK_spectral.f90
2011-11-02 20:08:42 +05:30
select case (warning_ID)
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (34_pInt)
made DAMASK to work with gfortran: -removed to long lines -restructured f2py modules and merged make_DAMASK2Python into setup processing -setup_code.py now sets library path in makefile and asks for compile switches for spectral code -substituted \ in format strings with $ restructured DAMASK_spectral: -more logical output and structure of code -better input for spectral debug parameters
2011-12-06 22:28:17 +05:30
msg = 'invalid restart increment given'
added _pInt where appropriate
2012-02-02 18:49:02 +05:30
case (35_pInt)
N^2 initialization loop (former spectralPictureMode) rewritten in material.f90 additional output in DAMASK_spectral_interface.f90 132 character cut off in constitutive_nonlocal.f90 rounding error in math.f90 complex number initialization (1.0_pReal)*2.0_pReal*pi new $DAMASK_NUM_THREADS warning in numerics.f90 / IO.f90 polishing in DAMASK_spectral.f90
2012-01-30 19:22:41 +05:30
msg = 'could not get $DAMASK_NUM_THREADS'
restructured algorithm: moved into loop to reallocate fields and replan FFTW in case resolution changes during runtime ==> regridding introduced parameters for selective debugging of spectral code and partly introduced the advanced divergence calculation again which is controlled by debug.config added switch in numerics to control divergence behavior (uncorrected and corrected by phenomenological factor) added precision directive to all values I found
2011-11-15 23:24:18 +05:30
case (47_pInt)
msg = 'No valid parameter for FFTW given, using FFTW_PATIENT'
removed deadlock when calling IO_warning(33) polished IO_error and IO_warning checked OpenMP critical statements in DAMASK_spectral.f90 corrected writing of headed in DAMASK_spectral.f90
2011-11-02 20:08:42 +05:30
case (101_pInt)
lattice: (re)introduced _symmetryType function to replace unsafe lookup array numerics: polishing mpie_cpfem_marc: polishing ..powerlaw: aware of symmetryType function crystallite: aware of symmetryType function, smaller leapfrog acceleration IO: new warning 101 CPFEM: range of odd stress is now -1e15...+1e15, H_sym is used for stiffness
2010-11-03 20:28:11 +05:30
msg = '+ crystallite debugging off... +'
removed deadlock when calling IO_warning(33) polished IO_error and IO_warning checked OpenMP critical statements in DAMASK_spectral.f90 corrected writing of headed in DAMASK_spectral.f90
2011-11-02 20:08:42 +05:30
case (600_pInt)
lattice: (re)introduced _symmetryType function to replace unsafe lookup array numerics: polishing mpie_cpfem_marc: polishing ..powerlaw: aware of symmetryType function crystallite: aware of symmetryType function, smaller leapfrog acceleration IO: new warning 101 CPFEM: range of odd stress is now -1e15...+1e15, H_sym is used for stiffness
2010-11-03 20:28:11 +05:30
msg = '+ crystallite responds elastically +'
removed deadlock when calling IO_warning(33) polished IO_error and IO_warning checked OpenMP critical statements in DAMASK_spectral.f90 corrected writing of headed in DAMASK_spectral.f90
2011-11-02 20:08:42 +05:30
case (601_pInt)
* cpfem will now produce a warning if all components of the stiffness tensor are close to zero
2010-11-04 23:48:01 +05:30
msg = '+ stiffness close to zero +'
removed deadlock when calling IO_warning(33) polished IO_error and IO_warning checked OpenMP critical statements in DAMASK_spectral.f90 corrected writing of headed in DAMASK_spectral.f90
2011-11-02 20:08:42 +05:30
case (650_pInt)
lattice: (re)introduced _symmetryType function to replace unsafe lookup array numerics: polishing mpie_cpfem_marc: polishing ..powerlaw: aware of symmetryType function crystallite: aware of symmetryType function, smaller leapfrog acceleration IO: new warning 101 CPFEM: range of odd stress is now -1e15...+1e15, H_sym is used for stiffness
2010-11-03 20:28:11 +05:30
msg = '+ polar decomposition failed +'
removed deadlock when calling IO_warning(33) polished IO_error and IO_warning checked OpenMP critical statements in DAMASK_spectral.f90 corrected writing of headed in DAMASK_spectral.f90
2011-11-02 20:08:42 +05:30
case (700_pInt)
constitutive_nonlocal: - corrected flux term - multiplication is now aware of dislocation type - corrected change rate for "dipole size" dupper - corrected term for dipole dissociation by stress change - added transmissivity term in fluxes which accounts for misorientation between two neighboring grains (yet hardcoded transmissivity according to misorientation angle) - added more output variables constitutive: - 2 additional variables "previousDotState" and "previousDotState2", which are used to store the previous and second previous dotState (used in crystallite for acceleration/stabilization of state integration) - timer for dotState now measures the time for calls to constitutive_ collectState (used to reside in crystallite_updateState, which is not critical in terms of calculation time anymore) crystallite: - convergence check for nonlocal elments is now done at end of crystallite loop, not at the beginning; we simple set all elements to not converged if there is at least one nonlocal element that did not converge - need call to microstructure before first call to collect dotState for dependent states - stiffness calculation (jacobian): if there are nonlocal elements, we also have to consider changes in our neighborhood's states; so for every perturbed component in a single ip, we have to loop over all elements; since this is extremely time-consuming, we just perturb one component per cycle, starting with the one that changes the most during regular time step. - updateState gets a damping prefactor for our dotState that helps to improve convergence; prefactor is calculated according to change of dotState IO: - additional warning message for unknown crystal symmetry
2009-12-15 13:50:31 +05:30
msg = '+ unknown crystal symmetry +'
added subroutine to throw warnings (instead of terminal errors)
2009-03-31 14:51:57 +05:30
case default
lattice: (re)introduced _symmetryType function to replace unsafe lookup array numerics: polishing mpie_cpfem_marc: polishing ..powerlaw: aware of symmetryType function crystallite: aware of symmetryType function, smaller leapfrog acceleration IO: new warning 101 CPFEM: range of odd stress is now -1e15...+1e15, H_sym is used for stiffness
2010-11-03 20:28:11 +05:30
msg = '+ unknown warning number... +'
added subroutine to throw warnings (instead of terminal errors)
2009-03-31 14:51:57 +05:30
end select
!$OMP CRITICAL (write2out)
write(6,*)
removed deadlock when calling IO_warning(33) polished IO_error and IO_warning checked OpenMP critical statements in DAMASK_spectral.f90 corrected writing of headed in DAMASK_spectral.f90
2011-11-02 20:08:42 +05:30
write(6,'(a38)') '+------------------------------------+'
write(6,'(a38)') '+ warning +'
write(6,'(a38)') '+ +'
write(6,'(a17,i3,a18)') '+ ',warning_ID,' +'
write(6,'(a2,a)') '+ ', trim(msg)
if (present(ext_msg)) write(6,'(a2,a)') '+ ', trim(ext_msg)
added subroutine to throw warnings (instead of terminal errors)
2009-03-31 14:51:57 +05:30
if (present(e)) then
* cpfem will now produce a warning if all components of the stiffness tensor are close to zero
2010-11-04 23:48:01 +05:30
if (present(i)) then
if (present(g)) then
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
write(6,'(a12,1x,i6,1x,a2,1x,i2,1x,a5,1x,i4,a2)') '+ at element',e,'IP',i,'grain',g,' +'
* cpfem will now produce a warning if all components of the stiffness tensor are close to zero
2010-11-04 23:48:01 +05:30
else
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
write(6,'(a12,1x,i6,1x,a2,1x,i2,a13)') '+ at element',e,'IP',i,' +'
* cpfem will now produce a warning if all components of the stiffness tensor are close to zero
2010-11-04 23:48:01 +05:30
endif
added subroutine to throw warnings (instead of terminal errors)
2009-03-31 14:51:57 +05:30
else
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
write(6,'(a12,1x,i6,a19)') '+ at element',e,' +'
added subroutine to throw warnings (instead of terminal errors)
2009-03-31 14:51:57 +05:30
endif
endif
restructured stress integration. Double loop for stress integration (Lp and stress) with cutback step in IntegrateStress is replaced by single loop for Lp without cutback. loop counter and limits are accordingly renamed and implemented
2009-05-28 22:08:40 +05:30
write(6,'(a38)') '+------------------------------------+'
constitutive_nonlocal: - read in activation energy for dislocation glide from material.config - changed naming of dDipMin/Max to dLower/dUpper - added new outputs: rho_dot, rho_dot_dip, rho_dot_gen, rho_dot_sgl2dip, rho_dot_dip2sgl, rho_dot_ann_ath, rho_dot_ann_the, rho_dot_flux, d_upper_edge, d_upper_screw, d_upper_dot_edge, d_upper_dot_screw - poisson's ratio is now calculated from elastic constants - microstrucutre has state as first argument, since this is our output variable - periodic boundary conditions are taken into account for fluxes and internal stresses. for the moment, flag has to be set in constitutive_nonlocal. - corrected calculation for dipole formation by glide - added terms for dipole formation/annihilation by stress decrease/increase constitutive: - passing of arguments is adapted for constitutive_nonlocal model crystallite: - in stiffness calculation: call to collect_dotState used wrong arguments - crystallite_postResults uses own Tstar_v and temperature, no need for passing them from materialpoint_postResults homogenization: - crystallite_postResults uses own Tstar_v and temperature, no need for passing them from materialpoint_postResults IO: - changed error message 229 material.config: - changed example for nonlocal constitution according to constitutive_nonlocal all: - added some flush statements
2009-10-20 20:06:03 +05:30
call flush(6)
corrected two OMP CRITICAL terminations
2010-04-06 12:17:15 +05:30
!$OMP END CRITICAL (write2out)
added subroutine to throw warnings (instead of terminal errors)
2009-03-31 14:51:57 +05:30
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
endsubroutine
added subroutine to throw warnings (instead of terminal errors)
2009-03-31 14:51:57 +05:30
let's see whether this is the thing to do..!?!?
2007-03-20 19:25:22 +05:30
END MODULE IO