Martin Diehl
c54600fd1f
added kdtree2 source and changed makefile to compile it.
...
started to implement the nearest neighbor search for regridding
corrected calculation of divergence in real space.
corrected handling of maximum stress deviation (removed mask)
2012-01-04 17:43:26 +00:00
Martin Diehl
ee401520d7
corrected divergence calculation. Still some open questions, but improved understanding a lot
2011-12-23 12:30:35 +00:00
Martin Diehl
6dee2616f5
added help info how to use the executable
2011-12-23 12:29:04 +00:00
Martin Diehl
99fa0e0be8
corrected calculation of curl (last version was more a dummy function)
2011-12-23 12:23:13 +00:00
Philip Eisenlohr
e62d083f7a
removed (virtually) dead single-precision branch of the spectral solver.
2011-12-22 13:19:27 +00:00
Philip Eisenlohr
35bd460e39
added "core" package capability to processing
...
(renamed former f2py modules)
reverted pathinfo...
2011-12-22 10:36:59 +00:00
Krishna Komerla
f463c2da70
added addCurl.py, not sure if it works correctly. still in old (eq. to addDivergence.py) style
2011-12-21 17:25:31 +00:00
Philip Eisenlohr
b0ba5f7ba3
code setup now takes BLAS in order of precedence from IKML, ACML, and LAPACK.
...
Leave pathinfo blank for BLAS-options that are not available (like IKML typically)
2011-12-20 10:58:51 +00:00
Franz Roters
c4a4c5cdca
corrected error in apply_DAMASK_modifications
...
added support for Marc/Mentat2011 (classic)
2011-12-19 16:06:03 +00:00
Philip Eisenlohr
739c524848
makefile substitution now based on regexp match.
...
condensed makefile syntax.
can now use a system-defined LAPACK instead of ACML (required for OS X...).
fixed bug that 'ULTRA' was not using -O3 for most of the compiling...
2011-12-16 09:11:03 +00:00
Philip Eisenlohr
a8dece8694
adopted use of general damask package ("import damask" instead of "damask_tools")
2011-12-15 08:52:41 +00:00
Christoph Kords
b2547e0117
Math inversion used to return zero (math_inv3x3) or error (math_invert3x3) for negative determinant. Now checking whether the absolute(!) value of the determinant is close to zero to avoid singularities, negative determinants are very well allowed.
2011-12-14 08:55:24 +00:00
Martin Diehl
031a00237f
improved makefile, no directly linking to files of acml
...
corrected default pathinfo
improved test for intel compiler
2011-12-09 16:34:53 +00:00
Martin Diehl
7d1d2b6229
corrected small mistake in setup_code and makefile
...
test no. 1500, spectral compile test with intel compiler is now working. adjusted run_test and tamask_test
2011-12-09 15:13:27 +00:00
Philip Eisenlohr
2ceb78c337
introduced new damask_tools method "relPath" and adopted it in setup_processing.py and setup_code.py
2011-12-06 18:38:03 +00:00
Martin Diehl
ea0fe7b406
corrected 3D visualize (math.f90 part was not working)
...
small polishing of voronoi_tessellation.f90
2011-12-06 17:46:33 +00:00
Martin Diehl
5ebeb96e85
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 16:58:17 +00:00
Philip Eisenlohr
efadf9f728
mainly fixed error in output of spectral results (1:N,…) instead of (N,…)
...
rearranged some logic here and there.
(hopefully) improved readability of debug/standard output.
restarting logic would need some discussion with Martin/Krishna still…
2011-12-04 10:01:32 +00:00
Martin Diehl
ace6851389
moved all routines from postprocessingMath to math.90, renamed the module to DAMASK, changed scripts and interfaces accordingly.
...
polished math.f90 (mainly added _pInt/_pReal and intent(in/out))
curl_fft is still a dummy function
2011-12-01 12:01:13 +00:00
Franz Roters
e738bee1e3
moved todolist to private tree
2011-11-28 08:18:07 +00:00
Franz Roters
515986c4d6
renamed w0 and w0_slip to a and a_slip
2011-11-23 14:48:39 +00:00
Philip Eisenlohr
f685bbba0e
exchanged location of "graincount" in user block.
...
now matches wiki entry:
1. homogenization block
2. grain block
(grain count
sub-block[s] of crystallite and constitutive)
2011-11-23 09:09:00 +00:00
Martin Diehl
5a1e73b53d
added set_timelimit functionality to fortran-wrapper of FFTW and introcudec it to DAMASK_spectral.f90
2011-11-21 18:48:38 +00:00
Martin Diehl
566f16b6e9
implemented calculation of divergence in real space, polished spectral debugging
2011-11-21 18:12:40 +00:00
Martin Diehl
dc6c29a910
important bugfix for reading in results in case of restart
2011-11-17 22:11:05 +00:00
Martin Diehl
5ef73e164a
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 17:54:18 +00:00
Krishna Komerla
60c9293baf
restarting seems to work, spectral solver writes own defgrad to disk.
...
step counting rectified
added additional output of deformation gradient volume min and max
2011-11-11 14:17:43 +00:00
Christoph Kords
860594d9a4
forgot to multiply shearrate with timestep in order to get accumulated shear
2011-11-09 10:26:00 +00:00
Christoph Kords
6dc8a4d530
check for LFC only if there is any slip activity at all; also increased safety factor from 1.2 to 1.5
2011-11-09 09:55:39 +00:00
Christoph Kords
7033e84fc7
new outputs "accumulatedshear", "velocity_edge", and "velocity_screw" (the latter two replace the output "dislocationvelocity")
2011-11-09 09:22:52 +00:00
Martin Diehl
de96e99bca
corrected determination of restartRead from commandline (now done by FEsolving instead of DAMASK_spectral_interface)
...
recompiled fftw libraries with --enable-shared to make them work with f2py
added draft for compilation test of spectral solver
2011-11-07 18:25:10 +00:00
Christoph Kords
8bc39fcd48
immobile density output keeps its sign
2011-11-07 15:46:40 +00:00
Christoph Kords
34adec74e8
corrected wrong indices into state array for velocity
2011-11-07 15:17:58 +00:00
Krishna Komerla
49c0b8a3fd
restructured restarting capabilities
2011-11-07 11:04:57 +00:00
Christoph Kords
7dfb96a3da
constitutive_nonlocal:
...
dislocation velocities are stored in the state, so we actually now have three "parts" of the state, the basic states that are updated by "constitutive_dotState" come first, then the dependent states that are calculated by "constitutive_microstructure" follow, and finally we have a last part reserved for other variables that just use the memory reserved by the state array and are updated somewhere else.
constitutive:
LpAndItsTangent does not need the full state, but only the local state, so changed that at least for the nonlocal constitutive law
2011-11-04 13:12:17 +00:00
Christoph Kords
9d1bc584d0
corrected stress integration scheme: now use norm(Lpresiduum) as a target function for Armijo's rule instead of whole tensor Lp; also corrected the guess for the improvement in Lp
2011-11-04 12:57:12 +00:00
Christoph Kords
ca3d21a3b6
Depending on which state integrator one uses for the stiffness calculation, the initial state has to be chosen accordingly: e.g. for FPI choose last converged state, but for explicit RK choose converged state from start of increment (in case of explicit euler no state integration at all, but only stress integration). For this purpose we also need to remember "Fe" which now follows the cutbacking procedure as it is used for "Fp".
2011-11-04 12:44:50 +00:00
Martin Diehl
ee7022d8cf
correctet random seed generation to work with gfortran.
2011-11-04 10:29:35 +00:00
Martin Diehl
c2eac36b48
started to implement restart facilities for spectral solver.
...
restart write is on per default
restart read is switched on by using --restart or -r INT where INT gives step at which the calculation should restart
setting INT to a value <1 will turn restart write off
2011-11-03 19:32:11 +00:00
Philip Eisenlohr
b9647b2584
moved parts into new "lib" folder
2011-11-03 13:21:11 +00:00
Martin Diehl
680ba9082f
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 14:38:42 +00:00
Claudio Zambaldi
447e69019c
Introduced environment variable DAMASK_ROOT.
...
It is set in the DAMASK root folder by running damask_env.sh
Much needed in the testing routines so far (try ./testing/run_tests.py)
damask_env.sh could also trigger other scripts to get to a working setup after fresh checkout. For example it could run
./processing/setup/setup_processing.py
./code/setup/setup_code.py
if it finds that we have a fresh checkout.
2011-10-27 17:55:42 +00:00
Martin Diehl
15c356c3a7
some small improvements regarding the rotation of the loadcase frame
2011-10-25 13:38:24 +00:00
Martin Diehl
c13aa2a829
first try of implementing a rotation of the loadcase coordinate system. Seems to work so far for one loadcase and homogeneous materials.
2011-10-24 18:26:34 +00:00
Martin Diehl
0f34d14bee
again edited the makefile, no there is more flexibility in defining the libraries (FFTW and ACML)
...
added the source of fftw-3.3.tar.gz to fulfill the GPL
set_python_env.py was forgotten during last commit (for testing)
2011-10-24 16:57:51 +00:00
Martin Diehl
8ac1891b3c
moved DAMASK_NaN to prec.f90 respectively prec_single.f90 as it is precision dependent
2011-10-20 16:46:11 +00:00
Martin Diehl
ac342ad682
worked on the updating on increments, does not work properly until now
2011-10-19 09:05:02 +00:00
Krishna Komerla
0a46741a77
changing header to correct number of increments
2011-10-18 16:42:06 +00:00
Martin Diehl
9d3f74ce5f
polishing
2011-10-18 14:46:07 +00:00
Martin Diehl
7746390688
using keywords to indicate geometry and loadcase
...
enabled finetuning of FFTW
added some debugging options
reading in rotation of boundary conditions
using header in geometry file
corrected error in calculating tolerance for stress BC
polishing of output, variable declaration, and variable names
2011-10-18 14:45:32 +00:00