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
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
7642509a47
first working draft of automated code testing, run_tests.py assumes DAMASK is located under ~/DAMASK/
2011-10-24 16:24:15 +00:00
e7a60988a1
modified exit number extraction
2011-10-24 13:02:08 +00:00
93b2160c2f
extract exit number from MSC out-file. Other MSC utility functions can also go into MSC_TOOLS class.
2011-10-24 12:59:46 +00:00
414bacb096
deleted unused vars. changed linking library info/path
2011-10-20 12:47:57 +00:00
ed08dfc2be
added option to select between cell-centered (as before) or vertex-centered data
2011-10-20 12:29:59 +00:00
c78227cec1
added option to skip data points along x,y,z in case they are periodic images. relevant for vertex-centered data for which outer surface data is doubled, i.e. data(0,0,0) == data(Nx,Ny,Nz)
2011-10-20 12:28:39 +00:00
3a0911ca89
added option to output RGB palette on terminal
...
number of colors is now variable (default 32)
warning issued if MSC path not valid
new "redblue" sym colormap
2011-10-20 11:59:15 +00:00
96223ddf39
corrected for-loop index bug. two-dim output is now optional.
2011-10-14 23:43:45 +00:00
708b4d2eca
seed-file header format corrected, setup script now contains minimalSurface and vicinityOffset. removed .exe from compile output
2011-10-14 22:39:58 +00:00
ba0488638b
alter a spectral geom file such that any voxel that sees within a given vicinity a different microstructure than itself gets offset by a given value. this allows to identify/isolate grain boundary regions...
2011-10-11 17:35:53 +00:00
98c9bf9c15
generates spectral geom file for (currently) gyroid or diamond minimal surfaces. Inside and outside of a given threshold value are distinguished by microstructure 1 and 2.
2011-10-11 17:30:55 +00:00
ffc6442098
can now deal with arbitrarily long lines for microstructure definition. useful for 2D ASCII representation of grain morphology in the geom file.
...
additional internal code reworking.
2011-10-11 17:25:22 +00:00
c7c541f393
corrected typos and such
2011-09-13 15:57:58 +00:00
d235de1aa5
added two small (quick and dirty) tools to convert data from EBSD to input files for spectral method, put them together with patchFromReconstructedBoundaries into new folder.
2011-09-13 15:46:44 +00:00
b412239d9f
reflects new output format of DAMASK_spectral to find
...
*** Loadcase...
lines.
2011-09-05 13:39:44 +00:00
48ac3956ce
added '\n' escaping of ID string in VTK header
2011-09-05 13:38:10 +00:00
f458f026db
forgot to include string module, too
2011-09-05 08:54:29 +00:00
0ac5109c94
forgot to include string module
2011-09-05 08:50:36 +00:00
4eb8e25b6d
no fitting recent file renames
2011-08-26 12:25:11 +00:00
f3c2547fba
name change to better reflect what happens: comparison between local deformation field and a (necessarily compatbile) geometry reconstructed from the whole deformation field.
2011-08-26 07:17:35 +00:00
eb0ce3992b
debugging of addDebugInformation (now working, but not tested in detail)
...
Added Id to postprocessingMath
2011-08-25 18:18:38 +00:00
7aba05ed9f
reworked to match former script layouts and logics. (not yet tested, use at your own risk)
2011-08-25 17:55:36 +00:00
282d2d7b97
corrected wrong Id tag
2011-08-25 16:47:12 +00:00
bbee6c4658
forget to check in addDebugInformation
2011-08-25 16:44:36 +00:00
a5f176cf18
put divergence calculation using FDM to postprocessingMath, added divergence calculation using differentation in Fourier space to postprocessingMath
...
added functions to calculate the shape and volume mismatch of compatible and non-compatible geometry reconstruction to postprocessingMath. They can be accessed using addDebugInformation
2011-08-25 15:31:05 +00:00
cc3bc4a216
made handling of unknown dimension (i.e. resolution ==1) the same as for preprocessing (last commit)
2011-08-25 15:28:33 +00:00
272f7cc655
added auto scaling of dimension to voronoi_tessellation.f90. If dimension(i) == 0, it will get the smallest step of the remaining dimensions in that direction. Same is now done by patchFromReconstructedBoundaries
...
setup_processing.py now gives executables the extension *.exe and removes the old executables before compiling
2011-08-25 12:16:06 +00:00
39b97f7b67
added new script "spectral_scaleGeom"
2011-08-19 13:32:27 +00:00
2dd8a353bb
now also works for single stepped loadcases, too
2011-08-19 10:03:42 +00:00
483234047c
renamed to be consistent with other file names
2011-08-19 07:53:35 +00:00
878ea443d6
added options to control written output (mesh(default), box, points) and to control reconstruction. --undeformed does not calculate average deformation of ve, --scaling controls the scaling of the fluctuation. Default is 1.0, with 0.0 no fluctuations are shown
2011-08-19 07:33:22 +00:00
589666d19e
added python script to simply (without antialising) scale a given *.geom file
2011-08-19 07:29:50 +00:00
def056f30f
added $Id$ escaping
2011-08-18 08:00:19 +00:00
10493ea8b0
added $Id$ login name escaping
2011-08-18 07:53:07 +00:00
164da22730
added some warnings if Mentat or PIL cabability missing
2011-08-18 07:52:07 +00:00
ff3bac2496
does not double size of 2D data for spectral method any more
2011-08-17 16:06:50 +00:00
3298e6821b
added forgotten "averageDown" script to list of shortcuts
2011-08-16 18:20:52 +00:00
bedce496a8
deleted unused option "dimension" since this is figured out by parsing the respective spectral or geom input files anyway.
2011-08-15 08:34:55 +00:00
e851c934ee
interrogate output from DAMASK_spectral for "error divergence:" evolution.
...
results in an ASCII datafile containing the increment counter with corresponding divergence value. the evolution within an increment i is stretched between i and i+1 (and likely jumps then for the next increment i+1 at the value i+1)
2011-08-14 10:44:15 +00:00
a58158fc7c
bug fix related to hash key error when requesting --increments range
2011-07-25 12:53:54 +00:00
8cb8637d41
fixed running with default --range & using --separation options.
2011-07-22 11:59:03 +00:00
105a09f3df
only consider increments actually present in (spectral) result
...
added switch to change from range of positions to range of increments
polished help output
2011-07-21 15:45:41 +00:00
cb96559d99
mappings should rather be called "sumabs" and "avgabs" than "abssum" and "absavg"
2011-07-15 11:27:47 +00:00
63451f6edf
added "absavg" (=average of the absolute values) to the predefined mappings
2011-07-15 08:07:33 +00:00
b252b467e8
moved useful ASCII post processing scripts from my Code repository into the DAMASK post-proc.
2011-06-21 16:25:48 +00:00
daa45306f2
fixed a nasty bug in reading the header of spectral files
...
file type is guessed from extension
script tries to gracefully continue when data is missing but not crucial
ASCII file splitting now has correct increment (not position) in its name: poop_inc123.txt refers to true increment 123 not to position 123 in the result file anymore...
2011-06-21 12:38:58 +00:00
b6b02f6cf9
now reports over all loadcases present in the output file, not just a selected one...
2011-06-15 18:14:05 +00:00
dd11ef3072
added capability to deal with (Marc, spectral) output files that have less than every increment stored.
2011-06-15 17:49:59 +00:00
4d2f6782e7
removed "error" msg about missing /msc directory. if not present, now silently ignored (together with any requests for Mentat interaction).
2011-06-09 12:52:54 +00:00
d52fb1666e
"step" header is modified to "inc" to be consistent with postResults outputs.
2011-06-09 12:39:08 +00:00
a72a97f0ba
added nodalScalar capability.
...
--ns 'elements' reports the nodal connectivity (list of elements touching a node)
--prefix allows to prefix the output filename
2011-06-08 16:54:46 +00:00
f2c41aac83
do not use pwd library as it is unix only
...
due to the above path to subroutine is no longer set in mentat
spectral output uses different "grain numbers" for x-margin, y-margin, and margin edge
2011-06-08 15:06:37 +00:00
5740e13485
added new predefined mapping "abssum", which returns the sum of all absolute values
2011-06-08 14:07:15 +00:00
2cbfb69235
user defined mapping didn't work because of misspelled variable in "mapIncremental"
2011-06-08 13:56:21 +00:00
7aed307bfc
removed obsolet option --frame
...
renamed texture component frame to margin
2011-06-08 06:47:23 +00:00
8921f29dce
merged imaging, Mentat, and spectral versions of generating a patch of grains from a reconstructed boundary file into one script.
2011-06-07 19:15:34 +00:00
9370d9a049
added 'spectral_geomCheck' to symlink table
2011-06-06 15:29:14 +00:00
9f5c9142e5
generate ParaView mesh from spectral geometry (*.geom) description file to check for correct grain assignment
2011-06-01 08:44:14 +00:00
4a2078e965
produces two-column data of step | iteration at which convergence was reached (or aborted). Source file is the redirected STDOUT of DAMASK_spectral
2011-05-31 12:31:35 +00:00
ac3a2ee4e1
alll grain have the same phase now
...
corrected small bug when Image library is not available
2011-05-28 09:48:39 +00:00
5e6de13962
removed useless grain selection option (--grain)
...
escaped MPIE\n.... in $ID$
fixed problem when writing output files into subdirectories
2011-05-26 09:44:10 +00:00
96a4d83b9b
fixed bug in user mapping function.
...
corrected examples and added incremental mapping explanation to usage hint.
2011-05-26 09:41:53 +00:00
3ae8b0f732
dropped grain index from output header (did not work anyhow)
...
introduced string escape for the likes of MPIE\n.jia to prevent newline literal
2011-05-24 17:23:22 +00:00
73b1dbc5f7
corrected printing bug
2011-05-23 15:11:01 +00:00
6daa8b710b
shortened output header to contain only one description line that includes the command line switches passed to postResults
2011-05-23 07:13:28 +00:00
2a8dc1bff1
Fixed a bug introduced in the last revision. Datatype keywords (SCALARS, VECTORS, ...) were not written in plural form into VTK file, thus had not been recognized by ParaView.
2011-05-22 14:01:18 +00:00
8879e03a60
better reporting and fixed a bug when requesting non-existing data
2011-05-21 21:43:40 +00:00
25e3c4baf8
nicer STDOUTput and can now properly deal with resolutions of 1
2011-05-13 16:53:27 +00:00
4658190951
now reports command line switches used in generating output files
2011-05-12 09:53:44 +00:00
25bf6bf7dd
reports command line switches used in generating output files
2011-05-12 08:25:34 +00:00
0a5e151699
extension to make content clear
2011-05-11 17:10:12 +00:00
1928fa816c
output now contains both node and ip number
...
relation between node and ip numbering is given in new function ipIDs
2011-05-06 10:00:27 +00:00
56b2b3e572
* quick fix in order to make it work for elements of type 57 (quadratic with reduced integration)
...
* length for results of type "elemental tensor" is exactly 7
2011-05-05 09:16:29 +00:00
2fa4a344b9
added sloppy option to circumvent increment range check. that is useful to reach the last inc in former FFT sim results, which have a wrong inc count in their header
2011-05-02 16:10:18 +00:00
97185fb30b
removed debug printing...
2011-04-14 14:03:06 +00:00
f598e3b6ed
corrected mixup of Fortran (1,2,3,...) indexing to Python (0,1,2,...) style.
...
fixed erroneous array indexing in neighborhood calculation.
2011-04-14 10:11:23 +00:00
2a04bfdbba
fixed error concerning mixed-up index starts (0 vs 1)
...
"illegal" mapping functions are not checked anymore, i.e., if you want to average Euler angles you get it..!
sorting is now "natural", i.e., separating or sorting for "x,y,z" returns x varying fast and z slow. "--sort" overrides "--separate" in terms of output sorting (as expected...)
2011-04-13 16:31:44 +00:00
482be626e0
major recoding.
...
now with tiny memory footprint and better guessing of remaining time.
2011-04-12 17:46:35 +00:00
e4cbdc5743
now also works for 2D quadratic triangular element type 125
2011-04-06 13:11:25 +00:00
a989aa979e
corrected increment offsetting for spectral file type. should now be consistent with new (0+N) mpie_spectral output scheme
2011-04-06 10:00:19 +00:00
fffff202b0
2D mesh option now also alters corresponding job option
2011-04-04 15:47:09 +00:00
73accf79b6
added option for 2D mesh (3D extrusion step is skipped)
2011-04-04 12:12:01 +00:00
a546f8dbd2
changed wording of input queries to (hopefully) clearer statements
2011-03-30 15:48:18 +00:00
8467c91dab
moved setup scripts back to subdirs where they are needed
2011-03-22 15:42:53 +00:00
6ac2b4cf88
deleted useless storage of data resulting in partial speed-up.
...
large datasets (>1 mio elements) cannot be processed at the moment---investigation required, why each element eats about 20kB of memory?
2011-03-11 16:32:01 +00:00
1f0569230b
new installation/setup dir to hold all top-level make_*-scripts
2011-03-11 09:20:37 +00:00
621c1198a4
Marc/Mentat related things are only executed for filetype marc. this has been too lax before, causing trouble.
...
filetype is now lowercased right from the start, keeping effort low downstream...
2011-03-10 09:45:57 +00:00
6b63d488cb
need py_post, not py_mentat
2011-02-28 10:00:08 +00:00
9312aed545
set executable bit
2011-02-27 21:51:56 +00:00
8cae4d609a
some polishing for single precision version.
...
Now only prec and mpie_spectral are needed in different versions
2011-02-25 12:41:46 +00:00
cd5407b08b
removed all math functions only for double precision by the more flexible counterpart, e.g. "dsqrt --> sqrt", "dsin --> sin". Should not cause any harm, as long as "implicit none" is used.
...
Now it is possible to compile a single precision spectral solver/crystal plasticity by replacing mesh.f90 and prec.f90 with mesh_single.f90 and prec_single.f90.
For the spectral method, just call "make precision=single" instead of "make". Use "make clean" evertime you switch precision
2011-02-25 09:25:53 +00:00
f525c02ded
added single precision libraries for FFTW
...
First try of implement single precision crystal plasticity, not working yet.
polishing text about geometry construction.
polishing postResults, still having problems concerning machines without MSC installation
2011-02-24 18:39:57 +00:00
3c944e79fa
did some things Philip told me to do (postResults for spectral files is now working without MSC installation)
2011-02-22 16:18:47 +00:00
003b0def46
fixed a stupido error
2011-02-22 16:06:24 +00:00
333445be8b
changed meshing routine from python to fortran (f2py)
2011-02-22 16:00:22 +00:00
2dc3b35725
no need anymore for msc python interface when dealing with spectralOut files
2011-02-22 15:57:27 +00:00
8a5e28d5a6
now using fft reconstruction in 3Dvisualize, linear python code for reconstruction is removed
2011-02-21 21:03:21 +00:00
c399a06c97
tried to implement an inverse reconstruction, meaning to calculate an average deformation gradient from the 8 cornes of a node in reference and current configuration. Not working yet.
...
polishing, mainly in fft reconstruction.
added fftw library file and fortran file for type specification
added folder references with literature for visualization/geometry reconstruction
2011-02-21 16:30:18 +00:00
34336bc659
fixed bug when no argument was supplied
2011-02-18 09:02:29 +00:00
438bf95105
added ft-based reconstruction of deformed configuration to postprocessingMath.f90 and postprocessingMath.pyf
...
also added function to calculate average of tensor
some polishing in mpie_spectral.f90, added sanity check to see im resolution is power of 2
2011-02-14 17:21:31 +00:00
874d5067c7
switches drawing of servo links temporarily off to speed up processing.
2011-02-11 11:29:58 +00:00
3d7fad6ba9
mpie_spectral.f90: changed fourier transormation, now using the advanced interface to do the tranformation both ways with just one call. now also using the same variable for pk-stress in both domains and the change of deformation gradient in both domains.
...
postprocessing: renamed name of python/f2py modul from "reconstruct" to "postprocessingMath", added some numerical operations to use for postprocessing.
2011-02-09 17:47:28 +00:00
65ae7de42b
fixed bug for missing file name
2011-02-09 09:18:21 +00:00
7a7ca1aab7
mpie_spectral and mpie_interface: changed workingdir from pathToLoadFile to pathToGeomFile
...
mpie_spectral and numerics: added switch to prevent pre calculation of gamma_hat. slower, but saves memory
3Dvisualize: started to add support for gmsh (not fully working yet)
reconstruct: new version of f2py/Fortran subroutines for output of results from spectral method
2011-02-07 14:35:42 +00:00
6511b4b5a2
added possibility to select predefined colorscheme from config file
2011-02-04 09:36:26 +00:00
05350d573e
fixed a bug in new --sort
...
--separation option now keeps data point location as intact as meaningful after applying averaging
2011-02-01 18:32:20 +00:00
669893b397
binary \n will not longer cause problems
2011-02-01 18:25:40 +00:00
14d3d0ba99
less output but can now handle 'n/a' in data files without crashing
2011-02-01 16:43:00 +00:00
ebf201e89c
sorry---wrong file version uploaded before.
...
now correct: added "--sort" option extending the existing sorting according due to "--separation". This is, for instance, useful for separating /elementhomogeneous/ results by element and later listing the results along x,y,z...
slight improvement in file-error reporting
2011-02-01 13:24:19 +00:00
958c51fc26
added "--sort" option extending the existing sorting according due to "--separation". This is, for instance, useful for separating /elementhomogeneous/ results by element and later listing the results along x,y,z...
2011-02-01 13:13:05 +00:00
d6edb64929
produce VTK visualizations from postResults files
2011-02-01 10:48:44 +00:00
fec2c14a4e
removed hard-coded tolerances from mpie_spectral and put them to numerics/IO
...
removed storage of full cauchy stres field from mpie_spectral.f90, only average is stored now
added cauchy stress and von mises equivalent calculation to spectral post.
2011-01-31 17:07:42 +00:00
0cf1404db0
changed from element 117 (hexahedral with reduced integration) to element 7 (hexahedral with full integration). one-IP elements give trouble in the ping-pong scheme by being frequently/always outdated. this is likely connected to their finite element technology such that even high stiffness values do not prevent the solver from altering displacements in subsequent cycles. remedy is to go for full integration combined with /elementhomogeneous/ phase...
2011-01-27 14:32:29 +00:00
f84ac44f9a
now compiles Fortran sources on the fly
2011-01-26 15:17:03 +00:00
17b8205e3f
reconstruction of geometry is now working. It is directly implemented in python (file spectral_post.py).
...
reconstruction in fortran is not working (file reconstruct.f90) due to some problems with f2py
2011-01-26 12:56:52 +00:00
38409ebc71
included both new voronoi_xxx scripts
2011-01-24 16:21:49 +00:00
ee7a8ad52a
general speed up
...
spectral format records physical coordinates
tessellation now based on phys coords instead of discretization.
2011-01-24 16:21:17 +00:00
e721f2b654
sight restructuring of algorithm -- same functionality
2011-01-20 19:22:22 +00:00
52f06009be
rewrote voronoi tessellation. the VE generation is now divided into two parts:
...
-generation of random positions (in a 1x1x1 parameter space) and random angels
-actual voronoi tessellation
for the generation of the positions, a resolution can be specified to ensure that the tessellation works at that resolution. The tessellation will always work on larger resolutions but not necessarily at smaller ones.
2011-01-20 17:21:38 +00:00
5990d94ca3
added description and $ID$
2011-01-19 15:14:05 +00:00
581955f68d
parsePostFile screwed up the position of outputs for a grain number larger than 1
2011-01-14 13:15:19 +00:00
a662a28a97
added support for spectral method result files
2011-01-12 16:55:56 +00:00
322cbd2597
changed tools for voronoi tessellation to match new specification (file extension and resolution). Also renamed the files, because the space characters cause trouble
...
improved reconstruct.f90 and spectral_post.py, both files work now under linux
added make_reconstruct.py, small shell script for using f2py
2011-01-07 14:37:05 +00:00
20c00288b8
---had some problems with svn, hope everything is ok now---
...
added new tools to generate colormaps for paraview and gmsh, written in python
removed old fortran colormap generator.
removed test.py (not longer needed) and the python module reconstruct.pyd (not running under linux)
2011-01-05 14:53:31 +00:00
3b0eeb9956
added files for geometry reconstruction from given deformation gradient.
...
reconstruct.f90 is a fortran source file with comments to use with f2py
reconstruct.pyd is the compiled python module
test.py is a test file to check if the reconstruction works (unfortunately i does not at the moment, probably because there are type conflicts between fortran and python)
2010-12-22 10:59:54 +00:00
8d6ddcf28b
restructured (or structured at all) reading of header, first try of reading in deformation gradient
2010-12-09 15:53:52 +00:00
ac676623ca
reading of header is now working
2010-12-03 14:57:28 +00:00
2f47d4b09b
added conversion of 'resolution' in result file to integer array
2010-12-02 14:53:17 +00:00
ab82dbc99c
continued working on reading header of output file
2010-11-30 14:26:30 +00:00
9d32bc7ae6
added output of complete mesh to voronoi fast.f90
...
first version of post processing for spectral method: spectral_post.py
2010-11-29 10:08:53 +00:00
b2cc34349f
corrected mistake in nodeCoordinate assignment, bug-fix for option " --filter 'x,y,z' "
2010-11-12 08:56:04 +00:00
be265aef37
slight improvement of remaining time estimate...
2010-11-02 16:55:56 +00:00
796bffee2e
1) speed-up of elementTensor output
...
2) added estimate of remaining time
2010-11-02 15:45:23 +00:00
e80e055c75
added three small programs useful for spectral method (fortran sourececode)
...
voronoi fast.f90 and voronoi small memory.f90 are two variants to do a voronoi tessellation and write the result to a mesh file that can be interpreted by mpie_spectral.f90. Difference is the memory management resulting in one faster and one memory saving version. voronoi fast.f90 has also the ability to write out a file for Ricardo Lebensohns spectral code.
colormap.f90 is a simple code that can be used to generate colormaps for gmsh.
2010-10-31 16:06:10 +00:00
74743c0856
little typo fixed...
2010-10-26 18:26:04 +00:00
487466dd88
scripts are now aware of symbolic links to themselves...
2010-10-26 16:26:55 +00:00
18344b4e65
forgot linebreaks and executable flag...
2010-10-26 15:32:03 +00:00
5f64cacc74
generate alias in ~/bin for processing scripts
2010-10-26 15:02:20 +00:00
5be6604fd6
included definitions of "material" and "geometry" properties after mesh generation...
2010-10-26 13:45:23 +00:00
d965f14f90
generate hexahedral mesh in Mentat based on a data file of Ricardo Lebensohn's "txfft" format (phi1, Phi, phi2, x, y, z, grainId, phaseId)
...
Only grainId is used, physical dimension and subdivisions need to be provided on command line.
improvement possible by parsing x, y, z columns to get those directly.
2010-10-25 17:57:57 +00:00
3a4743f1d3
fixed bug in grain naming scheme. occurred for gaps in consecutive numberings...
2010-09-28 16:56:49 +00:00
23746e7e78
chopped line beginnings
...
introduced check for PIL, if this fails, no image will be output
2010-09-23 07:59:52 +00:00
d8a4cd1739
added optional argument --baseindex to specify what colorset should be influenced. identify colors, for instance, start at idx 214...
2010-09-13 07:31:27 +00:00
006dce7d8a
printing euler angles (and defgrads and orientations) for individual ips now works
2010-09-08 11:55:57 +00:00
4a31b175da
1) introduced default colors
...
2) complaints are more verbose
2010-09-06 16:51:11 +00:00
3b671faae7
added some flexibility (and bug-fixing) to the pre/post processing. the file 'MSCpath' contains the LOCAL path to the MSC installation (/msc in our case)
2010-08-17 18:21:22 +00:00
94efd57663
started to put Marc/Mentat related scripts to corresponding post/pre processing folders. We should move development from my Code folder to here (sorry, lost history then...)
2010-08-16 20:47:27 +00:00
Martin Diehl