Martin Diehl
724ec040a2
added pInts and pReals
2012-02-09 15:55:04 +00:00
Philip Eisenlohr
3228cf563c
fixed bug when column was referenced multiple time in a formula.
...
changed column tagging character to '#' ('$' gave some trouble on the shell...)
2012-02-09 15:46:25 +00:00
Martin Diehl
f8013ce4bf
test now running, but proof not done yet
2012-02-09 14:42:59 +00:00
Martin Diehl
7960600dbe
improved updating of reference results
2012-02-09 14:41:18 +00:00
Martin Diehl
a329a9b91d
added flags for ifort and gfortran to use preprocessor
2012-02-09 13:30:59 +00:00
Philip Eisenlohr
dced3f9430
changed FFT reporting of dimension and resolution to one line each.
2012-02-09 13:18:59 +00:00
Philip Eisenlohr
0f7c74e02d
set ID flag and line endings
2012-02-09 13:06:22 +00:00
Martin Diehl
195491aaaa
added missing pInts
2012-02-09 12:35:55 +00:00
Martin Diehl
a26456f9c0
now handling capitals in restart argument correctly
2012-02-09 12:20:14 +00:00
Martin Diehl
d2bb7c10e5
new test for addCurl
2012-02-09 08:58:45 +00:00
Martin Diehl
e909efdad6
removed accidently inserted comments
2012-02-09 08:54:09 +00:00
Claudio Zambaldi
aad69a6b7c
with the min # recycles set to 1 the compile-check (1101) runs now
...
started to get hex-elastic (2101) up and running again with new test class
postprocessing to be done for both
2012-02-08 16:38:29 +00:00
Claudio Zambaldi
eec2af6505
set min number of recycles to 1 for one-element Marc calculations
2012-02-08 15:19:47 +00:00
Philip Eisenlohr
5a6cd1fc8f
slightly prettified handling of finite difference accuracy requesting
2012-02-08 12:01:38 +00:00
Philip Eisenlohr
56ce779f57
generalized slip system idea of Christoph to field data of "special" dimension (i.e. not 3 or 9).
...
dropped "norm"-suffix from input names of norms.
ASCIItable output is now called "normType" with type being Abs, Frobenius, or Max...
2012-02-08 12:00:31 +00:00
Philip Eisenlohr
282e4a0360
scripts now figure out dimension and resolution on their own...
2012-02-08 11:18:15 +00:00
Philip Eisenlohr
9c35e4e148
small polishing of help output.
2012-02-08 09:12:08 +00:00
Christoph Kords
51e8981ba3
not specifying the location of python executable directly, but rather getting it from $path, otherwise the scripts won't work on every machine
2012-02-07 13:09:10 +00:00
Christoph Kords
2f07faa0e2
norm capable of using different p-norms (abs norm, frobenius norm, and max norm). also able to handle "slipsystem" outputs (use option -s). number of slip systems can be set by option -i, if not 12 is default.
2012-02-07 13:02:56 +00:00
Christoph Kords
e7d407c2f6
crystallite_Fp has to be initialized in crystallite_init, cause it is used by constitutive_microstructure
2012-02-03 14:52:09 +00:00
Christoph Kords
5ddade27fc
changed Peierls stress to more realistic values for fcc
2012-02-03 13:12:45 +00:00
Christoph Kords
d62eddc0cd
corrected kinetics law and changed parameters. For solid solution hardening there are 3 parameters: the activation energy, the concentration of obstacles that determines the activation length and meanfreepath, and the obstacle size that determines the activation volume. For the Peierls mechanism there is: the width of doublekinks that determines the activation volume and the Peierls stress for edge and screw.
...
Still testing needed to check whether the current formulation makes sense or not.
2012-02-03 12:50:54 +00:00
Christoph Kords
bbf4f25898
renamed subroutine "mesh_marc_get_mpieOptions" to "mesh_get_damaskOptions" which is now also used for abaqus and spectral input files. The only available option so far is "periodic", which allows for periodic fluxes in the specified direction. Examples for usage (add the following lines to the respective input file, in case of spectral file this goes into the header):
...
- periodicity in x and z direction for marc:
$damask periodic x z
- periodicity in y direction for abaqus:
**damask periodic y
- periodicity in x and y direction for spectral:
periodic y x
2012-02-03 12:37:52 +00:00
Philip Eisenlohr
a3e47b2d65
added method to return data as list of floats.
...
safer check whether ASCII file has ended.
2012-02-02 17:13:51 +00:00
Philip Eisenlohr
bffc22fbe1
major restructuring. packing stencil can be shifted to allow for element or nodal value averaging.
2012-02-02 17:12:48 +00:00
Krishna Komerla
7e23f84d2c
new subroutine find_nearest_neigbor, basically a wrapper function for our purpose to kdtree2
...
split up kdtree2.f90 again from math.f90 to make the files better accessible
2012-02-02 15:57:22 +00:00
Martin Diehl
72aeb2113a
minor improvement on the still not working long (polarization field based) algorithm
2012-02-02 13:20:09 +00:00
Martin Diehl
faed4ec232
added deallocation of arrays (fftw_free)
2012-02-02 13:19:35 +00:00
Martin Diehl
2e4f2f3d9c
added _pInt where appropriate
2012-02-02 13:19:02 +00:00
Martin Diehl
683384681a
doing average correction now in Fourier space, now sure that the constant term is correct.
...
changed order of stress BC calculation/spectral method to avoid average calculation of stress in real space
2012-02-01 20:30:27 +00:00
Martin Diehl
c6fb2122be
added some comments for further debuggin/compliance with standard checking
2012-02-01 20:28:58 +00:00
Martin Diehl
c5128e94eb
some f2008 standard related correction i didn't check in last time
2012-02-01 20:20:05 +00:00
Claudio Zambaldi
7dd67e3ebb
improved creation of material.config files, mostly done by Philip
2012-02-01 19:44:36 +00:00
Claudio Zambaldi
590b10ee3e
not much
2012-02-01 19:43:11 +00:00
Martin Diehl
800e291240
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-01-31 19:18:55 +00:00
Franz Roters
93f311f2b9
now skips environment variables in comments, i.e. after '#'
2012-01-31 15:43:39 +00:00
Martin Diehl
80583fefb5
added compilation_info.f90, to store compilation date,time and compiler info in each file
...
added #include statement at each init() routine
2012-01-31 14:54:49 +00:00
Martin Diehl
e87296b10a
added debug output for tests, corrected bug reporting number of successful tests
2012-01-31 13:10:14 +00:00
Christoph Kords
9b679c3c0b
switched to damask include instead of outdated damask_tools
2012-01-31 13:08:25 +00:00
Martin Diehl
8cf67e1be6
corrected weighting of difference divergence measures, added some additional output and introduced exit code 0 for successful termination.
2012-01-30 20:25:04 +00:00
Martin Diehl
9464937db7
merge declaration and initialization with standard values
...
corrected handling of $DAMASK_NUM_THREADS
2012-01-30 20:16:19 +00:00
Martin Diehl
eeda357710
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 13:52:41 +00:00
Martin Diehl
7731926fec
two step renaming because otherwise svn get confused
2012-01-27 09:58:08 +00:00
Martin Diehl
2a89d9b81a
simplified compilation testing by introducing some small helper functions for common tasks, continued to rename
2012-01-27 09:55:19 +00:00
Martin Diehl
7c35e5a1cd
renamed to new numbering scheme, added first test cases for addDivergence
2012-01-26 17:29:26 +00:00
Martin Diehl
20a775e8d1
small improvements (more output to screen)
2012-01-26 15:26:21 +00:00
Martin Diehl
c6c9fd71de
marc compile test is compiling and running again, postprocess still not working
2012-01-26 15:23:41 +00:00
Christoph Kords
4edd824d8c
forgot homogenization.f90 in last commit
2012-01-26 13:51:01 +00:00
Christoph Kords
1330576a01
added new math function "math_deviatoric33" which returns the deviatoric part of a 3x3 tensor
...
renamed some math functions, so that we have a universal naming scheme: for matrix multiplications use an "x" (e.g. math_mul33x3); don't use the "x" to describe the shape of the tensor that the function is applied to (e.g. math_invert33 instead of math_invert3x3)
2012-01-26 13:50:00 +00:00
Christoph Kords
df931890e0
corrected typos
2012-01-26 12:50:04 +00:00