9e926f1545
centralize test for error
2020-04-01 07:41:55 +02:00
570c042100
update of dependent state always in conjunction with state integration
2020-04-01 07:18:37 +02:00
54e3455bd4
RKCK45 store intermediate state only per point
2020-04-01 07:10:02 +02:00
6bfa51c307
LHS allocation does not work for pointers
2020-04-01 07:02:23 +02:00
e46220cd8a
OMP bugfix for FPI integrator, memory-efficient RK4
2020-03-31 22:02:53 +02:00
9c95ce36f4
automatic LHS (re)-allocation does not work for pointers
...
group_float has pointers, not allocatables
2020-03-31 20:57:09 +02:00
20a3ebecea
Merge branch 'state-integration-one-loop' into memory-efficient-state-storage
2020-03-31 20:29:01 +02:00
e810b6b23b
not needed here
2020-03-31 12:39:23 +02:00
00be291fa0
todo will be reset after state integration
2020-03-31 12:14:54 +02:00
2de4b87c61
bugfix for FPI: loop over sources!
...
no memory waste for adaptive Euler
2020-03-31 11:24:13 +02:00
5b29af8473
clearer
2020-03-31 10:16:14 +02:00
505c1432b1
no need for temp storage
2020-03-31 10:12:25 +02:00
2cc0c746d3
no variable needed
2020-03-31 10:07:01 +02:00
28dadb4422
no need to check multiple times
2020-03-31 09:49:41 +02:00
604bcd1229
memory efficient FPI state integrator
2020-03-31 09:38:08 +02:00
18f60a94a9
one variable suffices
2020-03-31 09:04:08 +02:00
c0d03a8c26
not needed
2020-03-31 08:37:34 +02:00
a7b602f4ca
Merge branch 'development' into state-integration-one-loop
2020-03-31 08:35:07 +02:00
9507101969
can be done in one loop
2020-03-31 08:33:35 +02:00
5f2759c03c
no need to check for converged
2020-03-31 08:23:38 +02:00
2080749189
does not change during the iteration
2020-03-31 08:09:41 +02:00
369682aa03
Merge branch 'fix-RK4' into 'development'
...
Fix rk4
See merge request damask/DAMASK!150
2020-03-30 20:05:09 +02:00
81c8ffe5c1
integrate into the same loop
2020-03-30 14:16:23 +02:00
28cf06ca98
extra loop not needed
2020-03-30 14:13:04 +02:00
cff62219b6
Merge branch 'development' into adaptive-euler-substate0
2020-03-27 14:43:55 +01:00
75fff766e1
not needed anymore
2020-03-26 22:02:28 +01:00
fdbbb94aac
takeover from RKCK45
2020-03-26 21:46:28 +01:00
c490b4bea4
standard names
2020-03-26 20:19:05 +01:00
adbf6739a9
Merge branch 'Adaptive-Euler-one-loop' into RKCK45-one-loop
2020-03-26 10:56:51 +01:00
fe8ee825d6
bugfix: F, not Fp for explicit nonlocal
2020-03-26 10:52:10 +01:00
8725c2ceba
Merge remote-tracking branch 'origin/development' into Adaptive-Euler-one-loop
2020-03-26 10:49:34 +01:00
a4674a6461
whitespace adjustments
2020-03-25 20:55:37 +01:00
3ce8902245
no need to split up
2020-03-25 20:47:10 +01:00
aa2d440107
finally one loop left
2020-03-25 20:37:37 +01:00
599de26dad
further integration
2020-03-25 19:44:51 +01:00
939d2af1d9
bugfix: explicit nonlocal requires F, not Fe
2020-03-25 18:58:58 +01:00
9900bf0be7
Merge branch 'FPI-one-loop' into RKCK45-one-loop
2020-03-25 18:56:51 +01:00
2d5d48288c
bugfix, related making nonlocal state update explicit
2020-03-25 13:20:14 +01:00
ddef241dcc
Merge branch 'FPI-one-loop' into Euler-one-loop
2020-03-25 13:19:42 +01:00
a6f94efa51
bugfix, introduced when making nonlocal explicit
2020-03-25 13:17:48 +01:00
86abba477a
use same loop
2020-03-25 11:27:03 +01:00
b54a109d99
do plasticState and sourceState at once
2020-03-25 11:20:39 +01:00
652846cdc9
no need for extra loop
2020-03-25 11:16:26 +01:00
67f64cf7e1
correct iteration skipping
2020-03-25 11:13:46 +01:00
32b9a5ab15
all in one loop
2020-03-25 11:07:47 +01:00
9a188784e2
no need for an extra loop
2020-03-25 10:54:17 +01:00
02774a89d9
2 space indentation
2020-03-25 10:33:41 +01:00
52ca7cc43c
only one stateJump per integration
2020-03-25 10:27:55 +01:00
2dd3ccdad1
no need to transpose
2020-03-25 10:22:21 +01:00
8f6dc054d0
move up the nonlocal skip
2020-03-24 22:25:40 +01:00
e7c585a02e
loop order that allows more memory efficient code
2020-03-24 20:47:41 +01:00
66aa20ad39
extra check for nonlocal needed
2020-03-24 20:35:27 +01:00
0740c9f339
dot state in loop
2020-03-24 20:33:26 +01:00
44e24a9c4f
merging into one loop
2020-03-24 20:25:29 +01:00
eb6fe8a3a2
merge into one loop
2020-03-24 16:14:14 +01:00
369ea31a4b
name unification for simple copy and paste
2020-03-24 15:02:55 +01:00
93abf2b6d8
this SEEMS to be the way how it was planned
2020-03-24 13:06:29 +01:00
424fcabb90
move into the same loop
2020-03-24 12:54:31 +01:00
871241e31b
first round of loop removal
2020-03-24 12:43:41 +01:00
b88ffb8d4f
converged(g,i,e) matters only for FPI
2020-03-24 12:30:43 +01:00
d16af3bfb3
one loop for Euler integrator
2020-03-24 12:19:36 +01:00
d45f1f1c6f
further integration into one loop
2020-03-24 11:37:00 +01:00
0e5f0a3068
no need for two loops
2020-03-24 11:27:53 +01:00
106cc1de92
complete iteration of each materialpoint step by step
2020-03-24 10:45:38 +01:00
cde558e736
convergence check on the fly
2020-03-24 10:34:15 +01:00
4c485f1af0
state jump in same loop
2020-03-24 10:18:17 +01:00
3a6269e802
wrong omp pragmas
2020-03-24 10:06:08 +01:00
8d6dcd779c
dotstate/residual/state calculation in the same loop
2020-03-24 06:34:42 +01:00
830ea61739
better readable
2020-03-24 06:27:43 +01:00
9c134e68fc
in-line dot state
2020-03-23 08:15:33 +01:00
e5743c7d27
same logic
2020-03-23 06:53:23 +01:00
027818c942
run in one loop
2020-03-23 00:32:41 +01:00
0b321bd9d4
bugfix for nonlocal
...
logic is hard to understand ...
2020-03-23 00:18:57 +01:00
86db8f8ca4
avoid flush
2020-03-23 00:15:00 +01:00
e7d61e49fe
can be done in the same loop
2020-03-22 23:16:00 +01:00
91e728d065
polishing
2020-03-19 11:30:36 +01:00
9ed48f7e5f
getting rid of totalNslip in nonlocal
2020-03-16 10:09:58 +01:00
66302fa6da
rotational part is always of 3x3 tensor
2020-03-15 14:21:11 +01:00
ae49e6710d
documenting for doxygen
2020-03-15 12:54:35 +01:00
8cc16da53f
atol is part of the state structure, no need for suffix
2020-03-15 09:51:40 +01:00
eb08f9f0b2
polishing
2020-03-14 19:40:05 +01:00
3713810592
simplified allocation for same shape
2020-02-29 13:53:04 +01:00
bcccf06450
Merge remote-tracking branch 'origin/development' into less-public-variables
2020-02-29 13:18:02 +01:00
82dee9db0e
matching names
2020-02-25 17:53:15 +01:00
48604292e2
Merge remote-tracking branch 'origin/development' into MiscImprovements
2020-02-25 17:33:39 +01:00
2c952c3410
these variables can be at leat write protected
2020-02-25 09:50:21 +01:00
7d7eff0d94
crystallite should be responsible of crystallite variables
2020-02-25 09:42:07 +01:00
2cf89b1fdd
bugfix: output dataset was too large in case of Ngrains>1
...
material_phaseAt is defined per constituent/grain/component => no need
to multiply with homogenization_maxNgrains
2020-02-22 11:28:03 +01:00
a2e710c89c
alinged for better readability
2020-02-21 08:45:11 +01:00
c9c78aa90d
do not store invFp and invFi for all points
...
requires to explicitly calculate inverse of Fp and Fi for the tangent
calculation. Hence, classical tradeoff between memory consumption and
runtime.
2020-02-21 08:41:08 +01:00
d108d76a61
was never read
...
probably it was optimized away by the compiler, still confusing
2020-02-20 15:08:32 +01:00
a8e2ee0a86
[skip ci] is read only for other modules
2020-02-14 21:31:03 +01:00
6adb116712
[skip ci] whitespace adjustments
2020-02-14 06:24:17 +01:00
4f7bbb323e
not needed
2020-02-13 18:43:20 +01:00
64e86666c6
also set initial det(Fp)=0
...
Marc element lib test failed otherwise for type 117
2020-02-13 17:10:27 +01:00
0f70a19266
Fp matters, not Fp^-1
...
mathematically absolutely equivalent, but numerically not. Sometikes makes a
huge difference in convergence behavior, even though abs(det(Fp)-1) is
in the order of 1e-15
2020-02-13 14:48:32 +01:00
fc82ec9b1b
tiny deviations from det(Fp) = 1 cause convergence problems
2020-02-13 09:18:49 +01:00
b861ad11c1
tiny differences between orientation conversions result in more
...
iterations
Seen so far only for plasticityDetectChanges test (confirmed for dislotwin and
phenopowerlaw).
The max difference between the entries of the rotation matrix from
orientation0%asMatrix() and eu2om(Eulers)/eu2om(orientation0%asEulers) is
1e-15. This is the ratio of km/atom radius! Still, the number of
iterations is consistently higher.
Results are the same. I believe this is a strange coincidence where one
particular orientation causes problems. The current version recovers almost the 'good' behavior
of math_EulerToR(Eulers)
2020-02-12 20:49:02 +01:00
8770613e9c
better readable
2020-02-12 06:26:22 +01:00
ab475b7c6b
need 'error return'
...
revert from change in 5b72110d
2020-02-12 06:02:37 +01:00
Martin Diehl