Commit Graph

218 Commits

Author SHA1 Message Date
Martin Diehl 91a2748131 simplifying 2018-12-04 20:33:28 +01:00
Martin Diehl 94695f773e more verbose and works for arbitrary precision 2018-08-28 07:47:05 +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 0a6bde70c5 Merge branch 'development' into 19-NewStylePhenopowerlaw 2018-08-03 06:56:43 +02:00
Martin Diehl 449449b500 does the same as numpy.clip 2018-07-28 01:31:02 +02:00
Martin Diehl d7da70cefb Merge branch '30_parsePhasePartOnce' into 19-NewStylePhenopowerlaw 2018-06-19 19:30:01 +02:00
Martin Diehl 32d481020f Merge branch 'development' into 30_parsePhasePartOnce 2018-06-16 14:49:57 +02:00
Martin Diehl f8ae0ec925 list does tokenizing automatically 2018-06-10 11:07:17 +02:00
Vitesh Shah c903880d19 using on unitialized error caused segfault in some cases 2018-06-04 12:38:28 +02:00
Martin Diehl 804febe7f9 WIP: separate reading in and parsing 2018-05-31 18:31:44 +02:00
Martin Diehl 72b69959de Merge branch 'development' into 19-NewStylePhenopowerlaw 2018-05-31 17:02:45 +02:00
Martin Diehl f0309ac4dd Bessel functions not used 2018-05-24 07:51:16 +02:00
Martin Diehl 6b27460a3b Merge branch 'development' into 12-fixOrientationSampling 2018-05-05 19:03:51 +02:00
Martin Diehl 235b2fe3c8 ntwin = [0] cause problems
math_expand now works for sum(how) == 0
2018-05-02 16:44:27 +02:00
Martin Diehl 81cab02d7a simplified reading in and initialization
init of absolute tolerances and state0 done on the fly
reading in array type parameters in the param structure
:
2018-04-25 19:43:03 +02:00
Martin Diehl 8509b13f16 avoid numerical noise 2018-03-13 09:21:41 +01:00
Martin Diehl 302a5017d4 Merge branch 'development' into 12-fixOrientationSampling 2018-03-13 04:12:03 +01:00
Martin Diehl 2ceb3aeea0 not all variables were consistently renamed 2018-03-13 03:56:02 +01:00
Martin Diehl 44e5644e78 fixed random Gaussian sampling
sampling needs to be performed from unfiform misorientation, NOT uniformly distributed rotations
for Fiber, compute uniform tilt of Fiber axis
2018-03-13 03:36:24 +01:00
Zhuowen Zhao 361cf680a5 merged more recent development 2018-03-09 18:19:41 -05:00
Zhuowen Zhao f26fd1d1dc Fixed a physics issue: sense change with respect to difference between resolved stress and backstress 2018-03-09 18:01:09 -05:00
Martin Diehl 0ccce7facc fixed missing angle initialization, simplified and commented 2018-03-08 17:33:54 +01:00
Martin Diehl e716441787 indexing now follows description 2018-02-26 20:18:38 +01:00
Martin Diehl cd4115bada now using the correct rotation formula 2018-02-26 11:14:02 +01:00
Martin Diehl e1c4c7e6e8 wrong indices in matrix rotation
fixing this resolves issue 14: "Poor performance after load case rotation"
now inline with http://www.continuummechanics.org/rotationmatrix.html
C'_mnop = R_mi R_nj R_ok R_pl C_ijkl
2018-02-25 21:01:23 +01:00
Martin Diehl 4f4fa5daf8 simplified 2018-02-22 14:54:49 +01:00
Martin Diehl e51de7ffd8 explicitly select halton bases in call 2018-02-22 14:16:36 +01:00
Martin Diehl c6c66bb653 using proper Gauss sampling also for the fiber components
function most probably still contains a bu
2018-02-22 00:42:06 +01:00
Martin Diehl 5c908e44ec correct scaling for FWHM
FWHM was wrong by a factor of approx 2.
Analytic expression can be found in https://en.wikipedia.org/wiki/Gaussian_function
2018-02-22 00:10:38 +01:00
Martin Diehl 9173d12d14 correct algorithm for sampling of uniform orientations and fix for Halton series
Halton series gives strange results for large prime numbers, now always starting with 2 for first dimension, 3 for second etc.
Consecutive Halton numbers for rejection sampling seem to cause problems (i.e. introduce patterns).
Algorithm for uniformly distributed orientations with FWHM specified is taken from https://math.stackexchange.com/questions/131336.
WIP: Gauss filtering is currently not implemented!
2018-02-21 20:32:52 +01:00
Martin Diehl ae27660e86 simplified halto procedure (still needs testing) 2018-02-21 18:47:39 +01:00
Martin Diehl 98df2d1427 better description and names 2018-02-21 16:49:29 +01:00
Martin Diehl a52f54a9a0 fixed typo in prime number list and extended to 1600 values
source https://people.sc.fsu.edu/~jburkardt/f_src/halton/halton.f90
2018-02-21 14:37:02 +01:00
Martin Diehl 6b1c1af609 no early return required here 2018-02-21 08:55:33 +01:00
Martin Diehl 9b1e72e7bc more sensible cutoff for scatter 2018-02-15 11:26:19 +01:00
Martin Diehl 230ffe8f35 Bessel functions for Gaussian scatter 2018-02-14 13:03:50 +01:00
Martin Diehl fe07a34c69 source code should not be executable 2018-02-12 12:42:18 +01:00
Martin Diehl b6abbfca9d early return not needed here 2018-02-12 11:45:42 +01:00
Martin Diehl 2254746177 Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development 2018-02-12 11:23:05 +01:00
Zhuowen Zhao 7f487bb77b merge development into kinematic hardening branch 2018-02-07 11:30:57 -05:00
Fengbo Han 9249e7db2b Merge branch 'development' into HMS 2018-02-07 11:44:17 +01:00
Martin Diehl b455fcd0a3 acos has the same limits/stability as indirect asin calculation, enforing limits directly instead of scaling 2018-02-03 09:21:19 +01:00
Martin Diehl 9908734710 wrong random number, simplified 2018-02-03 09:04:49 +01:00
Martin Diehl d80e15bd76 compiler_options and compiler_version supported by new Intel compiler 2018-02-02 12:36:09 +01:00
Zhuowen Zhao 7c755a0241 Merge remote-tracking branch 'origin/spectralSolver-cutbackfix' into 3-adding-plastic-constitutive-law-with-kinematic-hardening 2018-01-18 13:59:42 -05:00
Philip Eisenlohr e4700cda25 changed fixed_seed to random_seed for clarity 2017-11-07 14:56:28 -05:00
Zhuowen Zhao ab11830e71 Merge branch '6-debug-replaces-openmp-preprocessor-switch' into 3-adding-plastic-constitutive-law-with-kinematic-hardening 2017-11-06 18:20:51 -05:00
Fengbo Han d599aec37f Merge branch 'development' into HMS 2017-10-24 11:24:54 +02:00
Fengbo Han 82758bd90f added yield criterion of facet potential 2017-10-24 11:15:34 +02:00
Martin Diehl eda8c24b7d Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development 2017-10-07 17:24:59 +02:00
Martin Diehl e32b9d9ca8 for comparison with de-facto stardard rotation definitions 2017-10-07 13:18:42 +02:00
Martin Diehl e9440ccd6f suppress warning when compiling with gfortran 2017-10-05 16:35:34 +02:00
Zhuowen Zhao 7d6d822ee0 Merge branch '5-flexible-dot-delta-state-relation' into 3-adding-plastic-constitutive-law-with-kinematic-hardening 2017-09-29 18:46:35 -04:00
Zhuowen Zhao 80bed8b8aa improved math_expand algorithm 2017-09-29 18:32:52 -04:00
Martin Diehl d8d42c32e7 cleaned + shortened, fixed handling of optional order parameter in math_exp33 2017-09-19 23:39:19 +02:00
Zhuowen Zhao ae868d3ada only one error for math_check; new "math_expand('what' by 'how')" 2017-09-14 15:25:22 -04:00
Martin Diehl f0d8d7c50a bug in eigenvalue calculation, caused wrong eigenvectors for sym33 matrices. Credits to Duancheng Ma 2017-09-06 16:20:24 +02:00
Martin Diehl 7e0a9b75f9 simplified prime fuction 2017-08-12 06:35:44 +02:00
Martin Diehl cf6894442b moved specific functions into the scope of the calling functions 2017-08-12 06:03:40 +02:00
Martin Diehl 9823f5f495 one time used variables with name a and x to not help, assignement of i and j via loop is safe 2017-08-12 04:51:10 +02:00
Fengbo Han 5cedba0721 implemented yield stop criteria 2017-07-27 16:21:02 +02:00
Martin Diehl b35ff67f99 using IEEE_is_NAN and IEEE_quiet_NaN instead of hand-written solution, will not work for gfortran < 5 2017-05-04 00:32:44 +02:00
Martin Diehl a5f8c5ba32 Merge branch 'development' into cmake 2017-04-15 16:17:50 +02:00
Martin Diehl a248a24250 consistent with other routines 2017-02-08 20:01:42 +01:00
Martin Diehl d86074829c Merge branch 'development' into cmake 2017-02-03 19:28:23 +01:00
zhangc43 3a561c9a8d remove obsolete j2 module 2016-04-26 14:24:43 -04:00
Chen Zhang 672c8d0606 Merge remote branch 'origin/development' into cmake 2016-04-18 09:11:29 -04:00
Chen Zhang e33594cb44 change code folder to src 2016-03-04 14:43:29 -05:00