Commit Graph

14100 Commits

Author SHA1 Message Date
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
Martin Diehl 7fe3e3b406 easier to read 2019-09-20 17:38:57 -07:00
Martin Diehl 44125f3f09 convenience function for symmetric 4th order tensors
unfortunately, array size does not qualify to distinguish generics
2019-09-20 17:36:37 -07:00
Martin Diehl ad83c8541d same names as in python 2019-09-20 17:18:09 -07:00
Martin Diehl 7e6265b128 make use of consistency checks 2019-09-20 16:52:55 -07:00
Martin Diehl ef4204da14 Merge branch 'MiscImprovements' of magit1.mpie.de:/damask/DAMASK into MiscImprovements 2019-09-20 16:48:50 -07:00
Martin Diehl 6c0e92d5c1 functionality to rotate 4-tensor 2019-09-20 16:47:44 -07:00
Martin Diehl 4dfb1f6450 next try for MARC 2019 2019-09-21 01:07:04 +02:00
Martin Diehl 93c8aa8f75 bugfix: used uninitialized vector 2019-09-21 00:08:14 +02:00