Martin Diehl
d7fa3a9791
just sorting
2018-10-18 21:34:26 +02:00
Martin Diehl
59399f2cd2
sanity checks
2018-10-18 19:06:21 +02:00
Martin Diehl
33d816a5d0
polishing
2018-10-17 23:56:57 +02:00
Martin Diehl
121d146a35
using more central functions
2018-10-17 23:20:39 +02:00
Martin Diehl
4320061554
using central functionality
2018-10-17 23:13:47 +02:00
Martin Diehl
8017fccbc3
adjusting to phenopowerlaw
2018-10-17 22:18:33 +02:00
Martin Diehl
2fd348f251
just re-ordering
2018-10-17 21:47:50 +02:00
Martin Diehl
6b4b5b511f
don't let dipole formation be undefined
2018-10-17 21:41:21 +02:00
Martin Diehl
b7bc5afa90
Merge remote-tracking branch 'origin/development' into 20-NewStyleDislotwin
2018-10-17 20:10:49 +02:00
Martin Diehl
7ac96bd630
polishing
2018-10-17 00:01:21 +02:00
Martin Diehl
d92ab04892
Merge branch 'development' into 19-NewStylePhenopowerlaw
2018-10-15 22:38:55 +02:00
Martin Diehl
be8d7e19fe
missing use statement caused compilation error
2018-10-15 22:38:34 +02:00
Martin Diehl
ad1a64c338
rename was missing
2018-10-15 06:01:38 +02:00
Martin Diehl
57ca93c3c9
added missing allocation
2018-10-15 05:28:52 +02:00
Martin Diehl
e5ec5ed54c
check prerequisites as early as possible
2018-10-15 05:03:53 +02:00
Martin Diehl
767ca0edd4
avoid code repetition
2018-10-14 20:16:30 +02:00
Martin Diehl
a311cc7673
Merge remote-tracking branch 'origin/development' into 19-NewStylePhenopowerlaw
2018-10-14 15:53:50 +02:00
Martin Diehl
6599aa487c
kinetics_slip does not seem to be the problem
2018-10-14 15:53:24 +02:00
Martin Diehl
1839228add
sorting
...
new mappings. Needs more work (don't forget https://xkcd.com/927 )
2018-10-14 12:27:42 +02:00
Martin Diehl
943fa764c0
not needed
2018-10-14 11:15:07 +02:00
Martin Diehl
4ac12d1dc5
removed hybridIA
2018-10-14 10:11:26 +02:00
Martin Diehl
4c780226d1
polishing
2018-10-14 09:27:34 +02:00
Martin Diehl
513faa2218
investigating the reason for the poor performance
2018-10-13 11:29:07 +02:00
Franz Roters
d9fbf2ab50
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-10-13 10:52:03 +02:00
Franz Roters
905b5a33eb
added link for Marc2018
2018-10-13 10:50:20 +02:00
Franz Roters
1a458108bc
added support for Marc2018
2018-10-13 10:26:42 +02:00
Martin Diehl
fa9cbe13c9
do concurrent not safe here
2018-10-13 08:35:36 +02:00
Martin Diehl
6207781eb6
[skip ci] [skip sc] streamlining
2018-10-12 17:24:46 +02:00
Franz Roters
d1561e867a
Merge branch '44-column-major-access-to-interaction-matrices' into 'development'
...
Resolve "column-major access to interaction matrices"
Closes #44
See merge request damask/DAMASK!43
2018-10-12 13:03:22 +02:00
Martin Diehl
5597ee338e
simplifying
...
- whole type is private, no need to make the components private
- target not needed when using associate
2018-10-12 08:26:24 +02:00
Martin Diehl
854d3ab774
clearer logic for controlling dipole formation
2018-10-12 08:23:56 +02:00
Vitesh Shah
38f9807c42
cleaning up the codes
2018-10-11 18:00:01 +02:00
Martin Diehl
51a3308a3d
bugfix: same tau was used for all systems
2018-10-11 17:28:00 +02:00
Vitesh Shah
cd2fb9f417
adding prefixes to writing of datasets under groups
2018-10-10 19:08:44 +02:00
Vitesh Shah
fa13e71493
removed typo(&)
2018-10-10 18:09:10 +02:00
Vitesh Shah
78db0ecb3c
reading and writing of subroutines of datasets of different ranks
2018-10-10 17:54:55 +02:00
Martin Diehl
7217cdac1b
Merge branch '29-rename-mesh_element-array' into development
2018-10-10 15:55:23 +02:00
Martin Diehl
ef1899cfd0
improved error reporting
...
- tell user specifically whether number of systems is too high or too
low
- c/a sanity check can be done centrally
2018-10-10 07:42:07 +02:00
Philip Eisenlohr
f0b9c0caf7
polishing and encapsulating of Abaqus-specific local variables
2018-10-09 17:57:10 -04:00
Franz Roters
1d0584dcdf
typo in message
2018-10-09 16:44:47 +02:00
Martin Diehl
d81c9f744e
avoid errors by consistently closing objects
...
general bugfixes, works now
2018-10-09 16:17:21 +02:00
Martin Diehl
2064ed80fd
more flexible file open routine
2018-10-09 15:15:08 +02:00
Martin Diehl
8626ba1ce6
adjusted variable names and added comments
2018-10-09 14:46:57 +02:00
Martin Diehl
48403c37b7
initializing and closing the HDF5 library
2018-10-09 14:13:51 +02:00
Martin Diehl
e7e0cf7118
preprocessor statements are case sensitive..
...
... and the PETSc variable is PETSc, not PETSC
2018-10-09 13:57:29 +02:00
Martin Diehl
e3943d282a
Merge branch 'development' into 47-save-restart-information-in-hdf5
2018-10-09 13:51:58 +02:00
Vitesh Shah
09ba95b029
non-executable
2018-10-09 10:58:45 +02:00
Vitesh Shah
70c746a8f1
Interfacing for subroutines to identify real and integer data
2018-10-09 10:57:06 +02:00
Martin Diehl
e71ec2057e
hotfix
...
does not seem to cause trouble here, but better be safe
2018-10-09 08:37:03 +02:00
Martin Diehl
4da3923440
more straight-forward definition of array shape
2018-10-08 23:54:45 +02:00
Martin Diehl
05aae4f7eb
removed forgotten debug statement
2018-10-08 23:34:05 +02:00
Martin Diehl
50f0362aa1
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2018-10-08 23:33:31 +02:00
Martin Diehl
6064daceda
do error check earlier
...
array out of bounds happened already before the check
2018-10-08 23:32:46 +02:00
Martin Diehl
b51c3b5a07
handling files without LF at last line correctly
2018-10-08 23:27:06 +02:00
Philip Eisenlohr
225a48021d
foramtting of random_seed output to align with rest
2018-10-08 17:01:56 -04:00
Martin Diehl
78b5fa466d
fixed hick-up in lattice definitions
2018-10-08 19:32:18 +02:00
Martin Diehl
71d80cc203
sanity check too strong
...
still makes sense to have, but some tests need to be updated before
setting it active again
2018-10-08 13:27:15 +02:00
Martin Diehl
2a79f7ea61
taking care of c/a for hex and bct
2018-10-08 13:25:40 +02:00
Martin Diehl
5efcad952a
avoid type conversion
...
more systematic checking of correct input parameters
2018-10-08 08:28:56 +02:00
Martin Diehl
302cf9b6c2
bugfix, now inline with code from initializeStructure
2018-10-07 22:16:18 +02:00
Martin Diehl
facb1c2407
unified error checking
2018-10-07 20:03:45 +02:00
Martin Diehl
0e17b17b90
using full precision
2018-10-07 18:40:02 +02:00
Martin Diehl
34e0aca564
using full precision
2018-10-07 18:18:24 +02:00
Martin Diehl
a53488d666
calculation of slip/twin/trans/damage-coordinate system was wrong
2018-10-07 15:44:13 +02:00
Martin Diehl
032c35a499
WIP: Calculation of elasticity matrices for twin and trans
2018-10-06 10:42:25 +02:00
Martin Diehl
3a2f86df1c
non-schmid is always bcc
2018-10-05 23:04:06 +02:00
Martin Diehl
81a77b508a
check for array size was off by 1
2018-10-05 22:46:53 +02:00
Martin Diehl
3a39c2d68c
[skip ci] strings will be converted to lowercase
2018-10-05 17:54:24 +02:00
Martin Diehl
0cddf36106
typo, using value of p_slip for q
2018-10-05 17:47:46 +02:00
Martin Diehl
2a51e528bd
Merge branch 'development' into 44-column-major-access-to-interaction-matrices
2018-10-05 11:11:20 +02:00
Martin Diehl
40cf7f134d
Merge branch '44-column-major-access-to-interaction-matrices' of magit1.mpie.de:damask/DAMASK into 44-column-major-access-to-interaction-matrices
2018-10-05 11:08:29 +02:00
Martin Diehl
aa34b83978
avoid array out of bounds
...
allows to handle un- or partly defined non schmid coefficients
2018-10-05 10:36:44 +02:00
Martin Diehl
a560fff2ac
[skip ci] [skip sc] generic interfaces for HDF5 operations
...
most existing HDF5 functions are only designed for writing output files
and more general functionality is needed for storing the restart data
2018-10-05 09:03:20 +02:00
Martin Diehl
47a748831f
Cannot write out resolved stress including non schmid
...
Would need to write out tau_pos and tau_neg separately, so limit to
resolved stress in the sense of Schmid
2018-10-05 07:26:13 +02:00
Martin Diehl
6ccbc5e4e6
characteristing shear is defined per system
2018-10-05 07:20:51 +02:00
Martin Diehl
e74a852ef7
new style is noticable slower
...
runtime (polycrystal example) increased from 70s to 90s.
Need to figure out how to write readable AND fast code
2018-10-05 06:38:20 +02:00
Martin Diehl
3a00fcde16
WIP: characteristic shear calculated in lattice
2018-10-05 05:09:31 +02:00
Martin Diehl
fda2676365
calculate only if required
2018-10-05 04:54:47 +02:00
Martin Diehl
79b0efc678
also using nonSchmid matrix from lattice
2018-10-04 21:20:42 +02:00
Martin Diehl
35956ce0c8
needs to be multiplikative
2018-10-04 21:06:57 +02:00
Vitesh Shah
a6ebdfc1ff
Reading and writing not very efficient method
2018-10-04 17:00:24 +02:00
Martin Diehl
edb94b8bc9
using the same kinetics for slip everywhere
2018-10-04 16:12:14 +02:00
Martin Diehl
12d053d123
strange parameter caused simulation not to start if not defined
...
probably a leftover. now only bothers you if twin systems are defined
2018-10-04 14:54:00 +02:00
Martin Diehl
f1750b4eef
avoiding use of unitialized arrays
2018-10-04 14:51:32 +02:00
Martin Diehl
cd7f285686
prevent use of unallocated array
...
interaction matrices are used even if no slip/twin systems defined
2018-10-04 07:17:00 +02:00
Martin Diehl
dc289a278b
clearer code
...
homogenization is defined per element, not per IP
hence, use material_homogenizationAt instead of deprecated
material_homog which pretends a dependency on the integration point
2018-10-04 06:48:33 +02:00
Martin Diehl
df473302f4
consistent (and better understandable) names
2018-10-04 06:03:48 +02:00
Martin Diehl
f494803c76
fixed tangent calculation
2018-10-03 16:40:40 +02:00
Martin Diehl
ec77d4450b
checking for required things
2018-10-03 09:17:06 +02:00
Martin Diehl
db645efda7
function to calculate nonSchmid matrix
2018-10-03 08:49:23 +02:00
Martin Diehl
103b46ecbb
don't repeate code over different plastic_xxx modules
2018-10-03 07:58:02 +02:00
Martin Diehl
89e03be699
also building interactions for slip-twin and twin-slip
2018-10-03 07:51:11 +02:00
Martin Diehl
0cac6ce0d2
allocating to wrong size
2018-10-03 07:36:57 +02:00
Martin Diehl
ba383336f5
only one projection matrix required for all (non)schmid projections
...
following the example from plastic_nonlocal
2018-10-03 07:23:18 +02:00
Martin Diehl
d718521890
polishing
2018-10-03 07:13:30 +02:00
Martin Diehl
67c55965e9
segmentation fault, tried LHS with scalar
...
Nxx per system is stored in LATTICE_yyy_NxxxSYSTEM
Using CAPITALS for parameters improves readbility
2018-10-03 06:45:42 +02:00
Martin Diehl
97d2c2b353
mesh_element(1) has a size of 4
...
fixed sigsegv
2018-10-02 22:54:44 +02:00
Martin Diehl
740325165d
Merge branch 'development' into 20-NewStyleDislotwin
2018-10-02 20:37:11 +02:00
Martin Diehl
0047e6c14c
correct calculation of dGdot_dTau
...
needs to take drag coefficient into account for the derivative
2018-10-02 18:09:57 +02:00
Martin Diehl
9cc3a77e01
mean velocity of dislocation corrected for higher temperatures
...
- v_wait: obstacle spacing over waiting time (effective velocity)
- v_run: 'true' velocity of a dislocation moving between obstacles
for higher temperatures v_run becomes important and bounds the velocity
2018-10-02 15:56:18 +02:00
Martin Diehl
5d1e648d17
drag coefficient B introduced and read in
2018-10-02 15:25:40 +02:00
Martin Diehl
384a785805
more sevice functions use
...
still trying
2018-10-01 23:02:31 +02:00
Martin Diehl
ff2614b757
using 'service functions' from lattice
...
get rid of many global array on the long run
2018-10-01 22:48:14 +02:00
Martin Diehl
a943940760
Merge branch 'development' into 44-column-major-access-to-interaction-matrices
2018-10-01 22:32:38 +02:00
Martin Diehl
98296a978e
Merge branch 'development' into 29-rename-mesh_element-array
2018-10-01 22:29:31 +02:00
Martin Diehl
e7b8adadd0
passing in instance and of is much easier
...
also using always the full Mandel stress tensor to avoid repeated
conversion
2018-10-01 22:14:54 +02:00
Martin Diehl
ddfbe19f02
Merge branch 'development' into 20-NewStyleDislotwin
2018-10-01 21:16:11 +02:00
Martin Diehl
27bde05529
Merge branch '38-introduce-rudimentary-PETSc-based-FEM-solver' into 'development'
...
Resolve "Introduce Rudimentary PETSc based FEM solver"
Closes #38
See merge request damask/DAMASK!41
2018-10-01 21:10:03 +02:00
Franz Roters
d5cf7d2144
Merge branch '19-NewStylePhenopowerlaw' into 'development'
...
Resolve "New coding style for phenopowerlaw"
See merge request damask/DAMASK!39
2018-10-01 16:08:05 +02:00
Martin Diehl
49a4202e26
not needed here
2018-09-29 10:31:45 +02:00
Martin Diehl
14509876d4
error == 0 is normal termination
2018-09-29 09:57:53 +02:00
Martin Diehl
637544509e
rewind needed
...
the file is read twice from beginning to end
2018-09-28 10:06:21 +02:00
Vitesh Shah
8619b49e20
made it nonexecutable
2018-09-28 08:47:00 +02:00
Vitesh Shah
8bd0cb3b25
Made it non executable
2018-09-28 08:45:52 +02:00
Vitesh Shah
6efde5f4ce
Combined the dataset creation and writing dataset in single subroutine. Created group for looped writing.
2018-09-28 08:41:29 +02:00
Martin Diehl
b141252f1f
doing the same things in the same way in both solvers
2018-09-28 08:06:22 +02:00
Martin Diehl
d1f614991e
merging good style mutually FEM <-> Spectral
2018-09-28 07:49:52 +02:00
Martin Diehl
6780217193
adopted from DAMASK_spectral
2018-09-28 07:25:32 +02:00
Martin Diehl
e8f687a99c
typo, Abaqus/Marc did not compile
2018-09-27 20:18:37 +02:00
Martin Diehl
822bc4dbbf
Merge branch 'development' into 29-rename-mesh_element-array
2018-09-27 20:12:09 +02:00
Martin Diehl
4c057ba529
both solvers can share quit
2018-09-27 20:09:59 +02:00
Chuanlai Liu
b391ad304d
tensor input/output is of the form 11,12,13, 21,22,23, 31,32,33. Thus row index i is slow, while column index j is fast
2018-09-27 11:53:45 +02:00
Martin Diehl
6509775d41
cleaning
...
avoiding compiler warnings and confusion of the user
2018-09-27 08:23:30 +02:00
Vitesh Shah
5b52ac91e5
Revert "made it non executable"
...
This reverts commit 05c784aefc
.
2018-09-26 09:28:47 +02:00
Martin Diehl
c3380aa9f3
using wrong (uninitialized) stress
2018-09-26 09:22:12 +02:00
Vitesh Shah
2edf8f7f53
Need to commit this before reverting the changes in other files
2018-09-26 09:21:53 +02:00
Vitesh Shah
3408282755
windows and linux editor incompatibility
2018-09-25 16:49:25 +02:00
Vitesh Shah
05c784aefc
made it non executable
2018-09-25 16:48:02 +02:00
Vitesh Shah
b71896d4b6
writing of multidimensional dataset into hdf5 does not appear as multiple pages in hdfviewer
2018-09-25 16:42:43 +02:00
Martin Diehl
ceb385ef39
calcMode not needed for spectral and FEM
2018-09-23 21:31:30 +02:00
Martin Diehl
c42eb87a33
using arrays with new names
2018-09-23 20:53:35 +02:00
Martin Diehl
72b87b0a9b
better name
2018-09-23 19:58:43 +02:00
Martin Diehl
f701910c71
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-23 19:32:19 +02:00
Martin Diehl
ebef12e446
syntax errors in declaration
2018-09-23 19:04:17 +02:00
Martin Diehl
67483487ea
more descriptive names
2018-09-23 18:50:54 +02:00
Martin Diehl
854d99250c
mesh_element(1,:) only used for debug output
...
set to -1 at the moment to indicate that it is not used.
Re-implementation should be done for MSC.Marc and Abaqus only.
2018-09-23 18:42:49 +02:00
Martin Diehl
4b14cc5560
calcmode only needed for Abaqus and MSC.Marc
2018-09-23 18:06:18 +02:00
Martin Diehl
cf6d388a6b
consistent order of solver-specific functionality
2018-09-23 17:57:48 +02:00
Martin Diehl
5936397ae7
introducing better names
...
allows further simplifications as we do not store max and per elem
values any more for number of integration points and number of cell
nodes
2018-09-23 17:37:57 +02:00
Martin Diehl
2fe2c4ca45
leaner syntax with sourced allocation
2018-09-23 17:26:13 +02:00
Martin Diehl
51390b1acf
Nelems /= NcpElems only in case of Abaqus/MSC.Marc
2018-09-23 17:05:01 +02:00
Martin Diehl
5814e07021
simplified
2018-09-23 16:55:03 +02:00
Martin Diehl
e3f2ad34b2
not needed
2018-09-23 16:32:13 +02:00
Martin Diehl
bd60018513
not needed for spectral
2018-09-23 15:57:21 +02:00
Martin Diehl
c1b8854132
only needed for commercial solvers
2018-09-23 15:31:19 +02:00
Martin Diehl
4862aca340
grouping solver specific variables for better readability
2018-09-23 15:27:51 +02:00
Martin Diehl
db45b7615a
drop support for heterogeneous meshes
...
heterogeneous meshes are neither advisable nor typically used
2018-09-23 15:24:14 +02:00
Martin Diehl
519cd840bd
cleaning
2018-09-23 15:19:23 +02:00
Pratheek Shanthraj
109ed4308f
change in tabulation order. should now be working correctly
2018-09-23 03:14:23 +02:00
Martin Diehl
9eb318b9bf
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 14:54:58 +02:00