Commit Graph

2182 Commits

Author SHA1 Message Date
Philip Eisenlohr 1df9be233d changed default threshold to 0.0 2012-11-08 17:08:41 +00:00
Philip Eisenlohr 682a98b704 forgot to alter name of 'geom_fromAng' in setup script... 2012-11-08 16:29:26 +00:00
Philip Eisenlohr 025d6c9048 renamed script to new convention
aligned output behavior to that of fromVoronoiTessellation, i.e., standard output is geom, --config gives associated material.config stud
2012-11-08 15:44:51 +00:00
Philip Eisenlohr af2aa8500d added option to specify crystallite index 2012-11-08 15:43:38 +00:00
Christoph Kords c775edaa6d better avoid initialization of residuum_old with huge, since it posed problems in Abaqus 2012-11-08 13:26:22 +00:00
Martin Diehl 07fd2681c6 slight changes in results required update of reference 2012-11-08 10:22:54 +00:00
Martin Diehl 3775e478cc added tests for gfortran (will crash, but might help to find the bug) 2012-11-08 09:52:24 +00:00
Martin Diehl f5c6b4e0a6 added extra compile debug options for more detailed information on crash 2012-11-08 09:49:20 +00:00
Philip Eisenlohr c4c27860b2 renamed script
updated to recent standard in io processing
introduced "primitive" bicontinuous structure in addition to "diamond" and "gyroid"
2012-11-07 16:52:10 +00:00
Philip Eisenlohr 7cf2024be3 fixed buggy parameter checking 2012-11-07 16:51:17 +00:00
Philip Eisenlohr e11d107956 fixed reporting to STDOUT/STDERR
fixed buggy parameter checking
2012-11-07 16:49:47 +00:00
Christoph Kords dad9922f54 fixed bug in crystallite_FPI: stateDamper always has to be defined for each grain
added some OMP FLUSH statements were necessary 
replaced openmp do by forall construct where possible; this is much safer and perhaps even as fast for small loops
2012-11-07 15:43:29 +00:00
Martin Diehl bbcffa668b some small changes in including PETSc *.h90 files to prevent warnings.
It seems that a bug in old gfortran 4.5.0 causes the problems, at least with 4.7.3 it seems to work after fixing the bugs in IO.f90
2012-11-07 13:11:41 +00:00
Martin Diehl 9e20149c6f updated reference results, differ because changes in crystallite (matrix inversion removed) 2012-11-07 09:57:34 +00:00
Martin Diehl 9b99825ac6 preventing array index out of bounds in case of empty line, removed to goto statements 2012-11-07 09:31:46 +00:00
Martin Diehl 5214b93342 fixed bug when reading in geometry for spectral solver 2012-11-06 17:16:01 +00:00
Martin Diehl 3ada4897fb reordered PETSc initialization (now first part done in interface, setting of parameters done in numerics), removed unnecessary includes for PETSc
still not running with gfortran, use at own risk!
2012-11-06 16:00:51 +00:00
Martin Diehl a86d528a4a replaced call flush(6) by flush(6), slightly improved formating for output of gfortran version 2012-11-06 15:50:20 +00:00
Philip Eisenlohr e0dc2710b5 now works in "sourcing" mode.
still have to worry about figuring out the lib/pathinfo details and setting LD_LIBRARY_PATH accordingly...
2012-11-06 15:40:45 +00:00
Christoph Kords 5b6baa7c0d introduced distinction between ip coordinates and cell center coordinates; the former can be globally calculated by subroutine mesh_build_ipCoordinates, the latter locally by the function mesh_cellCenterCoordinates; renamed mesh_ipCenterOfGravity to mesh_ipCordinates 2012-11-06 14:37:13 +00:00
Christoph Kords bb033c5fe7 Simplified algorithm of crystallite_integrateStress while preserving (almost) same functionality.
Removed "leapfrogging" (increase of step for next guess, when last guess was ok); Replaced Armijo rule testing for step size by simple check if the residuum got better, since the former virtually did not have any effect; consistently using the 2-norm of the residuum rather than infinity-norm for the convergence check throughout the function
2012-11-06 12:35:45 +00:00
Martin Diehl 1583ae74c3 adopted new naming scheme 2012-11-06 08:49:48 +00:00
Philip Eisenlohr cfc8823328 missed those three in former commit... 2012-11-05 21:28:01 +00:00
Philip Eisenlohr 1bd4262cfb renaming of (loosely related) "spectral" scripts to (more tightly related) "geom" versions. 2012-11-05 21:19:12 +00:00
Christoph Kords 639d6e0655 inserted missing flush(crystallite_todo) in parallel section; might have caused trouble in parallel calculations before 2012-11-05 14:45:52 +00:00
Martin Diehl 1fc59107fa removed debug statement in averageDown.py, deleted obsolete spectral_iterationCount.py (new solver will report statistics anyway and most probably the current count does not work as the screen output was modified) 2012-10-31 10:06:11 +00:00
Martin Diehl 689b64e3f5 updated geom pack/unpack test to new definition 2012-10-31 10:04:53 +00:00
Martin Diehl 086fe138b1 substituted matrix inversion to solve equation by direct solution routine from LAPACK 2012-10-31 09:56:26 +00:00
Philip Eisenlohr c0b83bd554 exported PDF versions as 300dpi PNG 2012-10-31 09:28:46 +00:00
Philip Eisenlohr afaf1f4e8c made errors report to STDERR.
maxMicrostructure remains silent if not present.
2012-10-30 16:32:00 +00:00
Philip Eisenlohr f2229883b5 exchanged "maxGrainCount" for "maxMicrostructure", since this is what it really means... 2012-10-30 16:15:10 +00:00
Philip Eisenlohr ff164c9209 changed packing syntax from "N copies of x" to "N of x" (now consistent with IO.f90) 2012-10-30 15:53:46 +00:00
Philip Eisenlohr 2d36d52cf7 changed packing syntax from "N copies of x" to "N of x" (now consistent with geomPack and IO.f90) 2012-10-30 15:47:11 +00:00
Test User be5d072e2d Added report of SUCCESSFUL test for Revision 1848 2012-10-30 03:59:17 +00:00
Christoph Kords 55dc286be2 corrected limit for upper stable dipole distance; fixed a bug in the postResults where division by zero might happen 2012-10-29 13:02:01 +00:00
Christoph Kords f666f8dcf3 added possibility to scale the plastic shearrate by the ratio of mobile to total density; keyword in material.config is "deadZone" 2012-10-29 12:49:28 +00:00
Martin Diehl b9aa50c59a improved performance for large packing values 2012-10-26 12:59:43 +00:00
Martin Diehl 4e5ae4056b did a mistake during last update of this test, now it should run again 2012-10-25 13:16:48 +00:00
Martin Diehl 78aa010565 now just comparing substring to figure out if setup_processing was importing the core module 2012-10-25 12:40:06 +00:00
Martin Diehl ca596bc7af error messages are now printed, added svn properties 2012-10-25 09:16:17 +00:00
Martin Diehl 339205c890 added missing quotes 2012-10-25 07:56:58 +00:00
Martin Diehl 868984cb89 import of damask.core will not complain any more if done in setup_processing.py, however setup_processing.py will check if core.so exist after compilation and delete it before compilation 2012-10-25 07:55:27 +00:00
Claudio Zambaldi b212da5928 this version of the env script maybe works 2012-10-24 15:22:30 +00:00
Claudio Zambaldi d6f69568ad Dynamically create DAMASK_ROOT and DAMASK_BIN
by sourcing this script.

