Martin Diehl
ffcb7a0a24
aim of the calculation clearer
...
don't know if brackets are honored here
2019-09-23 09:31:14 -07:00
Martin Diehl
da041f43bc
floating point arithmetic refinement
...
Intel 18.4 could not detect 0.0
2019-09-23 09:29:38 -07:00
Martin Diehl
96811f77ca
not used
2019-09-23 06:50:25 -07:00
Martin Diehl
0620eae797
less code
2019-09-23 06:48:39 -07:00
Martin Diehl
611cc9bcae
preprocessor does not allow spaces
2019-09-23 06:48:20 -07:00
Martin Diehl
48e0a7ad97
no floating point comparison for equality
2019-09-23 06:35:02 -07:00
Martin Diehl
b8204015f5
wrong type (math_range returns integers)
2019-09-23 06:29:32 -07:00
Martin Diehl
02fef16b38
used only by MARC/Abaqus
2019-09-23 06:28:10 -07:00
Martin Diehl
78344c01ac
avoid floating point comparison
2019-09-23 05:37:36 -07:00
Martin Diehl
da15ff85f3
tests now operational
2019-09-23 00:11:45 -07:00
Martin Diehl
285dab4451
self consistent solutions
2019-09-22 23:53:56 -07:00
Martin Diehl
bbb2441cc2
preparing more tests
2019-09-22 22:08:19 -07:00
Martin Diehl
0655ef2c90
small precision adjustments required
...
one in a Mio might have degenerated precision...
Also now finally testing core functionality
2019-09-22 21:58:18 -07:00
Martin Diehl
c1398e5fa4
more robust agains overflow + tests
2019-09-22 21:41:00 -07:00
Martin Diehl
dc56556075
trying to resolve MPI restart issues
2019-09-22 19:19:16 -07:00
Martin Diehl
40099fd581
typo
2019-09-22 17:46:37 -07:00
Martin Diehl
9fd65b0617
need to relax tolerance for case 3
2019-09-22 16:57:16 -07:00
Martin Diehl
d1566aa5bf
GNU 8 has issues with IEEE_infinite
...
errors occur when running aggressively optimized code. GNU 9.1 seems to
be ok
2019-09-22 16:05:41 -07:00
Martin Diehl
c0210659b6
statUnit was not opened, stat file was written to STDERR
2019-09-22 15:56:37 -07:00
Martin Diehl
e3de324fee
FEM tolerance was to strict
...
not clear why it pops up now, but there is no real significance beyond
these numbers anyway
2019-09-22 15:35:21 -07:00
Martin Diehl
26bdcd64b9
this might resolve the occasional MPI restart test failuers
2019-09-22 15:34:05 -07:00
Martin Diehl
4ff292ba70
IEEE infinite causes problems with older GNU
2019-09-22 14:46:27 -07:00
Martin Diehl
171a93ca30
also small negative values cause problems
...
Logic differs from the original formulation now
2019-09-22 14:17:17 -07:00
Martin Diehl
6b5b0fae22
mixed up real and aimag part in quaternion
...
+ some tests
2019-09-22 12:10:39 -07:00
Martin Diehl
25c9bb1cd7
testing some corner cases
2019-09-22 11:29:34 -07:00
Martin Diehl
62174a8978
testing all conversions
...
using norm(Delta q) ensures that the quality of the conversion is
tested, not the numerical properties of any orientation representation.
2019-09-22 11:15:27 -07:00
Martin Diehl
1e9d600068
polishing
2019-09-22 08:57:54 -07:00
Martin Diehl
6e5b9ca665
unused variables
2019-09-22 08:54:59 -07:00
Martin Diehl
77011a5dba
this fix restores the previous behavior of dislotwin
...
very strange behavior in plasticity_detect changes for dislotwin. This
fix restores the old behavior with respect to number of iterations.
Using orientation%fromMatrix() results in much more iterations and an
increase in runtime for the test from approx 6 min to 40 min. Results
still match in the end. Also, crystallite_Fp0 differs by approx 1e-15
only between the two methods.
I assume that something is wrong with either dislotwin or the
state/stress integration
2019-09-22 07:46:30 -07:00
Martin Diehl
de632ec85c
better names
2019-09-22 07:22:24 -07:00
Martin Diehl
372536d57e
unit test for rotation class
2019-09-22 06:53:03 -07:00
Martin Diehl
5fa1ecb170
determinant calculations give slightly different results
2019-09-22 06:41:55 -07:00
Martin Diehl
30afaf2a95
bugfic: Blocksize parameter was too small
2019-09-21 23:45:54 -07:00
Martin Diehl
51bd67fa29
findloc is in gfortran 9.1
2019-09-21 23:03:53 -07:00
Martin Diehl
25d91c79af
larger block size seems favorable
2019-09-21 23:03:27 -07:00
Martin Diehl
8b908fb350
some more testing
2019-09-21 21:25:55 -07:00
Martin Diehl
a963f1d2c3
better explicit
2019-09-21 09:16:58 -07:00
Martin Diehl
f5caee8026
issue now clearer
2019-09-21 08:23:56 -07:00
Martin Diehl
7ef4aca170
vectorized
...
there seems to be a conflict in the definition of the projection for
edge and screw.
Nonlocal uses the transpose compared to dislotwin/disloUCLA.
2019-09-21 08:21:39 -07:00
Martin Diehl
8013dec45f
safer
2019-09-21 07:14:45 -07:00
Martin Diehl
9e16f27fe9
proper capitalization, no more synonyms
2019-09-20 19:01:44 -07:00
Martin Diehl
fd4a0eb12e
cleaning:
...
- porosity/vacancy does not exist any more
- texture becomes orientation
- crystallite is getting removed
2019-09-20 19:00:20 -07:00
Martin Diehl
7837a6dc01
not used ...
...
... and equivalentStress33 seemed wrong
... and equivalentStrain33 was confusing
2019-09-20 18:45:54 -07:00
Martin Diehl
a17046cfdf
full merger
2019-09-20 18:44:23 -07:00
Martin Diehl
a16af9a7ad
single source of truth
2019-09-20 18:33:12 -07:00
Martin Diehl
97127f9034
following standard: Output first
2019-09-20 18:28:46 -07:00
Martin Diehl
53283d5c01
using newer interface
2019-09-20 18:20:33 -07:00
Martin Diehl
09c1150e3c
more error checking
2019-09-20 18:16:08 -07:00
Martin Diehl
c1d1c83088
better check input data
2019-09-20 17:49:55 -07:00
Martin Diehl
a8b72c3f53
even better to read
2019-09-20 17:42:55 -07:00