Martin Diehl
669893b397
binary \n will not longer cause problems
2011-02-01 18:25:40 +00:00
Philip Eisenlohr
9d2b1aa942
specified line endings and activated $Id$
2011-02-01 16:43:58 +00:00
Philip Eisenlohr
14d3d0ba99
less output but can now handle 'n/a' in data files without crashing
2011-02-01 16:43:00 +00:00
Philip Eisenlohr
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
Philip Eisenlohr
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
Philip Eisenlohr
d6edb64929
produce VTK visualizations from postResults files
2011-02-01 10:48:44 +00:00
Martin Diehl
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
Christoph Kords
b72d75ed05
corrected indices in flux density output
2011-01-28 12:49:58 +00:00
Philip Eisenlohr
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
Philip Eisenlohr
c00a42a9a4
fixed last remaining tensor outputs to their transposed versions
2011-01-27 14:29:06 +00:00
Philip Eisenlohr
f84ac44f9a
now compiles Fortran sources on the fly
2011-01-26 15:17:03 +00:00
Martin Diehl
aa7d5afd94
changed ifort options back to previous ones as -fast is not working
2011-01-26 13:33:44 +00:00
Martin Diehl
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
Christoph Kords
05d4d5fef2
* kinetics now according to Mohles; enables solid solution hardening
...
* changed example in material.config accordingly
2011-01-26 10:17:42 +00:00
Philip Eisenlohr
6ac06af517
sorry, forgot to alter the most important part: post_results..! eventually, the t16 has defgrads and friends in correct format: du_i/dx_j is listed as linear array (1,1),(1,2),(1,3),(2,1),...
2011-01-24 16:23:37 +00:00
Philip Eisenlohr
38409ebc71
included both new voronoi_xxx scripts
2011-01-24 16:21:49 +00:00
Philip Eisenlohr
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
Philip Eisenlohr
4ee40df5ba
transposed writing into t16 file (and screen/debug) did not work with just stating array(1:3,:)...
...
fixed by using math_transpose3x3 (and similar) instead.
2011-01-24 15:02:19 +00:00
Philip Eisenlohr
794aeb253b
crystallite and CPFEM now print/output tensors in ij notation (used to be, particularly in the t16, ji, i.e. transposed...)
...
math reports random seed used for rnd() generation
2011-01-20 19:25:45 +00:00
Philip Eisenlohr
e721f2b654
sight restructuring of algorithm -- same functionality
2011-01-20 19:22:22 +00:00
Martin Diehl
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
Philip Eisenlohr
5990d94ca3
added description and $ID$
2011-01-19 15:14:05 +00:00
Christoph Kords
581955f68d
parsePostFile screwed up the position of outputs for a grain number larger than 1
2011-01-14 13:15:19 +00:00
Martin Diehl
b546d3bb91
some corrections regarding Multicore-support on mpie_spectral.f90 and makefile
...
slight changes on mpie_spectral.f90 regarding file output
2011-01-12 17:02:42 +00:00
Martin Diehl
a662a28a97
added support for spectral method result files
2011-01-12 16:55:56 +00:00
Christoph Kords
062c6b3d54
* allow zero cutoff radius
...
* added missing error codes
2011-01-12 12:36:48 +00:00
Christoph Kords
d92f992eb7
if nonlocal ip has local neighbor assume perfect transparency of the (perhaps virtual) phase boundary. this is only a quick fix! need to rethink the way constitutive and real boundaries are treated!
2011-01-11 18:00:29 +00:00
Christoph Kords
6a3dd78259
fluxes are disabled when keyword "/nonlocal/" is omitted in the material.config (this is a way to use the constitutive_nonlocal model as a local dislocation based model)
2011-01-11 14:55:36 +00:00
Martin Diehl
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
Martin Diehl
ccc6aac10b
did not delete mpie_spectral2.f90 properly
2011-01-07 13:11:32 +00:00
Martin Diehl
71fb5eedf2
changed makefile: material subroutines are now compiled using openmp, enabled linux threads for fftw3 (also did some slight changes to enable that in the fortran source code)
...
renamed mpie_spectral2.f90 to mpie_spectral2d.f90 (testing file, not properly working at the moment)
changed file extension and variable names in mpie_spectral.f90 and mpie_spectral_interface.f90 from "mesh" to "geom". Removed direct output from mpie_spectral.f90, all output is now base on materialpoint_results(:,1,:)
2011-01-07 12:56:45 +00:00
Franz Roters
9447330ab0
ATTENTION repository has been moved to new server!
...
Now at musvn1:/home/svn/repos/cpfem
2011-01-05 16:11:51 +00:00
Martin Diehl
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
Martin Diehl
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
Martin Diehl
8d6ddcf28b
restructured (or structured at all) reading of header, first try of reading in deformation gradient
2010-12-09 15:53:52 +00:00
Martin Diehl
ac676623ca
reading of header is now working
2010-12-03 14:57:28 +00:00
Martin Diehl
2f47d4b09b
added conversion of 'resolution' in result file to integer array
2010-12-02 14:53:17 +00:00
Franz Roters
096a4803b0
some polishing concerning openMP initialization
...
number of threads is now written to *.out file
2010-12-02 11:04:29 +00:00
Martin Diehl
ab82dbc99c
continued working on reading header of output file
2010-11-30 14:26:30 +00:00
Martin Diehl
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
Denny Tjahjanto
36f27edd75
1. Adding more comments
...
2. Changing 'subroutine' into 'function' for optimizing the computation time
2010-11-26 11:50:20 +00:00
Christoph Kords
6b199e7830
allow 2digit values to be read in for MPIE_NUM_THREADS
2010-11-25 10:34:55 +00:00
Franz Roters
4216a3e066
corrected name of include file
2010-11-24 07:55:08 +00:00
Franz Roters
c07a91d19e
removed -save option from the compile statement using OpenMP
2010-11-23 16:25:58 +00:00
Christoph Kords
ea6a6c9503
* environment variable MPIE_NUM_THREADS now determines number of threads used in parallel regions
...
* default value of the OMP_NUM_THREADS variable has to be restored at the end of mpie subroutine, since marc also seems to use and change(!) this
* usage: "export MPIE_NUM_THREADS=<number of threads>" to set variable in shell, then restart mentat and compile with option 3 (at the moment this does only work on ws 6, since all other workstations use compiler option "-save"; this puts all local variables by default in static memory, which is a killer for parallelization!)
2010-11-19 17:45:27 +00:00
Christoph Kords
b9b226548c
* enclose wall time measurement for LpAndItsTangent and count of leapfrogbreaks by CRITICAL construct
...
* better use SINGLE (having an implicit barrier at the end) instead of MASTER construct
* deleted all explicit BARRIERs after do loops since parallel loop construct implies barrier at the end
2010-11-19 17:29:29 +00:00
Christoph Kords
15728d3305
enclose wall time measurement for dotState and dotTemperature function by CRITICAL construct
2010-11-19 17:19:03 +00:00
Philip Eisenlohr
b2cc34349f
corrected mistake in nodeCoordinate assignment, bug-fix for option " --filter 'x,y,z' "
2010-11-12 08:56:04 +00:00
Christoph Kords
4f267f4375
parallelization now works properly (not yet tested extensively though):
...
* had to add some BARRIER constructs
* only the master thread is allowed to increase the state counter
yet parallelization seems not to give a significant decrease in calculation time with nonlocal model (because of too many CRITICAL statements?)
2010-11-11 16:16:05 +00:00
Christoph Kords
14594dc9e2
* OpenMP does not like unbounded array subscripts in function calls
...
* some write statements were still not enclosed by CRITICAL construct, some needed a MASTER construct
2010-11-11 13:14:53 +00:00