also documents the changes to .bashrc if one wants to 
'install' DAMASK permanently
2012-10-24 15:19:41 +00:00
Christoph Kords 93cc466749 third entry in mesh_ipNeighborhood stores the neighbor_index that points from each neighbor back to the central ip; needed in nonlocal model 2012-10-24 14:03:02 +00:00
Martin Diehl 13b55275b1 documented utilities and structured, worked on the restart capabilities of the new basic solver 2012-10-24 11:31:40 +00:00
Martin Diehl aefe8d7e32 added a test for compiling all spectral solver variants with gfortran 2012-10-24 09:26:21 +00:00
Martin Diehl 265a03d32a new references, as the results are slightly different for corrected phenopowerlaw 2012-10-23 16:25:36 +00:00
Philip Eisenlohr 5ad0eda1b6 switched saturation behavior!!
was \dot s_alpha = (1-s_alpha/s_inf)^a h_alphabeta \dot gamma_beta
now \dot s_alpha = h_alphabeta  (1-s_beta/s_inf)^a \dot gamma_beta

current form is consistent with the genmat implementation (and appears to make more physical sense). Kalidindi_etal1992 suggested this form, but altered it to the alpha-one in Bachu+Kalidindi1998... By now, it seems that some groups use alpha, others beta approach.

introduced two new absolute tolerance values for "shears" and "twinFrac" (default 1e-6).
2012-10-22 14:55:07 +00:00
Christoph Kords 2be331b74d In crystallite for state integrators 2 and 3: now that "stateJump" is called before integrating dotState, we have to do state = state + dotState *dt instead of state = subState0 + dotState * dt; otherwise the deltaState is lost 2012-10-22 13:34:15 +00:00