Martin Diehl
42bc1b468a
cleaning+adopting to petsc-3.10
2018-09-22 11:19:08 +02:00
Martin Diehl
a95e2b8083
Merge branch 'PETSc3.10' into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-09-22 10:36:31 +02:00
Martin Diehl
f590851a78
Merge branch '38-introduce-rudimentary-PETSc-based-FEM-solver' of magit1.mpie.de:damask/DAMASK into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-09-22 10:33:20 +02:00
Martin Diehl
80914470cc
most probably pInt and PetscInt is the same
2018-09-22 10:32:26 +02:00
Jaeyong Jung
4ddb8d843c
modified DMGetDefaultGlobalSection to DMGetGlobalSection
2018-09-21 14:08:44 +02:00
Jaeyong Jung
9150844c96
Merge branch '38-introduce-rudimentary-PETSc-based-FEM-solver' of magit1.mpie.de:damask/DAMASK into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-09-21 10:13:27 +02:00
Martin Diehl
180105d3ce
external statements will not work if interfaces exist
...
we're currently providing a patch for PETSc with interfaces for all
functions that we use
2018-09-21 08:25:35 +02:00
Martin Diehl
ff241ae4a6
only use required includes
2018-09-21 08:19:36 +02:00
Jaeyong Jung
abe5c1d825
Merge branch '38-introduce-rudimentary-PETSc-based-FEM-solver' of magit1.mpie.de:damask/DAMASK into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-09-20 09:57:30 +02:00
Martin Diehl
b775d5084b
more (temporarly) cleaning
2018-09-20 08:26:59 +02:00
Martin Diehl
1545a53ea9
ID do not exist at the moment
2018-09-20 08:00:39 +02:00
Martin Diehl
cc262ae198
Merge branch 'development' into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-09-20 07:29:19 +02:00
Martin Diehl
3e7b80a3ef
debug only available if compiled in debug mode
2018-09-20 07:27:12 +02:00
Martin Diehl
fc016bbc6e
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-09-20 07:16:03 +02:00
Martin Diehl
310ea62964
only print out the essential information
2018-09-20 06:58:31 +02:00
Martin Diehl
fcff6b908a
can be easily computed during post processing
2018-09-20 06:40:23 +02:00
Martin Diehl
901355d2ae
don't use unnecessarily long names
2018-09-20 06:27:53 +02:00
Martin Diehl
df0464c31b
use (import) debug variables only when needed
2018-09-20 06:24:03 +02:00
Martin Diehl
6aa4dd842a
define debug variables only if needed
2018-09-20 06:09:02 +02:00
Martin Diehl
1623a33b48
cleaning (mainly OMP FLUSh)
2018-09-19 21:46:26 +02:00
Martin Diehl
191ad9df09
not used
2018-09-19 20:40:35 +02:00
Martin Diehl
1bf1e13b46
some debug statements
2018-09-19 20:23:28 +02:00
Martin Diehl
e7e959af47
using explicit interface from PETScDT
2018-09-19 20:06:35 +02:00
Martin Diehl
11d4c28d88
flushes not needed + further cleaning
2018-09-19 19:45:57 +02:00
Martin Diehl
b84476e681
cleaning and debugging
2018-09-19 18:33:37 +02:00
Martin Diehl
c313dc1675
only read access
2018-09-19 17:04:12 +02:00
Martin Diehl
0bf64645a1
should be done by the plasticity laws (for the moment)
2018-09-19 14:22:35 +02:00
Martin Diehl
a8fb7d7ade
not needed
...
but I'm under the impression that the compiler removes such things anyway
2018-09-19 14:21:10 +02:00
Martin Diehl
20f0bee459
fallback dPdF not needed
...
save a lot of memory
2018-09-19 06:19:40 +02:00
Martin Diehl
18db5f5652
not needed
2018-09-18 05:20:08 +02:00
Martin Diehl
82cdc551a5
MPI Communicator only needed for PETSc 3.10.
...
thx to Jaeyong for figuring this out
2018-09-18 05:17:48 +02:00
Martin Diehl
245676728d
not needed
2018-09-17 23:28:00 +02:00
Jaeyong Jung
3beab556a2
runtime polar decompoistion error comes up
2018-09-17 15:41:58 +02:00
Martin Diehl
e6fa3f3d35
correct stress in postResults
2018-09-16 22:57:50 +02:00
Martin Diehl
cea2fba063
don't repeat code that can be handled centrally
2018-09-16 22:38:57 +02:00
Martin Diehl
9570fb894a
correct names and no superflous conversions anymore
2018-09-16 22:31:55 +02:00
Martin Diehl
5f06a35900
Merge branch 'development' into 19-NewStylePhenopowerlaw
2018-09-16 22:03:05 +02:00
Martin Diehl
b54742cfa7
[skip sc] [skip ci] PETSc 3.10 is out
2018-09-15 20:09:27 +02:00
Martin Diehl
075d2e51b2
debugging, still crashing
2018-09-15 16:19:41 +02:00
Martin Diehl
3b1817a954
cleaning, only having mech at the moment
2018-09-15 15:14:48 +02:00
Martin Diehl
934ca56811
clean exit if required files not found
2018-09-15 14:54:56 +02:00
Martin Diehl
a92b945e3f
does not crash anymore
...
xx is nc, the number of components
needs to be 3 in the current case
2018-09-15 11:40:07 +02:00
Martin Diehl
5e33900664
Merge branch 'development' into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-09-15 08:38:30 +02:00
Martin Diehl
25a38ad438
initialization broken
...
- need to consider case of no twin/slip active
- state0 needs to be initialized
2018-09-15 05:45:10 +02:00
Martin Diehl
9f16cefd9f
renaming in accordance with the DAMASK paper
2018-09-14 11:17:14 +02:00
Martin Diehl
af32b3d85b
reordering for easier overview
2018-09-14 10:51:44 +02:00
Martin Diehl
3dd47eade9
just renaming, better readable...
2018-09-14 10:18:44 +02:00
Martin Diehl
07b0ddf711
polishing
2018-09-14 10:08:36 +02:00
Martin Diehl
18858301d5
using consistently i as first running index and j as second
2018-09-14 09:59:04 +02:00
Martin Diehl
55511782e8
incomplete transition/renaming
...
caused segmentation fault
2018-09-14 08:19:39 +02:00
Martin Diehl
f98243e4ac
Merge branch 'development' into 43-wrong-stress-in-dotstate-and-deltastate
2018-09-14 05:39:15 +02:00
Martin Diehl
e2a66f6ddb
should not differ from development
2018-09-14 05:33:39 +02:00
Martin Diehl
b165e48d9f
non-Schmid contribution were ignored in test
...
reverted, just to ensure that the results have not changed: Never
polishing and change physics at the same time
2018-09-14 05:21:05 +02:00
Martin Diehl
cc20b044a6
forgot to allocate array now used internally only
2018-09-13 09:27:54 +02:00
Martin Diehl
ae8efe2755
Merge branch 'development' into 44-column-major-access-to-interaction-matrices
2018-09-13 05:55:59 +02:00
Martin Diehl
8ecb019566
real must be pReal for MSC.Marc
2018-09-13 05:51:58 +02:00
Martin Diehl
4266872965
array out of bounds was possible for Ntwin =0
...
in postResults(c+1:c+prm%totalNtwin) c+1 will be out of bounds if
c is already size(postResults)
2018-09-13 05:41:44 +02:00
Martin Diehl
3068caa9a3
[skip sc] resolved stress not needed
...
using kinetics_xxx as in disloUCLA
compiles on gfortran but pre-receive hook with intel compiler (MSC.Marc)
fails
2018-09-12 21:05:40 +02:00
Martin Diehl
4e68d049a8
WIP: calculating rotated stiffness matrices for transformation
2018-09-12 20:37:55 +02:00
Martin Diehl
b95174a8b7
reducing global variables
2018-09-12 16:16:31 +02:00
Martin Diehl
192bb6453d
cleaning
2018-09-12 15:57:54 +02:00
Martin Diehl
0a621ece8b
WIP: polishing and implementing remaining functions
2018-09-12 14:29:22 +02:00
Martin Diehl
f1f8922ab8
more meaningful error messages
2018-09-12 14:03:45 +02:00
Martin Diehl
edebe4d1ed
vectorization error
2018-09-12 12:26:33 +02:00
Martin Diehl
c9b5b3fb7b
should be always set, even if no twinning is enabled
...
fixes NaN in dotState
2018-09-12 12:08:03 +02:00
Martin Diehl
253f318e55
preventing NaN in dotState
2018-09-12 11:47:35 +02:00
Martin Diehl
2337dde525
cleaning
2018-09-12 11:26:59 +02:00
Martin Diehl
8a406150f8
polishing
2018-09-12 10:06:04 +02:00
Martin Diehl
3352611267
cleaning
2018-09-12 09:59:09 +02:00
Martin Diehl
01196ca36b
Merge branch 'development' into 19-NewStylePhenopowerlaw
2018-09-12 09:50:24 +02:00
Martin Diehl
c9208315f5
avoid repeated calculations
...
does not save so much here, but avoids having inconsistent calculation
(e.g. nonSchmid effects) and serves as a template for more complex
models
2018-09-12 09:23:11 +02:00
Martin Diehl
ce6e6679d5
causes array access out of bounds, needs further checks
2018-09-11 12:20:05 +02:00
Martin Diehl
3d49c70dbc
commenting
2018-09-08 20:15:17 +02:00
Philip Eisenlohr
9be2c084e4
Merge branch '39-simplify-obscure-numerics-integration-mode' into 'development'
...
Resolve "simplify obscure numerics integration mode"
Closes #39
See merge request damask/DAMASK!36
2018-09-07 19:34:35 +02:00
Philip Eisenlohr
7e330f1249
Merge branch '37-include-hdf5-routines' into 'development'
...
Resolve "Include HDF5 routines"
Closes #37
See merge request damask/DAMASK!37
2018-09-07 19:18:07 +02:00
Satyapriya Gupta
a16454f1f9
fixed missing propagation of dependent state variables
2018-09-06 16:48:28 -04:00
Martin Diehl
e70efd2608
inform about debug versions indenpendently of the solver
2018-09-06 16:49:39 +02:00
Martin Diehl
8324962123
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-09-06 16:26:46 +02:00
Martin Diehl
cf65aae92a
correct names
2018-09-01 06:49:03 +02:00
Philip Eisenlohr
6a64637ac3
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-08-31 11:23:25 -04:00
Philip Eisenlohr
3e81110830
tell user whether DEBUG version or not
2018-08-31 11:22:59 -04:00
Martin Diehl
3b88096cd7
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-08-31 11:09:52 +02:00
Martin Diehl
5330fb31ad
more verbose error message
2018-08-31 10:19:37 +02:00
Martin Diehl
f028e05298
guessing was not correctly set
2018-08-31 10:14:33 +02:00
Martin Diehl
3cb279b083
one more loop not needed
2018-08-31 10:03:15 +02:00
Martin Diehl
017563e061
one more loop not needed
2018-08-31 09:14:16 +02:00
Martin Diehl
c0ac05aa88
no need to use two loops
2018-08-31 09:12:50 +02:00
Martin Diehl
732022d4f0
functions from IO does not add any value here
...
only prevents the automated selection of a new unit
2018-08-31 08:58:13 +02:00
Martin Diehl
d18d74ae31
mechanics solver type is set only once
...
as all mech routines have the same interface, this leaner syntax can be
used
2018-08-31 08:20:23 +02:00
Martin Diehl
555ca8e5f1
debug statements not needed anymore
2018-08-30 17:21:22 +02:00
Martin Diehl
62413e8187
now compiles again
2018-08-30 17:18:56 +02:00
Martin Diehl
5f97ae9080
not used at the moment
2018-08-30 16:39:18 +02:00
Martin Diehl
cbecb71000
[skip sc] using explicit interface helps to get correct type
...
still some issues with different ranks and datatypes
2018-08-30 13:10:18 +02:00
Martin Diehl
b930c68887
PETSc 3.9 style (no SIGSEGVs anymore)
2018-08-30 12:38:44 +02:00
Martin Diehl
8b6bef159e
a little bit more verbose
2018-08-30 12:38:29 +02:00
Martin Diehl
3254ae79a5
collecting questions to Pratheek
2018-08-30 12:37:47 +02:00
Martin Diehl
d180af494e
adopted to new PETSc fortran interface
2018-08-30 11:55:13 +02:00
Martin Diehl
29e55d20fb
message better to understand and giving error instead of SIGSEGV
2018-08-30 09:42:45 +02:00
Martin Diehl
b96e01c128
Merge branch 'development' into lattice-builds-interaction
2018-08-30 05:41:40 +02:00
Martin Diehl
105c4077ad
Merge branch 'development' into 19-NewStylePhenopowerlaw
2018-08-30 00:58:54 +02:00
Martin Diehl
b163a8aaa0
Non-schmid activated again
...
internally, no need for long name plastic_phenopowerlaw_postResults
2018-08-30 00:57:15 +02:00
Martin Diehl
b884349e7b
only renaming
...
3333 not needed for dX_dY if X and Y are 3x3 tensors
PK2 stress is S not T according to the DAMASK paper
2018-08-29 13:16:37 +02:00
Martin Diehl
6d28883db5
starting to include non-Schmid terms (again)
2018-08-29 12:07:51 +02:00
Martin Diehl
922273f230
does not make sense to store and use the 6-vector version of the Schmid
...
matrix
2018-08-29 12:03:31 +02:00
Martin Diehl
baeb449e07
WIP: debugging
2018-08-29 11:44:39 +02:00
Martin Diehl
f2acc14882
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-08-29 08:27:15 +02:00
Martin Diehl
680c9e11d4
segmentation fault in cause of empty list
2018-08-29 08:26:28 +02:00
Martin Diehl
b5e4ad247d
LpAndItsTangent naming convention following DAMASK paper
2018-08-28 15:40:24 +02:00
Martin Diehl
fc3ce54667
return more than one tangent
2018-08-28 15:07:39 +02:00
Martin Diehl
15e2d4a7cd
Tstar is S
...
according to the DAMASK paper
2018-08-28 15:02:58 +02:00
Martin Diehl
4c14f988a3
rename according to paper
2018-08-28 14:58:17 +02:00
Martin Diehl
06e7156351
did not work for values <0
2018-08-28 12:41:21 +02:00
Martin Diehl
42f8b0a063
labels of slip and twin systems for more self-explanatory output
2018-08-28 12:27:56 +02:00
Martin Diehl
94695f773e
more verbose and works for arbitrary precision
2018-08-28 07:47:05 +02:00
Franz Roters
38fd517c44
Merge remote-tracking branch 'remotes/origin/36-faster-file-handling-for-material-config-use-stream-access-instead-of-line-wise-reading' into development
2018-08-27 15:13:16 +02:00
Martin Diehl
38a855243d
Merge branch 'development' into 43-wrong-stress-in-dotstate-and-deltastate
2018-08-27 06:31:45 +02:00
Martin Diehl
dd81fe806d
naming was not completely corrected
2018-08-27 06:30:55 +02:00
Martin Diehl
b11666effa
Merge branch 'development' into 38-introduce-rudimentary-PETSc-based-FEM-solver
2018-08-26 13:59:50 +02:00
Martin Diehl
a5f139b786
unused variables
2018-08-26 11:53:18 +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
0041d21777
output was screwed up
2018-08-26 06:20:48 +02:00
Martin Diehl
4f1becb503
cleaning and fixing bugs (wrong indices)
2018-08-25 21:27:14 +02:00
Martin Diehl
f458de82fa
simplified
2018-08-25 21:13:32 +02:00
Martin Diehl
3ff7c9c0eb
extendend to Schmid_twin
2018-08-25 20:55:15 +02:00
Martin Diehl
fbaac484ea
extended to postResults
2018-08-25 20:32:01 +02:00
Martin Diehl
780699740d
internal (non)-Schmid definitions also in dotState
2018-08-25 20:19:02 +02:00
Martin Diehl
ba9330c8dd
early rename
2018-08-25 20:08:46 +02:00
Martin Diehl
d51b00c2d3
no need to loop over families and systems independently. no need for 3333
2018-08-25 20:06:19 +02:00
Martin Diehl
d055ef665a
using only internal schmid and non schmid matrices
2018-08-25 20:01:39 +02:00
Martin Diehl
cc5d04ff2a
newly introduced parameters for non-schmid were confusing
2018-08-25 19:28:09 +02:00
Martin Diehl
5b3bed3e48
functionality to calculate SchmidMatrices
2018-08-25 18:32:55 +02:00
Martin Diehl
d3db0d816c
ordered by dimension and set private where possible
2018-08-25 17:29:20 +02:00
Martin Diehl
1c76e156d9
new symmetry handling will be handled in a different branch
2018-08-25 17:19:48 +02:00
Martin Diehl
bde242902b
Merge branch 'development' into 43-wrong-stress-in-dotstate-and-deltastate
2018-08-25 16:02:12 +02:00
Martin Diehl
b8b5bac684
dotState and deltaState parse Mstar instead of Tstar
...
requires to parse in Fi
2018-08-25 15:59:34 +02:00
Martin Diehl
22a232ad08
bug: memory access out of bounds
...
introduced when moderninzing reading in of parameters
2018-08-25 15:50:43 +02:00
Martin Diehl
7e002d1bfb
missing rename
2018-08-25 14:42:44 +02:00
Martin Diehl
e46605f0ef
forall is deprecated
2018-08-25 14:36:21 +02:00
Martin Diehl
17d88184a7
typos fixed
2018-08-25 14:35:49 +02:00
Martin Diehl
35688a6acf
temp33 not needed (compiler should be smart enough)
2018-08-25 14:29:03 +02:00
Martin Diehl
953fff79ac
prepared LpAndItsTangent to remove superflous forward-backward
...
conversion
updating the individual plastic laws will be done in the respective
branches
2018-08-25 14:21:06 +02:00
Martin Diehl
d585deee7e
unified syntax
...
_v for vector representation of tensor conflicts with 3333 suffix for
4-th order tensors.
General idea: Mark symmetric second order tensors in vector notation with
'6' and fourth order tensors in second order matrix notation with '99'.
Append nothing for 'natural' representation, i.e. F and NOT F33, Lp_dTstar and NOT
Lp_dTstar3333
2018-08-25 14:07:23 +02:00
Martin Diehl
025cbddd00
generic function and interfaces for slipSlip, twinTwin, transTrans
2018-08-25 13:08:32 +02:00
Martin Diehl
7ecb7689f1
Intel compiler failed with SIGSEV
...
derived types, pointers, finalize .... altogether seems to bring both Compilers to their limits.
I cannot see what was wrong before, but now it works and might be a little faster
2018-08-23 00:13:57 +02:00
Martin Diehl
271b9ba76b
intersting note ...
2018-08-22 23:28:47 +02:00
Martin Diehl
8c5f3d4e07
only needed once
2018-08-22 22:52:12 +02:00
Martin Diehl
51b4ef319a
HDF5 was not compiled due to missing dependency
2018-08-22 20:12:16 +02:00
Martin Diehl
dc596e6789
zero termination does not work
2018-08-22 18:09:17 +02:00
Martin Diehl
1a943df97e
small flaws
2018-08-22 17:52:00 +02:00
Martin Diehl
ab45818d51
seems to work now
...
anyway, nicer code
2018-08-22 14:30:51 +02:00