Commit Graph

218 Commits

Author SHA1 Message Date
Martin Diehl 53283d5c01 using newer interface 2019-09-20 18:20:33 -07:00
Martin Diehl 6c0e92d5c1 functionality to rotate 4-tensor 2019-09-20 16:47:44 -07:00
Martin Diehl 42fba28fa1 getting rid of conversions with unclear behavior 2019-09-20 07:37:58 -07:00
Martin Diehl 5d0a67e5dc texture class replaces manual conversions 2019-09-19 18:25:51 -07:00
Martin Diehl b41d6182a5 safe use of optional argument 2019-09-19 14:20:02 -07:00
Martin Diehl f6c616c84b consistent orientation/quaternion handling is in own submodule 2019-09-19 13:59:18 -07:00
Martin Diehl 7fa979f800 removed unused function
were based on/used for deprecated orientation handling
2019-09-19 13:51:37 -07:00
Martin Diehl be0d961954 cleaning 2019-09-19 13:40:03 -07:00
Martin Diehl e30478127d cleaned 2019-06-15 20:42:16 +02:00
Martin Diehl 9dfe71aa06 better readable 2019-06-10 09:42:23 +02:00
Martin Diehl 7ac0013271 more consistent private/public declarations 2019-05-17 05:24:36 +00:00
Martin Diehl 7b620e3ce9 [skip ci] consistent with rest of the module 2019-05-17 06:49:25 +02:00
Martin Diehl 436dae8dd5 cleaning 2019-05-11 11:52:05 +02:00
Martin Diehl 93ded62e32 better readable 2019-05-10 16:38:47 +02:00
Martin Diehl b8c8c4ade6 follow numpy naming convention 2019-05-10 15:19:00 +02:00
Martin Diehl 1cf6cc217d one implicit none is enough 2019-05-10 15:12:40 +02:00
Martin Diehl 811fc5c30c don't put implementation details in function name 2019-05-10 15:03:54 +02:00
Martin Diehl c81a90350c not needed anymore 2019-05-10 15:00:48 +02:00
Martin Diehl dc57f6aba0 better use transpose 2019-05-10 14:53:22 +02:00
Martin Diehl e25b4d6f17 math for quaterion / orientation have own classes 2019-05-10 14:52:01 +02:00
Martin Diehl 23cf134d6c cleaning 2 2019-05-08 22:41:09 +02:00
Martin Diehl 453eb538f7 preparing for PGI compiler 2019-04-03 16:24:15 +00:00
Martin Diehl 3887d5631a handling of PGI compiler was to cumbersome 2019-04-03 13:04:03 +02:00
Martin Diehl 4604e65a42 use matmul instead of hand-written functions
- performance is the same
- leaner code
- matmul works (was buggy a few years ago)
2019-04-03 08:22:04 +02:00
Martin Diehl 66aa9399a7 don'r rely on halton sequence
not really trustworthy as we saw some pattern for higher dimensional
data
2019-03-17 18:06:41 +01:00
Martin Diehl c7f33f4696 no pInt 2019-03-13 06:46:59 +01:00
Martin Diehl 55d55a156b following numpy syntax 2019-03-09 16:19:56 +00:00
Martin Diehl a4cdbab7cf not needed 2019-03-07 22:44:33 +01:00
Martin Diehl 4a636dbd99 not used at all 2019-03-07 22:41:07 +01:00
Martin Diehl 8224797e75 simplified 2019-03-07 22:31:54 +01:00
Martin Diehl f0eeb3d9b4 using naming convention of numpy 2019-03-07 22:22:49 +01:00
Martin Diehl d83847446d Merge branch '32_NewStyleNonlocal-4' into development 2019-02-25 20:06:20 +01:00
Martin Diehl 93ae5cbd07 Merge remote-tracking branch 'origin/development' into 10-consistent-orientation-conversions-3 2019-02-22 23:05:36 +01:00
Martin Diehl cb2d2b02dc re-ordered and cleaned 2019-02-21 00:24:35 +01:00
Martin Diehl 7ad866b90f not needed any more 2019-02-17 12:16:12 +01:00
Martin Diehl 87f3e3f621 more flexible and user friendly 2019-02-11 10:08:34 +01:00
Martin Diehl 09859f1b12 wrong variable rename (was forgotten) 2019-02-09 16:53:05 +01:00
Martin Diehl de26e41684 some first steps to support debugging with the PGI compiler
norm2 and sum for initialization are not supported yet, need fixes
2019-02-09 15:02:07 +00:00
Martin Diehl a260bd2d2b Merge branch 'development' into 10-consistent-orientation-conversions 2019-02-01 07:33:52 +01:00
Martin Diehl 7484849b22 Merge remote-tracking branch 'origin/development' into 42-new-coding-style-for-homogenization-NEW 2019-01-27 08:44:46 +01:00
Martin Diehl 221c587362 using separate functions for stress and tangent
extensively tested in 46-simplification-of-crystallite-f90-NEW3 already
2019-01-18 15:30:50 +01:00
Martin Diehl 54bb6a5eb6 more verbose name
credits to Franz
2019-01-14 16:36:08 +01:00
Martin Diehl adb68ff792 avoid converstions and use of global variables 2019-01-14 12:45:07 +01:00
Martin Diehl f82a422e21 do concurrent causes problems on matesting/Intel 18.0.1 2019-01-14 09:34:14 +01:00
Martin Diehl 01a2fffd3c clearer naming for vector <-> tensor conversion
and a bunch of other improvements
2019-01-14 07:27:18 +01:00
Martin Diehl 619baefe19 some comments
calling a conversion "Mandel" that does not follow the Mandel convention
(at least according to wikipedia) is not really intuitive
2019-01-13 18:54:20 +01:00
Martin Diehl fe28e0d739 better follow the same conventions everywhere 2019-01-13 08:47:01 +01:00
Franz Roters cd8ee4503b Merge branch '22-NewStyle_disloUCLA-2' into 'development'
Resolve "New coding style for plastic_disloUCLA"

See merge request damask/DAMASK!51
2018-12-18 15:02:12 +01:00
Martin Diehl ccdf1e5e8e polishing 2018-12-08 13:49:42 +01:00
Martin Diehl 5768356642 PGI compatible 2018-12-06 01:12:15 +01:00
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