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
b72d75ed05
corrected indices in flux density output
2011-01-28 12:49:58 +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
c00a42a9a4
fixed last remaining tensor outputs to their transposed versions
2011-01-27 14:29:06 +00:00
f84ac44f9a
now compiles Fortran sources on the fly
2011-01-26 15:17:03 +00:00
aa7d5afd94
changed ifort options back to previous ones as -fast is not working
2011-01-26 13:33:44 +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
05d4d5fef2
* kinetics now according to Mohles; enables solid solution hardening
...
* changed example in material.config accordingly
2011-01-26 10:17:42 +00:00
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
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
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
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
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
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
a662a28a97
added support for spectral method result files
2011-01-12 16:55:56 +00:00
062c6b3d54
* allow zero cutoff radius
...
* added missing error codes
2011-01-12 12:36:48 +00:00
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
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
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
ccc6aac10b
did not delete mpie_spectral2.f90 properly
2011-01-07 13:11:32 +00:00
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
9447330ab0
ATTENTION repository has been moved to new server!
...
Now at musvn1:/home/svn/repos/cpfem
2011-01-05 16:11:51 +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
096a4803b0
some polishing concerning openMP initialization
...
number of threads is now written to *.out file
2010-12-02 11:04:29 +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
36f27edd75
1. Adding more comments
...
2. Changing 'subroutine' into 'function' for optimizing the computation time
2010-11-26 11:50:20 +00:00
6b199e7830
allow 2digit values to be read in for MPIE_NUM_THREADS
2010-11-25 10:34:55 +00:00
4216a3e066
corrected name of include file
2010-11-24 07:55:08 +00:00
c07a91d19e
removed -save option from the compile statement using OpenMP
2010-11-23 16:25:58 +00:00
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
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
15728d3305
enclose wall time measurement for dotState and dotTemperature function by CRITICAL construct
2010-11-19 17:19:03 +00:00
b2cc34349f
corrected mistake in nodeCoordinate assignment, bug-fix for option " --filter 'x,y,z' "
2010-11-12 08:56:04 +00:00
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
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
f3a7e323ea
pretty printing of hex slip/twin systems
2010-11-09 13:20:33 +00:00
bca951325d
* cpfem will now produce a warning if all components of the stiffness tensor are close to zero
2010-11-04 18:18:01 +00:00
54748bb6a5
cycleCounter reset was 0 now -1 ==> occurence of jacobian calculation is now correct, i.e. always calc at start of new increment
2010-11-04 18:15:50 +00:00
e561358a90
presentation to illustrate Marc/Mentat interface for restarting capability (thx to Denny!)
2010-11-04 15:31:29 +00:00
Philip Eisenlohr