Martin Diehl
d451a3a7a0
calculation of shearrates etc in one function
2018-11-27 18:25:06 +01:00
Martin Diehl
bfad81848a
kinetics similar to phenopowerlaw
2018-11-26 22:36:32 +01:00
Yi-Chin Yang
c13b9718f6
Merge remote-tracking branch 'origin/development' into NOSlipFCC
2018-11-26 16:10:14 +01:00
Vitesh Shah
ef1e9cce0d
Functionality to avoid creating datasets in HDF5 of zero dimensions
2018-11-26 15:38:31 +01:00
Martin Diehl
cf4a0a69fd
Intel compiler detected use of unallocated variable
...
bug was in for a while
2018-11-26 09:15:44 +01:00
Martin Diehl
c63d297145
copied from phenopowerlaw
2018-11-26 07:10:43 +01:00
Martin Diehl
8f59a40f48
using new parameter structure for output
2018-11-26 01:34:38 +01:00
Martin Diehl
4c46f3daa5
adopting argument parsing to Phenopowerlaw
2018-11-26 00:14:27 +01:00
Martin Diehl
9f81fa8e9a
prevent segmentation fault
...
post_results better readable
2018-11-25 21:39:59 +01:00
Martin Diehl
6f93f8de04
starting to introduce parallel structure for new style parameter reading
2018-11-25 19:37:46 +01:00
Martin Diehl
2b9a34eac4
Merge remote-tracking branch 'origin/development' into 20-NewStyleDislotwin
2018-11-25 17:12:26 +01:00
Martin Diehl
7cc2892e64
no need to calculate twinned volume fraction as state
2018-11-25 16:14:46 +01:00
Martin Diehl
8a253856f1
no need to perform the state integration
2018-11-25 15:38:14 +01:00
Martin Diehl
8ccec40e35
Merge remote-tracking branch 'origin/development' into 19-NewStylePhenopowerlaw
2018-11-25 13:39:34 +01:00
Martin Diehl
1808b37357
can be easily computed during postprocessing
2018-11-25 11:17:12 +01:00
Martin Diehl
917453d191
polishing
2018-11-25 11:14:09 +01:00
Martin Diehl
c0ce95183c
can be easily calculated during post processing
...
and does not have to be a state
2018-11-25 10:59:36 +01:00
Martin Diehl
1d936fc5f2
simple variable rename causes problem with GCC 7.3
...
probalby a compiler bug
2018-11-24 09:53:18 +01:00
Yi-Chin Yang
efe90f86e2
missing comma
2018-11-23 19:00:37 +01:00
Yi-Chin Yang
f0359d36b6
adding FCC {110}<110> family
2018-11-23 18:53:50 +01:00
Vitesh Shah
e7581f06d9
Fixed error with parallel write
2018-11-23 17:54:02 +01:00
Vitesh Shah
0e55bd6140
started read parallelization but getting errors
2018-11-23 15:49:43 +01:00
Franz Roters
9760b138c2
added link lost during merge
2018-11-23 15:02:35 +01:00
Franz Roters
a498e5580a
Merge remote-tracking branch 'remotes/origin/54-populate-grains-deprecated' into development
2018-11-23 14:38:22 +01:00
Martin Diehl
c7fb868b40
state layout was broken
2018-11-23 07:01:04 +01:00
Martin Diehl
6df68d9428
no need do constantly convert 3x3 matrix <-> 6 vector
2018-11-23 06:04:39 +01:00
Martin Diehl
8b4781cf28
no need to repeat code
2018-11-23 05:37:31 +01:00
Vitesh Shah
a6bc2a16f1
Parallel works for groups and included write function for integer data
2018-11-21 19:35:37 +01:00
Vitesh Shah
5cc6d86c61
Added the parallelized functionality for integer datatypes (works for groups?)
2018-11-21 16:57:36 +01:00
Vitesh Shah
d00e3105ed
Made hdferr < 0 and removed the unnecessary call to h5getspace in hyperslab
2018-11-21 16:10:17 +01:00
Martin Diehl
e2e045f75b
using 'service function' for allocation of state
2018-11-20 19:46:55 +01:00
Martin Diehl
d0be1d61db
Merge branch 'development' into 20-NewStyleDislotwin
2018-11-20 19:41:30 +01:00
Martin Diehl
ecde7e5953
Merge remote-tracking branch 'origin/development' into 20-NewStyleDislotwin
2018-11-20 19:38:18 +01:00
Martin Diehl
baf56a0e94
respect dependencies of inclusion
2018-11-20 15:27:32 +01:00
Martin Diehl
50a1ba6287
hdferr needs to be integer type during compile time of HDF5 library
2018-11-20 15:26:49 +01:00
Vitesh Shah
c0481307ee
Changed the intent of input argument of dataset to inout
2018-11-20 11:54:51 +01:00
Martin Diehl
be8e6cbe57
Merge branch 'development' into 19-NewStylePhenopowerlaw
2018-11-19 20:27:10 +01:00
Martin Diehl
ba5b623f14
Merge branch 'development' into improved-HDF5
2018-11-18 19:44:36 +01:00
Martin Diehl
dd6f6bba9d
planning the forwarding of the results
2018-11-18 12:41:05 +01:00
Martin Diehl
dfd624e520
create results file
2018-11-18 11:58:49 +01:00
Martin Diehl
fb5b1bfd8e
compile results module
2018-11-18 11:32:53 +01:00
Martin Diehl
35211a8468
compilation order should reflect dependency
2018-11-18 11:24:55 +01:00
Martin Diehl
366c63e273
CONFIG is a keyword for Cmake
2018-11-18 11:21:43 +01:00
Martin Diehl
8b1785c05c
parallel writing for pReal
2018-11-18 09:47:50 +01:00
Martin Diehl
73ca289322
F_aim was wrong in case of MPI parallelization
2018-11-17 16:50:19 +01:00
Martin Diehl
d8a425b464
prevent error
...
h5close_f will give an error if h5open_f was not called.
according to the manual, calling h5open_f multiple time is not a problem
2018-11-17 15:00:51 +01:00
Martin Diehl
89504c3c52
will not work with python3
2018-11-17 08:08:27 +01:00
Martin Diehl
d5963b403a
setting the property when opening the file
2018-11-16 07:05:43 +01:00
Martin Diehl
9b32fe6dbd
MPI file access needed for output
2018-11-15 07:28:34 +01:00
Martin Diehl
b48bd3a082
keep only general functionality in HDF5_results
2018-11-15 07:07:59 +01:00
Martin Diehl
58f2a25ffd
clearer separation of tasks
...
1) general HDF5_utilities (for results, restart, ...)
2) results related helpers (based on HDF5_utilities)
2018-11-15 06:49:31 +01:00
Martin Diehl
e6c7a73886
first step to include HDF5 in MSC.Marc
2018-11-14 13:12:08 +01:00
Martin Diehl
fc5300817d
Merge remote-tracking branch 'origin/development' into 42-new-coding-style-for-homogenization-NEW
2018-11-13 23:42:39 +01:00
Franz Roters
174997563e
Merge branch '47-save-restart-information-in-hdf5' into 'development'
...
Resolve "save restart information in HDF5"
See merge request damask/DAMASK!44
2018-11-12 08:42:22 +01:00
Franz Roters
475ca29f3c
close all objects of type integer(HID_T) created by native hdf5 routines
2018-11-09 16:11:19 +01:00
Franz Roters
71cb635b10
missing link
2018-11-07 15:13:16 +01:00
Franz Roters
e18f39d64c
added support for Marc 2018.1
2018-11-07 14:11:10 +01:00
Martin Diehl
7386c6ff9c
no need for "use" statements
...
will be inherited from outer function
2018-11-04 20:16:12 +01:00
Martin Diehl
811a02eb51
better readable
2018-11-04 09:19:24 +01:00
Martin Diehl
c863419a0b
cleaning
2018-11-04 09:00:35 +01:00
Martin Diehl
7a37ea25f3
cleaning
2018-11-04 08:49:40 +01:00
Martin Diehl
78f4d4c5ee
simplified
2018-11-04 08:22:59 +01:00
Martin Diehl
ff2b693f1e
simplified
2018-11-04 07:56:27 +01:00
Martin Diehl
572576cf88
these functions have no side effects
2018-11-04 07:33:57 +01:00
Martin Diehl
4dc3761fea
setting of (random) cluster orientations remove
...
can be easily done in a pre processing step
2018-11-04 07:29:39 +01:00
Martin Diehl
f471911e01
fixed array out of bounds during initialization
...
happened for inactive homogenization
2018-11-04 07:27:25 +01:00
Martin Diehl
3c11905f63
cleaning
2018-11-03 23:43:25 +01:00
Martin Diehl
8127d85be1
simplifying
2018-11-03 23:13:20 +01:00
Martin Diehl
078729bfa3
using dependent state
2018-11-03 22:43:52 +01:00
Martin Diehl
70998f7f9a
cleaning
2018-11-03 22:00:40 +01:00
Martin Diehl
69079b6558
was not used at all
2018-11-03 21:16:36 +01:00
Martin Diehl
0aa21e507a
functions only used within updatestate
2018-11-03 21:11:43 +01:00
Martin Diehl
c16fdec749
misplaced debug statement
...
debugging the individual constituents should be done by the constitutive
laws
2018-11-03 21:00:54 +01:00
Martin Diehl
85a2f19b18
wrong averaging
2018-11-03 20:03:23 +01:00
Martin Diehl
0447ea9d74
starting to introduce state pointers
2018-11-03 20:00:02 +01:00
Martin Diehl
fdbc20b739
cleaned and unified
2018-11-03 19:29:00 +01:00
Martin Diehl
08c692bfc1
passing in instance simplifies things a lot
2018-11-03 16:40:17 +01:00
Martin Diehl
a0c3a65b39
Noutput not needed anymore
2018-11-03 15:41:54 +01:00
Martin Diehl
a678e9b94f
bugfixes
...
flaws introduced in 42-xx branch and detected by new test
- outputsize was wrong
- nconstituents is product along the 3 directions, not the sum
2018-11-03 15:13:11 +01:00
Martin Diehl
c5b98e544a
Merge remote-tracking branch 'origin/development' into 42-new-coding-style-for-homogenization
2018-11-03 11:30:12 +01:00
Martin Diehl
f88b78195a
not need to repeat shape definition
2018-10-26 10:20:45 +02:00
Martin Diehl
c413df7ed2
a little cleaning
2018-10-18 22:54:33 +02:00
Martin Diehl
49c7a6c524
sorting
2018-10-18 22:20:26 +02:00
Martin Diehl
f71e067f5e
using functionality from lattice
2018-10-18 21:41:17 +02:00
Martin Diehl
d7fa3a9791
just sorting
2018-10-18 21:34:26 +02:00
Martin Diehl
59399f2cd2
sanity checks
2018-10-18 19:06:21 +02:00
Martin Diehl
33d816a5d0
polishing
2018-10-17 23:56:57 +02:00
Martin Diehl
121d146a35
using more central functions
2018-10-17 23:20:39 +02:00
Martin Diehl
4320061554
using central functionality
2018-10-17 23:13:47 +02:00
Martin Diehl
8017fccbc3
adjusting to phenopowerlaw
2018-10-17 22:18:33 +02:00
Martin Diehl
2fd348f251
just re-ordering
2018-10-17 21:47:50 +02:00
Martin Diehl
6b4b5b511f
don't let dipole formation be undefined
2018-10-17 21:41:21 +02:00
Martin Diehl
b7bc5afa90
Merge remote-tracking branch 'origin/development' into 20-NewStyleDislotwin
2018-10-17 20:10:49 +02:00
Martin Diehl
7ac96bd630
polishing
2018-10-17 00:01:21 +02:00
Martin Diehl
d92ab04892
Merge branch 'development' into 19-NewStylePhenopowerlaw
2018-10-15 22:38:55 +02:00
Martin Diehl
be8d7e19fe
missing use statement caused compilation error
2018-10-15 22:38:34 +02:00
Martin Diehl
2ca7807438
cleavage systems avaialable as function
2018-10-15 22:29:23 +02:00
Martin Diehl
ad1a64c338
rename was missing
2018-10-15 06:01:38 +02:00
Martin Diehl
57ca93c3c9
added missing allocation
2018-10-15 05:28:52 +02:00
Martin Diehl
e5ec5ed54c
check prerequisites as early as possible
2018-10-15 05:03:53 +02:00
Martin Diehl
767ca0edd4
avoid code repetition
2018-10-14 20:16:30 +02:00
Martin Diehl
a311cc7673
Merge remote-tracking branch 'origin/development' into 19-NewStylePhenopowerlaw
2018-10-14 15:53:50 +02:00
Martin Diehl
6599aa487c
kinetics_slip does not seem to be the problem
2018-10-14 15:53:24 +02:00
Martin Diehl
1839228add
sorting
...
new mappings. Needs more work (don't forget https://xkcd.com/927 )
2018-10-14 12:27:42 +02:00
Martin Diehl
943fa764c0
not needed
2018-10-14 11:15:07 +02:00
Martin Diehl
4ac12d1dc5
removed hybridIA
2018-10-14 10:11:26 +02:00
Martin Diehl
4c780226d1
polishing
2018-10-14 09:27:34 +02:00
Martin Diehl
17c21dfc92
mesh_element should not be used anymore
2018-10-13 18:21:13 +02:00
Martin Diehl
06d71d9d2c
Merge remote-tracking branch 'origin/development' into 42-new-coding-style-for-homogenization
2018-10-13 17:27:25 +02:00
Martin Diehl
513faa2218
investigating the reason for the poor performance
2018-10-13 11:29:07 +02:00
Franz Roters
d9fbf2ab50
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-10-13 10:52:03 +02:00
Franz Roters
905b5a33eb
added link for Marc2018
2018-10-13 10:50:20 +02:00
Franz Roters
1a458108bc
added support for Marc2018
2018-10-13 10:26:42 +02:00
Martin Diehl
fa9cbe13c9
do concurrent not safe here
2018-10-13 08:35:36 +02:00
Martin Diehl
6207781eb6
[skip ci] [skip sc] streamlining
2018-10-12 17:24:46 +02:00
Franz Roters
d1561e867a
Merge branch '44-column-major-access-to-interaction-matrices' into 'development'
...
Resolve "column-major access to interaction matrices"
Closes #44
See merge request damask/DAMASK!43
2018-10-12 13:03:22 +02:00
Martin Diehl
5597ee338e
simplifying
...
- whole type is private, no need to make the components private
- target not needed when using associate
2018-10-12 08:26:24 +02:00
Martin Diehl
854d3ab774
clearer logic for controlling dipole formation
2018-10-12 08:23:56 +02:00
Vitesh Shah
38f9807c42
cleaning up the codes
2018-10-11 18:00:01 +02:00
Martin Diehl
51a3308a3d
bugfix: same tau was used for all systems
2018-10-11 17:28:00 +02:00
Vitesh Shah
cd2fb9f417
adding prefixes to writing of datasets under groups
2018-10-10 19:08:44 +02:00
Vitesh Shah
fa13e71493
removed typo(&)
2018-10-10 18:09:10 +02:00
Vitesh Shah
78db0ecb3c
reading and writing of subroutines of datasets of different ranks
2018-10-10 17:54:55 +02:00
Martin Diehl
7217cdac1b
Merge branch '29-rename-mesh_element-array' into development
2018-10-10 15:55:23 +02:00
Martin Diehl
ef1899cfd0
improved error reporting
...
- tell user specifically whether number of systems is too high or too
low
- c/a sanity check can be done centrally
2018-10-10 07:42:07 +02:00
Philip Eisenlohr
f0b9c0caf7
polishing and encapsulating of Abaqus-specific local variables
2018-10-09 17:57:10 -04:00
Franz Roters
1d0584dcdf
typo in message
2018-10-09 16:44:47 +02:00
Martin Diehl
d81c9f744e
avoid errors by consistently closing objects
...
general bugfixes, works now
2018-10-09 16:17:21 +02:00
Martin Diehl
2064ed80fd
more flexible file open routine
2018-10-09 15:15:08 +02:00
Martin Diehl
8626ba1ce6
adjusted variable names and added comments
2018-10-09 14:46:57 +02:00
Martin Diehl
48403c37b7
initializing and closing the HDF5 library
2018-10-09 14:13:51 +02:00
Martin Diehl
e7e0cf7118
preprocessor statements are case sensitive..
...
... and the PETSc variable is PETSc, not PETSC
2018-10-09 13:57:29 +02:00
Martin Diehl
e3943d282a
Merge branch 'development' into 47-save-restart-information-in-hdf5
2018-10-09 13:51:58 +02:00
Vitesh Shah
09ba95b029
non-executable
2018-10-09 10:58:45 +02:00
Vitesh Shah
70c746a8f1
Interfacing for subroutines to identify real and integer data
2018-10-09 10:57:06 +02:00
Martin Diehl
e71ec2057e
hotfix
...
does not seem to cause trouble here, but better be safe
2018-10-09 08:37:03 +02:00
Martin Diehl
4da3923440
more straight-forward definition of array shape
2018-10-08 23:54:45 +02:00
Martin Diehl
05aae4f7eb
removed forgotten debug statement
2018-10-08 23:34:05 +02:00
Martin Diehl
50f0362aa1
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-10-08 23:33:31 +02:00
Martin Diehl
6064daceda
do error check earlier
...
array out of bounds happened already before the check
2018-10-08 23:32:46 +02:00
Martin Diehl
b51c3b5a07
handling files without LF at last line correctly
2018-10-08 23:27:06 +02:00
Philip Eisenlohr
225a48021d
foramtting of random_seed output to align with rest
2018-10-08 17:01:56 -04:00
Martin Diehl
78b5fa466d
fixed hick-up in lattice definitions
2018-10-08 19:32:18 +02:00
Martin Diehl
71d80cc203
sanity check too strong
...
still makes sense to have, but some tests need to be updated before
setting it active again
2018-10-08 13:27:15 +02:00
Martin Diehl
2a79f7ea61
taking care of c/a for hex and bct
2018-10-08 13:25:40 +02:00
Martin Diehl
5efcad952a
avoid type conversion
...
more systematic checking of correct input parameters
2018-10-08 08:28:56 +02:00
Martin Diehl
302cf9b6c2
bugfix, now inline with code from initializeStructure
2018-10-07 22:16:18 +02:00
Martin Diehl
facb1c2407
unified error checking
2018-10-07 20:03:45 +02:00
Martin Diehl
0e17b17b90
using full precision
2018-10-07 18:40:02 +02:00
Martin Diehl
34e0aca564
using full precision
2018-10-07 18:18:24 +02:00
Martin Diehl
a53488d666
calculation of slip/twin/trans/damage-coordinate system was wrong
2018-10-07 15:44:13 +02:00
Martin Diehl
032c35a499
WIP: Calculation of elasticity matrices for twin and trans
2018-10-06 10:42:25 +02:00
Martin Diehl
3a2f86df1c
non-schmid is always bcc
2018-10-05 23:04:06 +02:00
Martin Diehl
81a77b508a
check for array size was off by 1
2018-10-05 22:46:53 +02:00
Martin Diehl
3a39c2d68c
[skip ci] strings will be converted to lowercase
2018-10-05 17:54:24 +02:00
Martin Diehl
0cddf36106
typo, using value of p_slip for q
2018-10-05 17:47:46 +02:00
Martin Diehl
2a51e528bd
Merge branch 'development' into 44-column-major-access-to-interaction-matrices
2018-10-05 11:11:20 +02:00
Martin Diehl
40cf7f134d
Merge branch '44-column-major-access-to-interaction-matrices' of magit1.mpie.de:damask/DAMASK into 44-column-major-access-to-interaction-matrices
2018-10-05 11:08:29 +02:00
Martin Diehl
aa34b83978
avoid array out of bounds
...
allows to handle un- or partly defined non schmid coefficients
2018-10-05 10:36:44 +02:00
Martin Diehl
a560fff2ac
[skip ci] [skip sc] generic interfaces for HDF5 operations
...
most existing HDF5 functions are only designed for writing output files
and more general functionality is needed for storing the restart data
2018-10-05 09:03:20 +02:00
Martin Diehl
47a748831f
Cannot write out resolved stress including non schmid
...
Would need to write out tau_pos and tau_neg separately, so limit to
resolved stress in the sense of Schmid
2018-10-05 07:26:13 +02:00
Martin Diehl
6ccbc5e4e6
characteristing shear is defined per system
2018-10-05 07:20:51 +02:00
Martin Diehl
e74a852ef7
new style is noticable slower
...
runtime (polycrystal example) increased from 70s to 90s.
Need to figure out how to write readable AND fast code
2018-10-05 06:38:20 +02:00
Martin Diehl
3a00fcde16
WIP: characteristic shear calculated in lattice
2018-10-05 05:09:31 +02:00
Martin Diehl
fda2676365
calculate only if required
2018-10-05 04:54:47 +02:00
Martin Diehl
79b0efc678
also using nonSchmid matrix from lattice
2018-10-04 21:20:42 +02:00
Martin Diehl
35956ce0c8
needs to be multiplikative
2018-10-04 21:06:57 +02:00
Vitesh Shah
a6ebdfc1ff
Reading and writing not very efficient method
2018-10-04 17:00:24 +02:00
Martin Diehl
edb94b8bc9
using the same kinetics for slip everywhere
2018-10-04 16:12:14 +02:00
Martin Diehl
12d053d123
strange parameter caused simulation not to start if not defined
...
probably a leftover. now only bothers you if twin systems are defined
2018-10-04 14:54:00 +02:00
Martin Diehl
f1750b4eef
avoiding use of unitialized arrays
2018-10-04 14:51:32 +02:00
Martin Diehl
cd7f285686
prevent use of unallocated array
...
interaction matrices are used even if no slip/twin systems defined
2018-10-04 07:17:00 +02:00
Martin Diehl
dc289a278b
clearer code
...
homogenization is defined per element, not per IP
hence, use material_homogenizationAt instead of deprecated
material_homog which pretends a dependency on the integration point
2018-10-04 06:48:33 +02:00
Martin Diehl
df473302f4
consistent (and better understandable) names
2018-10-04 06:03:48 +02:00
Martin Diehl
f494803c76
fixed tangent calculation
2018-10-03 16:40:40 +02:00
Martin Diehl
ec77d4450b
checking for required things
2018-10-03 09:17:06 +02:00
Martin Diehl
db645efda7
function to calculate nonSchmid matrix
2018-10-03 08:49:23 +02:00
Martin Diehl
103b46ecbb
don't repeate code over different plastic_xxx modules
2018-10-03 07:58:02 +02:00
Martin Diehl
89e03be699
also building interactions for slip-twin and twin-slip
2018-10-03 07:51:11 +02:00
Martin Diehl
0cac6ce0d2
allocating to wrong size
2018-10-03 07:36:57 +02:00
Martin Diehl
ba383336f5
only one projection matrix required for all (non)schmid projections
...
following the example from plastic_nonlocal
2018-10-03 07:23:18 +02:00
Martin Diehl
d718521890
polishing
2018-10-03 07:13:30 +02:00
Martin Diehl
67c55965e9
segmentation fault, tried LHS with scalar
...
Nxx per system is stored in LATTICE_yyy_NxxxSYSTEM
Using CAPITALS for parameters improves readbility
2018-10-03 06:45:42 +02:00
Martin Diehl
97d2c2b353
mesh_element(1) has a size of 4
...
fixed sigsegv
2018-10-02 22:54:44 +02:00
Martin Diehl
740325165d
Merge branch 'development' into 20-NewStyleDislotwin
2018-10-02 20:37:11 +02:00
Martin Diehl
0047e6c14c
correct calculation of dGdot_dTau
...
needs to take drag coefficient into account for the derivative
2018-10-02 18:09:57 +02:00
Martin Diehl
9cc3a77e01
mean velocity of dislocation corrected for higher temperatures
...
- v_wait: obstacle spacing over waiting time (effective velocity)
- v_run: 'true' velocity of a dislocation moving between obstacles
for higher temperatures v_run becomes important and bounds the velocity
2018-10-02 15:56:18 +02:00
Martin Diehl
5d1e648d17
drag coefficient B introduced and read in
2018-10-02 15:25:40 +02:00
Martin Diehl
384a785805
more sevice functions use
...
still trying
2018-10-01 23:02:31 +02:00
Martin Diehl
ff2614b757
using 'service functions' from lattice
...
get rid of many global array on the long run
2018-10-01 22:48:14 +02:00
Martin Diehl
a943940760
Merge branch 'development' into 44-column-major-access-to-interaction-matrices
2018-10-01 22:32:38 +02:00
Martin Diehl
98296a978e
Merge branch 'development' into 29-rename-mesh_element-array
2018-10-01 22:29:31 +02:00
Martin Diehl
e7b8adadd0
passing in instance and of is much easier
...
also using always the full Mandel stress tensor to avoid repeated
conversion
2018-10-01 22:14:54 +02:00
Martin Diehl
ddfbe19f02
Merge branch 'development' into 20-NewStyleDislotwin
2018-10-01 21:16:11 +02:00
Martin Diehl
27bde05529
Merge branch '38-introduce-rudimentary-PETSc-based-FEM-solver' into 'development'
...
Resolve "Introduce Rudimentary PETSc based FEM solver"
Closes #38
See merge request damask/DAMASK!41
2018-10-01 21:10:03 +02:00
Franz Roters
d5cf7d2144
Merge branch '19-NewStylePhenopowerlaw' into 'development'
...
Resolve "New coding style for phenopowerlaw"
See merge request damask/DAMASK!39
2018-10-01 16:08:05 +02:00
Martin Diehl
49a4202e26
not needed here
2018-09-29 10:31:45 +02:00
Martin Diehl
14509876d4
error == 0 is normal termination
2018-09-29 09:57:53 +02:00
Martin Diehl
637544509e
rewind needed
...
the file is read twice from beginning to end
2018-09-28 10:06:21 +02:00
Vitesh Shah
8619b49e20
made it nonexecutable
2018-09-28 08:47:00 +02:00
Vitesh Shah
8bd0cb3b25
Made it non executable
2018-09-28 08:45:52 +02:00
Vitesh Shah
6efde5f4ce
Combined the dataset creation and writing dataset in single subroutine. Created group for looped writing.
2018-09-28 08:41:29 +02:00
Martin Diehl
b141252f1f
doing the same things in the same way in both solvers
2018-09-28 08:06:22 +02:00
Martin Diehl
d1f614991e
merging good style mutually FEM <-> Spectral
2018-09-28 07:49:52 +02:00
Martin Diehl
6780217193
adopted from DAMASK_spectral
2018-09-28 07:25:32 +02:00
Martin Diehl
e8f687a99c
typo, Abaqus/Marc did not compile
2018-09-27 20:18:37 +02:00
Martin Diehl
822bc4dbbf
Merge branch 'development' into 29-rename-mesh_element-array
2018-09-27 20:12:09 +02:00
Martin Diehl
4c057ba529
both solvers can share quit
2018-09-27 20:09:59 +02:00
Chuanlai Liu
b391ad304d
tensor input/output is of the form 11,12,13, 21,22,23, 31,32,33. Thus row index i is slow, while column index j is fast
2018-09-27 11:53:45 +02:00
Martin Diehl
6509775d41
cleaning
...
avoiding compiler warnings and confusion of the user
2018-09-27 08:23:30 +02:00
Vitesh Shah
5b52ac91e5
Revert "made it non executable"
...
This reverts commit 05c784aefc
.
2018-09-26 09:28:47 +02:00
Martin Diehl
c3380aa9f3
using wrong (uninitialized) stress
2018-09-26 09:22:12 +02:00
Vitesh Shah
2edf8f7f53
Need to commit this before reverting the changes in other files
2018-09-26 09:21:53 +02:00
Vitesh Shah
3408282755
windows and linux editor incompatibility
2018-09-25 16:49:25 +02:00
Vitesh Shah
05c784aefc
made it non executable
2018-09-25 16:48:02 +02:00
Vitesh Shah
b71896d4b6
writing of multidimensional dataset into hdf5 does not appear as multiple pages in hdfviewer
2018-09-25 16:42:43 +02:00
Martin Diehl
ceb385ef39
calcMode not needed for spectral and FEM
2018-09-23 21:31:30 +02:00
Martin Diehl
c42eb87a33
using arrays with new names
2018-09-23 20:53:35 +02:00
Martin Diehl
72b87b0a9b
better name
2018-09-23 19:58:43 +02:00
Martin Diehl
f701910c71
Merge branch '38-introduce-rudimentary-PETSc-based-FEM-solver' of magit1.mpie.de:damask/DAMASK into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-09-23 19:32:19 +02:00
Martin Diehl
ebef12e446
syntax errors in declaration
2018-09-23 19:04:17 +02:00
Martin Diehl
67483487ea
more descriptive names
2018-09-23 18:50:54 +02:00
Martin Diehl
854d99250c
mesh_element(1,:) only used for debug output
...
set to -1 at the moment to indicate that it is not used.
Re-implementation should be done for MSC.Marc and Abaqus only.
2018-09-23 18:42:49 +02:00
Martin Diehl
4b14cc5560
calcmode only needed for Abaqus and MSC.Marc
2018-09-23 18:06:18 +02:00
Martin Diehl
cf6d388a6b
consistent order of solver-specific functionality
2018-09-23 17:57:48 +02:00
Martin Diehl
5936397ae7
introducing better names
...
allows further simplifications as we do not store max and per elem
values any more for number of integration points and number of cell
nodes
2018-09-23 17:37:57 +02:00
Martin Diehl
2fe2c4ca45
leaner syntax with sourced allocation
2018-09-23 17:26:13 +02:00
Martin Diehl
51390b1acf
Nelems /= NcpElems only in case of Abaqus/MSC.Marc
2018-09-23 17:05:01 +02:00
Martin Diehl
5814e07021
simplified
2018-09-23 16:55:03 +02:00
Martin Diehl
e3f2ad34b2
not needed
2018-09-23 16:32:13 +02:00
Martin Diehl
bd60018513
not needed for spectral
2018-09-23 15:57:21 +02:00
Martin Diehl
c1b8854132
only needed for commercial solvers
2018-09-23 15:31:19 +02:00
Martin Diehl
4862aca340
grouping solver specific variables for better readability
2018-09-23 15:27:51 +02:00
Martin Diehl
db45b7615a
drop support for heterogeneous meshes
...
heterogeneous meshes are neither advisable nor typically used
2018-09-23 15:24:14 +02:00
Martin Diehl
519cd840bd
cleaning
2018-09-23 15:19:23 +02:00
Pratheek Shanthraj
109ed4308f
change in tabulation order. should now be working correctly
2018-09-23 03:14:23 +02:00
Martin Diehl
9eb318b9bf
Merge branch '38-introduce-rudimentary-PETSc-based-FEM-solver' of magit1.mpie.de:damask/DAMASK into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-09-22 14:54:58 +02:00
Pratheek Shanthraj
d9bdf53628
quadrature weights defined twice for some reason
2018-09-22 12:55:23 +02:00
Pratheek Shanthraj
ffcadc1d86
use displacements instead of coordinates
2018-09-22 12:54:45 +02:00
Pratheek Shanthraj
4a25520828
removed debug write statements
2018-09-22 12:53:33 +02:00
Pratheek Shanthraj
6c21c60c70
corrected nbasis
2018-09-22 12:53:02 +02:00
Pratheek Shanthraj
becd959913
dmda not needed
2018-09-22 12:52:29 +02:00
Pratheek Shanthraj
81ced06393
order --> degree
2018-09-22 12:51:02 +02:00
Pratheek Shanthraj
22c1436af6
working with petsc 3.10
2018-09-22 12:49:56 +02:00
Pratheek Shanthraj
c9a4609290
reading in BCs correctly
2018-09-22 12:48:57 +02:00
Martin Diehl
42bc1b468a
cleaning+adopting to petsc-3.10
2018-09-22 11:19:08 +02:00
Martin Diehl
a95e2b8083
Merge branch 'PETSc3.10' into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-09-22 10:36:31 +02:00
Martin Diehl
f590851a78
Merge branch '38-introduce-rudimentary-PETSc-based-FEM-solver' of magit1.mpie.de:damask/DAMASK into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-09-22 10:33:20 +02:00
Martin Diehl
80914470cc
most probably pInt and PetscInt is the same
2018-09-22 10:32:26 +02:00
Jaeyong Jung
4ddb8d843c
modified DMGetDefaultGlobalSection to DMGetGlobalSection
2018-09-21 14:08:44 +02:00
Jaeyong Jung
9150844c96
Merge branch '38-introduce-rudimentary-PETSc-based-FEM-solver' of magit1.mpie.de:damask/DAMASK into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-09-21 10:13:27 +02:00
Martin Diehl
180105d3ce
external statements will not work if interfaces exist
...
we're currently providing a patch for PETSc with interfaces for all
functions that we use
2018-09-21 08:25:35 +02:00
Martin Diehl
ff241ae4a6
only use required includes
2018-09-21 08:19:36 +02:00
Martin Diehl
bd9a02bfe2
creating hdf5 file to contain restart information
...
first quick and dirty implementation to get started. needs polishing and
renaming
2018-09-20 19:42:58 +02:00
Jaeyong Jung
abe5c1d825
Merge branch '38-introduce-rudimentary-PETSc-based-FEM-solver' of magit1.mpie.de:damask/DAMASK into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-09-20 09:57:30 +02:00
Martin Diehl
b775d5084b
more (temporarly) cleaning
2018-09-20 08:26:59 +02:00
Martin Diehl
1545a53ea9
ID do not exist at the moment
2018-09-20 08:00:39 +02:00
Martin Diehl
cc262ae198
Merge branch 'development' into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-09-20 07:29:19 +02:00
Martin Diehl
3e7b80a3ef
debug only available if compiled in debug mode
2018-09-20 07:27:12 +02:00
Martin Diehl
fc016bbc6e
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-09-20 07:16:03 +02:00
Martin Diehl
310ea62964
only print out the essential information
2018-09-20 06:58:31 +02:00
Martin Diehl
fcff6b908a
can be easily computed during post processing
2018-09-20 06:40:23 +02:00
Martin Diehl
901355d2ae
don't use unnecessarily long names
2018-09-20 06:27:53 +02:00
Martin Diehl
df0464c31b
use (import) debug variables only when needed
2018-09-20 06:24:03 +02:00
Martin Diehl
6aa4dd842a
define debug variables only if needed
2018-09-20 06:09:02 +02:00
Martin Diehl
1623a33b48
cleaning (mainly OMP FLUSh)
2018-09-19 21:46:26 +02:00
Martin Diehl
191ad9df09
not used
2018-09-19 20:40:35 +02:00
Martin Diehl
1bf1e13b46
some debug statements
2018-09-19 20:23:28 +02:00
Martin Diehl
e7e959af47
using explicit interface from PETScDT
2018-09-19 20:06:35 +02:00
Martin Diehl
11d4c28d88
flushes not needed + further cleaning
2018-09-19 19:45:57 +02:00
Martin Diehl
b84476e681
cleaning and debugging
2018-09-19 18:33:37 +02:00
Martin Diehl
c313dc1675
only read access
2018-09-19 17:04:12 +02:00
Martin Diehl
0bf64645a1
should be done by the plasticity laws (for the moment)
2018-09-19 14:22:35 +02:00
Martin Diehl
a8fb7d7ade
not needed
...
but I'm under the impression that the compiler removes such things anyway
2018-09-19 14:21:10 +02:00
Martin Diehl
20f0bee459
fallback dPdF not needed
...
save a lot of memory
2018-09-19 06:19:40 +02:00
Martin Diehl
18db5f5652
not needed
2018-09-18 05:20:08 +02:00
Martin Diehl
82cdc551a5
MPI Communicator only needed for PETSc 3.10.
...
thx to Jaeyong for figuring this out
2018-09-18 05:17:48 +02:00
Martin Diehl
245676728d
not needed
2018-09-17 23:28:00 +02:00
Jaeyong Jung
3beab556a2
runtime polar decompoistion error comes up
2018-09-17 15:41:58 +02:00
Martin Diehl
c1a9d3fbf6
names in agreement with the DAMASK paper
2018-09-17 08:23:23 +02:00
Martin Diehl
4a5e55cce9
Merge branch 'development' into 20-NewStyleDislotwin
2018-09-17 07:59:46 +02:00
Martin Diehl
e6fa3f3d35
correct stress in postResults
2018-09-16 22:57:50 +02:00
Martin Diehl
cea2fba063
don't repeat code that can be handled centrally
2018-09-16 22:38:57 +02:00
Martin Diehl
9570fb894a
correct names and no superflous conversions anymore
2018-09-16 22:31:55 +02:00
Martin Diehl
5f06a35900
Merge branch 'development' into 19-NewStylePhenopowerlaw
2018-09-16 22:03:05 +02:00
Martin Diehl
b54742cfa7
[skip sc] [skip ci] PETSc 3.10 is out
2018-09-15 20:09:27 +02:00
Martin Diehl
075d2e51b2
debugging, still crashing
2018-09-15 16:19:41 +02:00
Martin Diehl
3b1817a954
cleaning, only having mech at the moment
2018-09-15 15:14:48 +02:00
Martin Diehl
934ca56811
clean exit if required files not found
2018-09-15 14:54:56 +02:00
Martin Diehl
a92b945e3f
does not crash anymore
...
xx is nc, the number of components
needs to be 3 in the current case
2018-09-15 11:40:07 +02:00
Martin Diehl
45dc59b1d8
kinetics also for trans
2018-09-15 10:43:05 +02:00
Martin Diehl
5e33900664
Merge branch 'development' into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-09-15 08:38:30 +02:00
Martin Diehl
ed570f0fe8
use of kinetics avoids different calculation of shearrates
2018-09-15 08:24:12 +02:00
Martin Diehl
3526c0a757
index missing, test failed
2018-09-15 07:49:14 +02:00
Martin Diehl
acc3294909
more verbose initialization
2018-09-15 06:39:02 +02:00
Martin Diehl
25a38ad438
initialization broken
...
- need to consider case of no twin/slip active
- state0 needs to be initialized
2018-09-15 05:45:10 +02:00
Martin Diehl
a3750738c0
type mismatch for MSC.Marc
2018-09-14 12:20:37 +02:00
Martin Diehl
9c5ca0e542
use of kinetics simplifies code
2018-09-14 12:17:35 +02:00
Martin Diehl
bd0b561772
better readable
2018-09-14 11:56:36 +02:00
Martin Diehl
9f16cefd9f
renaming in accordance with the DAMASK paper
2018-09-14 11:17:14 +02:00
Martin Diehl
af32b3d85b
reordering for easier overview
2018-09-14 10:51:44 +02:00
Martin Diehl
3dd47eade9
just renaming, better readable...
2018-09-14 10:18:44 +02:00
Martin Diehl
07b0ddf711
polishing
2018-09-14 10:08:36 +02:00
Martin Diehl
80a2062c4b
Merge remote-tracking branch 'origin/20-NewStyleDislotwin' into 20-NewStyleDislotwin
2018-09-14 10:00:11 +02:00
Martin Diehl
18858301d5
using consistently i as first running index and j as second
2018-09-14 09:59:04 +02:00
Martin Diehl
55511782e8
incomplete transition/renaming
...
caused segmentation fault
2018-09-14 08:19:39 +02:00
Martin Diehl
f98243e4ac
Merge branch 'development' into 43-wrong-stress-in-dotstate-and-deltastate
2018-09-14 05:39:15 +02:00
Martin Diehl
989393e6e8
should not differ from development
2018-09-14 05:37:17 +02:00
Martin Diehl
e2a66f6ddb
should not differ from development
2018-09-14 05:33:39 +02:00
Martin Diehl
b165e48d9f
non-Schmid contribution were ignored in test
...
reverted, just to ensure that the results have not changed: Never
polishing and change physics at the same time
2018-09-14 05:21:05 +02:00
Martin Diehl
cc20b044a6
forgot to allocate array now used internally only
2018-09-13 09:27:54 +02:00
Martin Diehl
ad3994c484
allocate all dependentStates/microstructures
2018-09-13 06:48:06 +02:00
Martin Diehl
e7f7fa74b7
don't do calculation of dependentStates/microstructure manually
2018-09-13 06:20:58 +02:00
Martin Diehl
8baa4f9c81
need to avoid array access out of bounds
2018-09-13 06:12:32 +02:00
Martin Diehl
ae8efe2755
Merge branch 'development' into 44-column-major-access-to-interaction-matrices
2018-09-13 05:55:59 +02:00
Martin Diehl
8ecb019566
real must be pReal for MSC.Marc
2018-09-13 05:51:58 +02:00
Martin Diehl
4266872965
array out of bounds was possible for Ntwin =0
...
in postResults(c+1:c+prm%totalNtwin) c+1 will be out of bounds if
c is already size(postResults)
2018-09-13 05:41:44 +02:00
Martin Diehl
66edad1cf8
avoid code doubling
...
shear rates needed multiple times. For now, just introduced but not
used. Will become active once the test passes again
2018-09-12 22:07:59 +02:00
Martin Diehl
3068caa9a3
[skip sc] resolved stress not needed
...
using kinetics_xxx as in disloUCLA
compiles on gfortran but pre-receive hook with intel compiler (MSC.Marc)
fails
2018-09-12 21:05:40 +02:00
Martin Diehl
ac3625afba
array access out of bounds
2018-09-12 20:42:57 +02:00
Martin Diehl
4e68d049a8
WIP: calculating rotated stiffness matrices for transformation
2018-09-12 20:37:55 +02:00
Martin Diehl
b95174a8b7
reducing global variables
2018-09-12 16:16:31 +02:00
Martin Diehl
192bb6453d
cleaning
2018-09-12 15:57:54 +02:00
Martin Diehl
b753a86d13
cleaning
2018-09-12 15:27:47 +02:00
Martin Diehl
0a621ece8b
WIP: polishing and implementing remaining functions
2018-09-12 14:29:22 +02:00
Martin Diehl
f1f8922ab8
more meaningful error messages
2018-09-12 14:03:45 +02:00
Martin Diehl
d7023096ad
introduced kinetics, not in use yet
2018-09-12 13:25:18 +02:00
Martin Diehl
edebe4d1ed
vectorization error
2018-09-12 12:26:33 +02:00
Martin Diehl
c9b5b3fb7b
should be always set, even if no twinning is enabled
...
fixes NaN in dotState
2018-09-12 12:08:03 +02:00
Martin Diehl
253f318e55
preventing NaN in dotState
2018-09-12 11:47:35 +02:00
Martin Diehl
2337dde525
cleaning
2018-09-12 11:26:59 +02:00
Martin Diehl
8a406150f8
polishing
2018-09-12 10:06:04 +02:00
Martin Diehl
3352611267
cleaning
2018-09-12 09:59:09 +02:00
Martin Diehl
01196ca36b
Merge branch 'development' into 19-NewStylePhenopowerlaw
2018-09-12 09:50:24 +02:00
Martin Diehl
c9208315f5
avoid repeated calculations
...
does not save so much here, but avoids having inconsistent calculation
(e.g. nonSchmid effects) and serves as a template for more complex
models
2018-09-12 09:23:11 +02:00
Martin Diehl
ad7c099fe3
Merge branch 'development' into 20-NewStyleDislotwin
2018-09-11 17:03:05 +02:00
Martin Diehl
ce6e6679d5
causes array access out of bounds, needs further checks
2018-09-11 12:20:05 +02:00
Martin Diehl
51a6f4b990
Merge branch 'development' into 20-NewStyleDislotwin
2018-09-11 02:39:53 +02:00
Martin Diehl
3d49c70dbc
commenting
2018-09-08 20:15:17 +02:00
Philip Eisenlohr
9be2c084e4
Merge branch '39-simplify-obscure-numerics-integration-mode' into 'development'
...
Resolve "simplify obscure numerics integration mode"
Closes #39
See merge request damask/DAMASK!36
2018-09-07 19:34:35 +02:00
Philip Eisenlohr
7e330f1249
Merge branch '37-include-hdf5-routines' into 'development'
...
Resolve "Include HDF5 routines"
Closes #37
See merge request damask/DAMASK!37
2018-09-07 19:18:07 +02:00
Satyapriya Gupta
a16454f1f9
fixed missing propagation of dependent state variables
2018-09-06 16:48:28 -04:00
Martin Diehl
e70efd2608
inform about debug versions indenpendently of the solver
2018-09-06 16:49:39 +02:00
Martin Diehl
8324962123
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-09-06 16:26:46 +02:00
Martin Diehl
8b3d30fdd4
don't know what happend to the copyright note
2018-09-05 23:47:01 +02:00
Martin Diehl
0f084a9844
cleaning output
2018-09-05 23:39:54 +02:00
Martin Diehl
918de24f83
using microstructure instead of 'dependent state'
2018-09-05 23:00:46 +02:00
Martin Diehl
1f9d268e29
renaming
2018-09-05 16:15:57 +02:00
Martin Diehl
97977f4940
all parameters should be stored in the constitutive laws
...
no need to know the 'phase number' anymore
2018-09-05 15:54:19 +02:00
Martin Diehl
038508aa11
'dependent states' do not need to be accessible globally
2018-09-05 11:07:00 +02:00
Martin Diehl
4fbe5811a3
too much whitespace
2018-09-01 11:08:38 +02:00
Martin Diehl
e9f738fade
cleaning + renaming
2018-09-01 10:58:50 +02:00
Martin Diehl
cf65aae92a
correct names
2018-09-01 06:49:03 +02:00
Martin Diehl
0b94b8085d
still just polishing
2018-08-31 18:17:35 +02:00
Philip Eisenlohr
6a64637ac3
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-08-31 11:23:25 -04:00
Philip Eisenlohr
3e81110830
tell user whether DEBUG version or not
2018-08-31 11:22:59 -04:00
Martin Diehl
34d1c011c7
long name does not help here
2018-08-31 16:38:17 +02:00
Martin Diehl
ea1fd621aa
more cleaning
2018-08-31 16:36:19 +02:00
Martin Diehl
7a67922c5f
general polishing
2018-08-31 16:08:01 +02:00
Martin Diehl
912926e604
more improvements
2018-08-31 15:07:22 +02:00
Martin Diehl
c0e7050867
more loops simplified
2018-08-31 14:54:33 +02:00
Martin Diehl
60e60e211c
also simplifying twin loops
2018-08-31 14:33:42 +02:00
Martin Diehl
dc91016729
simpler loop structures for better readability
2018-08-31 13:38:07 +02:00
Martin Diehl
42745c66f0
nested loops do not improve readability
2018-08-31 13:32:24 +02:00
Martin Diehl
ba65044ff5
looping over systems simplifies loop structure
2018-08-31 13:19:28 +02:00
Martin Diehl
3431086fdb
storing data per slip system simplifies calculation
2018-08-31 13:14:16 +02:00
Martin Diehl
c9bb051261
fill local Schmid Matrices
2018-08-31 11:52:43 +02:00
Martin Diehl
54fe0afb16
indentation wrong
2018-08-31 11:42:54 +02:00
Martin Diehl
fc4253ed57
conversion from 3x3 to 6 of stress in constitutive will be remove later
2018-08-31 11:37:14 +02:00
Martin Diehl
3b88096cd7
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-08-31 11:09:52 +02:00
Martin Diehl
5330fb31ad
more verbose error message
2018-08-31 10:19:37 +02:00
Martin Diehl
f028e05298
guessing was not correctly set
2018-08-31 10:14:33 +02:00
Martin Diehl
3cb279b083
one more loop not needed
2018-08-31 10:03:15 +02:00
Martin Diehl
017563e061
one more loop not needed
2018-08-31 09:14:16 +02:00
Martin Diehl
c0ac05aa88
no need to use two loops
2018-08-31 09:12:50 +02:00
Martin Diehl
732022d4f0
functions from IO does not add any value here
...
only prevents the automated selection of a new unit
2018-08-31 08:58:13 +02:00
Martin Diehl
d18d74ae31
mechanics solver type is set only once
...
as all mech routines have the same interface, this leaner syntax can be
used
2018-08-31 08:20:23 +02:00
Martin Diehl
0b049fc843
Merge branch 'development' into 42-new-coding-style-for-homogenization
2018-08-31 06:36:12 +02:00
Martin Diehl
26d18257d2
requested output is stored in prm%outputID
2018-08-31 06:35:30 +02:00
Martin Diehl
91a3b4ed69
almost done
2018-08-31 01:24:28 +02:00
Martin Diehl
6b45afa72f
using parameters from config.f90
2018-08-30 17:51:31 +02:00
Martin Diehl
555ca8e5f1
debug statements not needed anymore
2018-08-30 17:21:22 +02:00
Martin Diehl
62413e8187
now compiles again
2018-08-30 17:18:56 +02:00
Martin Diehl
5f97ae9080
not used at the moment
2018-08-30 16:39:18 +02:00
Martin Diehl
cbecb71000
[skip sc] using explicit interface helps to get correct type
...
still some issues with different ranks and datatypes
2018-08-30 13:10:18 +02:00
Martin Diehl
b930c68887
PETSc 3.9 style (no SIGSEGVs anymore)
2018-08-30 12:38:44 +02:00
Martin Diehl
8b6bef159e
a little bit more verbose
2018-08-30 12:38:29 +02:00
Martin Diehl
3254ae79a5
collecting questions to Pratheek
2018-08-30 12:37:47 +02:00
Martin Diehl
d180af494e
adopted to new PETSc fortran interface
2018-08-30 11:55:13 +02:00
Martin Diehl
1f9a614388
no need to have long prefixes for local variables and functions
2018-08-30 10:50:52 +02:00
Martin Diehl
88d776cad6
consistent renames
2018-08-30 10:46:30 +02:00
Martin Diehl
29e55d20fb
message better to understand and giving error instead of SIGSEGV
2018-08-30 09:42:45 +02:00
Martin Diehl
b96e01c128
Merge branch 'development' into lattice-builds-interaction
2018-08-30 05:41:40 +02:00
Martin Diehl
105c4077ad
Merge branch 'development' into 19-NewStylePhenopowerlaw
2018-08-30 00:58:54 +02:00
Martin Diehl
b163a8aaa0
Non-schmid activated again
...
internally, no need for long name plastic_phenopowerlaw_postResults
2018-08-30 00:57:15 +02:00
Martin Diehl
bcff95ddf8
using new values for initialization
2018-08-29 13:52:26 +02:00
Martin Diehl
0b2dd86bbf
handling cluster orientation
2018-08-29 13:36:46 +02:00
Martin Diehl
b884349e7b
only renaming
...
3333 not needed for dX_dY if X and Y are 3x3 tensors
PK2 stress is S not T according to the DAMASK paper
2018-08-29 13:16:37 +02:00
Martin Diehl
6d28883db5
starting to include non-Schmid terms (again)
2018-08-29 12:07:51 +02:00
Martin Diehl
922273f230
does not make sense to store and use the 6-vector version of the Schmid
...
matrix
2018-08-29 12:03:31 +02:00
Martin Diehl
baeb449e07
WIP: debugging
2018-08-29 11:44:39 +02:00
Martin Diehl
ef506b801e
fixed two memory faults (unallocated and wrong index)
2018-08-29 09:52:22 +02:00
Martin Diehl
f2acc14882
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-08-29 08:27:15 +02:00
Martin Diehl
680c9e11d4
segmentation fault in cause of empty list
2018-08-29 08:26:28 +02:00
Martin Diehl
b5e4ad247d
LpAndItsTangent naming convention following DAMASK paper
2018-08-28 15:40:24 +02:00
Martin Diehl
fc3ce54667
return more than one tangent
2018-08-28 15:07:39 +02:00
Martin Diehl
15e2d4a7cd
Tstar is S
...
according to the DAMASK paper
2018-08-28 15:02:58 +02:00
Martin Diehl
4c14f988a3
rename according to paper
2018-08-28 14:58:17 +02:00
Martin Diehl
63b939489b
reading in parameters using new style in parallel
2018-08-28 13:15:38 +02:00
Martin Diehl
f9214c8e1b
should not be handled by the individual model
...
avg(P) and avg(F) exist independently of RGC
2018-08-28 13:02:44 +02:00
Martin Diehl
7f00082d60
parameters easier to identify
2018-08-28 12:57:22 +02:00
Martin Diehl
06e7156351
did not work for values <0
2018-08-28 12:41:21 +02:00
Martin Diehl
42f8b0a063
labels of slip and twin systems for more self-explanatory output
2018-08-28 12:27:56 +02:00
Martin Diehl
94695f773e
more verbose and works for arbitrary precision
2018-08-28 07:47:05 +02:00
Martin Diehl
4112da183e
Merge branch 'development' into 42-new-coding-style-for-homogenization
2018-08-28 06:28:13 +02:00
Martin Diehl
9b3ddcd2c2
more explicit
2018-08-28 06:27:38 +02:00
Martin Diehl
331eaabaa4
standard style from plasticity
...
using derived type for parameters
explicit "use" to see dependencies
2018-08-28 06:24:26 +02:00
Franz Roters
38fd517c44
Merge remote-tracking branch 'remotes/origin/36-faster-file-handling-for-material-config-use-stream-access-instead-of-line-wise-reading' into development
2018-08-27 15:13:16 +02:00
Martin Diehl
38a855243d
Merge branch 'development' into 43-wrong-stress-in-dotstate-and-deltastate
2018-08-27 06:31:45 +02:00
Martin Diehl
dd81fe806d
naming was not completely corrected
2018-08-27 06:30:55 +02:00
Martin Diehl
11a7103675
more explicit
2018-08-27 00:00:21 +02:00
Martin Diehl
e6408e0ce3
corrected careless changes
...
using unallocated pointer and asigning wrong label
2018-08-26 21:05:59 +02:00
Martin Diehl
b11666effa
Merge branch 'development' into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-08-26 13:59:50 +02:00
Martin Diehl
a5f139b786
unused variables
2018-08-26 11:53:18 +02:00
Martin Diehl
241b2ade8b
more portable way to define PI
...
https://software.intel.com/en-us/forums/intel-visual-fortran-compiler-for-windows/topic/490432
and compiler will not complain about truncation
2018-08-26 11:40:38 +02:00
Martin Diehl
0041d21777
output was screwed up
2018-08-26 06:20:48 +02:00
Martin Diehl
4f1becb503
cleaning and fixing bugs (wrong indices)
2018-08-25 21:27:14 +02:00
Martin Diehl
f458de82fa
simplified
2018-08-25 21:13:32 +02:00
Martin Diehl
3ff7c9c0eb
extendend to Schmid_twin
2018-08-25 20:55:15 +02:00
Martin Diehl
fbaac484ea
extended to postResults
2018-08-25 20:32:01 +02:00
Martin Diehl
780699740d
internal (non)-Schmid definitions also in dotState
2018-08-25 20:19:02 +02:00
Martin Diehl
ba9330c8dd
early rename
2018-08-25 20:08:46 +02:00
Martin Diehl
d51b00c2d3
no need to loop over families and systems independently. no need for 3333
2018-08-25 20:06:19 +02:00
Martin Diehl
d055ef665a
using only internal schmid and non schmid matrices
2018-08-25 20:01:39 +02:00
Martin Diehl
cc5d04ff2a
newly introduced parameters for non-schmid were confusing
2018-08-25 19:28:09 +02:00
Martin Diehl
5b3bed3e48
functionality to calculate SchmidMatrices
2018-08-25 18:32:55 +02:00
Martin Diehl
d3db0d816c
ordered by dimension and set private where possible
2018-08-25 17:29:20 +02:00
Martin Diehl
1c76e156d9
new symmetry handling will be handled in a different branch
2018-08-25 17:19:48 +02:00
Martin Diehl
bde242902b
Merge branch 'development' into 43-wrong-stress-in-dotstate-and-deltastate
2018-08-25 16:02:12 +02:00
Martin Diehl
b8b5bac684
dotState and deltaState parse Mstar instead of Tstar
...
requires to parse in Fi
2018-08-25 15:59:34 +02:00
Martin Diehl
22a232ad08
bug: memory access out of bounds
...
introduced when moderninzing reading in of parameters
2018-08-25 15:50:43 +02:00
Martin Diehl
7e002d1bfb
missing rename
2018-08-25 14:42:44 +02:00
Martin Diehl
e46605f0ef
forall is deprecated
2018-08-25 14:36:21 +02:00
Martin Diehl
17d88184a7
typos fixed
2018-08-25 14:35:49 +02:00
Martin Diehl
35688a6acf
temp33 not needed (compiler should be smart enough)
2018-08-25 14:29:03 +02:00
Martin Diehl
953fff79ac
prepared LpAndItsTangent to remove superflous forward-backward
...
conversion
updating the individual plastic laws will be done in the respective
branches
2018-08-25 14:21:06 +02:00
Martin Diehl
d585deee7e
unified syntax
...
_v for vector representation of tensor conflicts with 3333 suffix for
4-th order tensors.
General idea: Mark symmetric second order tensors in vector notation with
'6' and fourth order tensors in second order matrix notation with '99'.
Append nothing for 'natural' representation, i.e. F and NOT F33, Lp_dTstar and NOT
Lp_dTstar3333
2018-08-25 14:07:23 +02:00
Martin Diehl
025cbddd00
generic function and interfaces for slipSlip, twinTwin, transTrans
2018-08-25 13:08:32 +02:00
Martin Diehl
6800a5a6f6
just adopting naming convention
2018-08-25 11:07:56 +02:00
Martin Diehl
49ae38d0f9
using parameters from linked list and removed output
...
homogenization models should only provide model specific output in
accordance with http://dx.doi.org/10.1007/s40192-017-0084-5
2018-08-25 11:03:43 +02:00
Martin Diehl
7ecb7689f1
Intel compiler failed with SIGSEV
...
derived types, pointers, finalize .... altogether seems to bring both Compilers to their limits.
I cannot see what was wrong before, but now it works and might be a little faster
2018-08-23 00:13:57 +02:00
Martin Diehl
271b9ba76b
intersting note ...
2018-08-22 23:28:47 +02:00
Martin Diehl
8c5f3d4e07
only needed once
2018-08-22 22:52:12 +02:00
Martin Diehl
51b4ef319a
HDF5 was not compiled due to missing dependency
2018-08-22 20:12:16 +02:00
Martin Diehl
dc596e6789
zero termination does not work
2018-08-22 18:09:17 +02:00
Martin Diehl
1a943df97e
small flaws
2018-08-22 17:52:00 +02:00
Martin Diehl
ab45818d51
seems to work now
...
anyway, nicer code
2018-08-22 14:30:51 +02:00
Martin Diehl
1b5623ad6c
avoid out of bound access, removed unneeded stuff
2018-08-22 14:22:37 +02:00
Martin Diehl
037ab3d081
getTag works now for tags with the same start and close tag
...
needed for '/echo/'. Still suggest to rather use /echo\
2018-08-22 13:10:59 +02:00
Martin Diehl
52088691d1
improved linked list and fixed solution for strange bug
...
Bug: Using automated LHS re-allocation for a string array that with global scope seems to cause trouble
Hence, "parse_file" works with a local string and assings only once to it
Linked_List: Now storing data in the list head also and last element is always empty.
Finalize allows simple handling of deallocation
2018-08-22 11:51:23 +02:00
Martin Diehl
52002f654e
to converge at one point to one (or two) string lenth values
2018-08-22 09:14:16 +02:00
Martin Diehl
0068dba40f
Merge branch 'development' into 36-faster-file-handling-for-material-config-use-stream-access-instead-of-line-wise-reading
2018-08-22 09:11:33 +02:00
Martin Diehl
4d432e5462
Merge branch 'development' into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-08-21 23:03:56 +02:00
Martin Diehl
0fa902100c
Merge branch 'development' into 39-simplify-obscure-numerics-integration-mode
2018-08-21 22:59:57 +02:00
Martin Diehl
61070fd908
Merge branch 'development' into 37-include-hdf5-routines
2018-08-21 22:46:30 +02:00
Martin Diehl
8de321382f
using final seems to be the better approach
...
http://www.training.prace-ri.eu/uploads/tx_pracetmo/AdvFTN_handout.pdf
still needs in-depth analysis, even though I cannot reproduce 'terminally ill'
2018-08-21 08:14:59 +02:00
Martin Diehl
d146417abe
hot fix for terminally ill with gfortran 7.3
...
might be a bug in the compiler or in the linked list.
waste some memory at the moment...
check linked list carefully before enabling again and blaming gfortran
2018-08-21 07:41:10 +02:00
Martin Diehl
465d950ab1
gfortran 7.3 and optimized code still gives 'terminally ill'
2018-08-21 07:09:50 +02:00
Martin Diehl
60f56255e4
[skip sc] [skip ci] simplified interfacing
2018-08-20 23:14:49 +02:00
Martin Diehl
69ad600916
more explicit file opening
...
still having trouble with Gfortran 7.3 and no numerics.config
2018-08-20 22:36:55 +02:00
Philip Eisenlohr
87a8a9536b
absTol for divergence was way too tight
2018-08-20 16:34:43 -04:00
Martin Diehl
c78396dd78
randomized FILEUNIT
...
matrix inversion error when numerics.config is not there or not present
very strange
2018-08-20 22:11:53 +02:00
Martin Diehl
e47677738a
more verbose error
2018-08-20 17:57:15 +02:00
Martin Diehl
031c59954f
working directory is changed automatically
2018-08-20 17:45:34 +02:00
Martin Diehl
29cd8fa5e2
Merge branch 'development' into 36-include-hdf5-routines
2018-08-20 17:35:32 +02:00
Martin Diehl
dbed7056e5
[skip sc] cleaning
2018-08-20 17:08:33 +02:00
Martin Diehl
93562d5142
mapping of elements etc not needed for PETSc-based FEM and spectral
...
solvers
2018-08-20 16:42:47 +02:00
Martin Diehl
e3e905938e
all elements are CP elements
2018-08-20 16:23:05 +02:00
Martin Diehl
3e4c878304
using shared interface for spectral and FEM solver
...
group_scalar seems inappropriate as integers are also scalars. renamed
to group_float (is actually usually of double precision).
think about better name, types should have a t prefix. tgroupFloat?
2018-08-20 16:09:40 +02:00
Martin Diehl
b8d56ae320
unfified interface for spectral and FEM solver
...
Note: extension to load case and geometry is not added automatically
anymore!
2018-08-20 15:59:13 +02:00
Martin Diehl
f29a5b3df3
not used at all
2018-08-20 13:15:44 +02:00
Martin Diehl
543c49426f
Merge branch 'development' into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-08-20 09:09:15 +02:00
Martin Diehl
3a097d1099
Merge branch 'development' into 36-faster-file-handling-for-material-config-use-stream-access-instead-of-line-wise-reading
2018-08-20 09:04:28 +02:00
Martin Diehl
8fb780ab42
now compiles with gfortran
2018-08-20 09:01:13 +02:00
Martin Diehl
d6c69e3a18
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-08-18 23:17:07 +02:00
Martin Diehl
85358bae1b
copy and paste error
2018-08-18 21:30:57 +02:00
Martin Diehl
97659c359a
MPI write only in case that PETSc is used
2018-08-18 21:01:25 +02:00
Martin Diehl
3b3e0bc068
forgotten use
2018-08-18 17:11:50 +02:00
Martin Diehl
67314fc204
leftover from reducing numerics_integrator to scalar
...
postponed simplification to scalar because of heavy modification in plastic
constitutive laws
2018-08-18 16:42:32 +02:00
Martin Diehl
2618593568
:Merge branch 'development' into Integration-Mode-Confusion
2018-08-18 16:37:25 +02:00
Martin Diehl
f7c20d74af
compiles now, but most likely does not work
2018-08-18 15:58:42 +02:00
Martin Diehl
0d8f17cbe6
adjusting to PETSc 3.9.x
2018-08-18 14:05:57 +02:00
Philip Eisenlohr
c03ea8f5c7
added forgotten C66(2,2) read-in from parameter database
2018-08-17 15:43:09 -04:00
Philip Eisenlohr
62d2806fe1
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-08-17 14:59:00 -04:00
Martin Diehl
d4bcfae82b
WIP: adopting to PETSc 3.9.x and modifications in development branch
2018-08-17 12:00:44 +02:00
Martin Diehl
821860987c
copied existing files
2018-08-17 00:14:25 +02:00
Martin Diehl
26b08bc21b
Merge branch 'development' into 36-faster-file-handling-for-material-config-use-stream-access-instead-of-line-wise-reading
2018-08-16 23:39:14 +02:00
Martin Diehl
b64b5b82de
integrationMode was always 1, very obscure code
...
was most likely the leftover of old functionality
2018-08-16 23:33:39 +02:00
Philip Eisenlohr
1493b33a03
fixed buggy state width assignment (endindex += %nSlip corrected to %totalNslip)
2018-08-14 01:07:06 +02:00
Philip Eisenlohr
66e68151ea
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-08-13 18:28:59 -04:00
Philip Eisenlohr
ba215ed9ea
small polish
2018-08-13 18:27:51 -04:00
Martin Diehl
5337fb1229
old runtime debugging leftovers
2018-08-05 11:08:17 +02:00
Martin Diehl
c0adb61350
outdated description
2018-08-05 10:41:01 +02:00
Martin Diehl
60fc47ca2a
unused import
2018-08-05 09:41:27 +02:00
Martin Diehl
fd578ab8eb
Merge branch 'development' into 19-NewStylePhenopowerlaw
2018-08-05 07:38:09 +02:00
Martin Diehl
386c3e7797
Merge branch 'development' into 36-faster-file-handling-for-material-config-use-stream-access-instead-of-line-wise-reading
2018-08-05 06:04:33 +02:00
Martin Diehl
661636d83e
default values needed to suppress error in case of missing key
2018-08-04 22:53:56 +02:00
Martin Diehl
f522f79703
Merge branch 'development' into 35-new-style-of-reading-in-parameters-in-lattice
2018-08-04 22:42:14 +02:00
Martin Diehl
6a4260ac07
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-08-04 22:29:52 +02:00
Martin Diehl
9360c171a1
polished for merge
...
access to array(n+m:) is safe for array of size n with m>1
2018-08-04 19:39:50 +02:00
Martin Diehl
20d1264d07
small improvements
...
default case of error handling, checking for recursion limit, some
comments to also understand it later
2018-08-04 13:58:01 +02:00
Philip Eisenlohr
8edeeaf213
renaming of p_vec and friends to names based on "group"
2018-08-03 19:39:14 -04:00
Martin Diehl
8a32ed8747
calculating non-schmid related data only once
2018-08-03 23:03:51 +02:00
Martin Diehl
8c07825d7e
Merge branch 'development' into 36-faster-file-handling-for-material-config-use-stream-access-instead-of-line-wise-reading
2018-08-03 22:21:37 +02:00
Martin Diehl
5b409fd6f7
type mismatch
2018-08-03 16:16:50 +02:00
Martin Diehl
f0645d84f5
not a target anymore
2018-08-03 15:10:38 +02:00
Martin Diehl
5cf2973715
naming reflect change from Piola-Kirchhoff to Mandel stress
...
not using the symmetrized stress anymore to avoid handling of
symmetrized Schmid matrizes. The time saved when calculating the double
contraction is probably anyway lost during the conversion from (3,3) to
(6) of Mstar
2018-08-03 14:53:40 +02:00
Martin Diehl
b8e8193001
many imported functions not used anymore, moving parameters to structure
2018-08-03 11:10:22 +02:00
Martin Diehl
38183fdbf6
Merge branch 'development' into 36-faster-file-handling-for-material-config-use-stream-access-instead-of-line-wise-reading
2018-08-03 08:57:44 +02:00
Martin Diehl
dc6855261a
number of helper functions for HDF5
2018-08-03 08:53:37 +02:00
Martin Diehl
c2f9709567
leaner, most of it also supported by older gcc
2018-08-03 08:42:26 +02:00
Martin Diehl
af931a78c7
finalizing reading in
2018-08-03 08:09:28 +02:00
Martin Diehl
68c476f690
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-08-03 07:56:02 +02:00
Martin Diehl
0a6bde70c5
Merge branch 'development' into 19-NewStylePhenopowerlaw
2018-08-03 06:56:43 +02:00
Martin Diehl
5d09e98e02
checking with Philip
2018-08-03 01:14:18 +02:00
Martin Diehl
3cceda5fe3
Merge branch 'development' into 31_consistent_wd
2018-08-03 00:08:16 +02:00
Pratheek Shanthraj
bc3f6ae97c
missing bits in getStrings function
2018-07-31 23:15:44 +02:00
Martin Diehl
a908e66306
WIP: reading in new style
2018-07-30 11:58:59 +02:00
Martin Diehl
ed97afb51c
WIP: inherent definition instead of check, using new reading in functionality
2018-07-29 21:18:04 +02:00
Martin Diehl
5eed7d06e9
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-07-29 20:29:55 +02:00
Martin Diehl
449449b500
does the same as numpy.clip
2018-07-28 01:31:02 +02:00
Martin Diehl
1c83f841aa
don't depend on any shear band parameters if sbVelocity is switched of
2018-07-20 17:52:15 +02:00
Martin Diehl
e7fd445816
WIP: Simplifying
...
using maxval gives inherently correct values and a consistency check is not needed any more
using config_phase capabilities to read in lattice related data
2018-07-20 14:13:13 +02:00
Martin Diehl
8554279351
outputID is not of default integer type (causes trouble with MSC.Marc)
2018-07-20 04:02:49 +02:00
Sharan Roongta
3f285d3a3f
polishing: prm replaces param(instance)
2018-07-17 19:29:39 +02:00
Martin Diehl
a3037a6064
cmpftwin was not read in
2018-07-17 18:02:55 +02:00
Martin Diehl
b3a6aabdc5
cthresholdtwin was not read in!
...
polished and commented
2018-07-17 12:32:57 +02:00
Martin Diehl
edcf97ea59
associate: clearer code and no performance drawbacks
2018-07-16 22:02:31 +02:00
Martin Diehl
df42e9ed88
using associate avoids performance issues related to pointers
2018-07-16 21:37:26 +02:00
Martin Diehl
e0e03827c5
simplified
2018-07-16 21:14:33 +02:00
Martin Diehl
fb1265db3d
checking for existing file
2018-07-16 20:08:31 +02:00
Martin Diehl
fc0499a9bb
bug fixes and polishing
...
bugs: wrong bracket in calculation of threshold_stress_twin
3**b_twin instead of 3*b_twin in calculation of threshold_stress_twin
ph instead of instance in calculation of mfp_twin
2018-07-16 19:53:26 +02:00
Martin Diehl
a1ad18c88a
256 characters is enough for material.config
...
larger values waste memory and decrease readability.
Still need to discuss how geom files are handled, for them longer limits make sense
2018-07-16 12:08:22 +02:00
Martin Diehl
7f05bf9c0a
leftover variables not needed anymore
2018-07-16 11:54:46 +02:00
Martin Diehl
b26c4a39ef
store raw material.config
2018-07-16 11:42:58 +02:00
Martin Diehl
3f7a1d1c07
function to read and store complete text file
...
reading as stream avoids costly repeated call to 'read'. Requires of course more memory, but that should be fine
also, recursion case ('{}') is internally handled. Old recursive was error prone and buggy when rewining (see 'reset' option)
2018-07-16 11:40:42 +02:00
Martin Diehl
bc5fcf2c14
leftover variable from cleaning
2018-07-16 10:51:42 +02:00
Martin Diehl
84c29bb255
Merge branch 'development' into 19-NewStyleDislotwin
2018-07-13 11:54:42 +02:00
Sharan Roongta
0c7c77987b
typo in plastic dislotwin
2018-07-12 23:22:55 +02:00
Sharan Roongta
eb7860ba22
Tolerances were different from original
2018-07-12 18:30:30 +02:00