Commit Graph

1897 Commits

Author SHA1 Message Date
Martin Diehl 6170209198 use statement for marc was wrong (preprocessor) 2015-04-11 05:09:15 +00:00
Martin Diehl d7b36c2c30 gfortran complaint about implicit casting,floating point comparison, and unused imported variables.
additionally changed use of math_invert33 to math_inv33 if det is not needed
2015-04-10 19:09:26 +00:00
Martin Diehl 7d0a1ebea6 math_transpose33 not used 2015-04-08 18:33:08 +00:00
Martin Diehl b6481c2513 introduced simpler multiplication and forall loops
matmul is ok for openmp, check in the web and run the state integration test.

Example program testing for new state update for rkck dot state:

program test

real, dimension(6,10) :: dotState=reshape(&
                      [1,1,1,1,1,1,1,1,1,1,&
                       2,2,2,2,2,2,2,2,2,2,&
                       3,3,3,3,3,3,3,3,3,3,&
                       4,4,4,4,4,4,4,4,4,4,&
                       5,5,5,5,5,5,5,5,5,5,&
                       6,6,6,6,6,6,6,6,6,6],[6,10])

real, dimension(10) :: residuum
real, dimension(6) :: B=2.5

integer :: i

residuum = B(1)*dotState(1,:)+& 
           B(2)*dotState(2,:)+& 
           B(3)*dotState(3,:)+& 
           B(4)*dotState(4,:)+& 
           B(5)*dotState(5,:)+& 
           B(6)*dotState(6,:)
do i =1,10
  print*,residuum(i)
enddo

residuum =  matmul(transpose(dotState),B)
do i =1,10
  print*,residuum(i)
enddo

end program test
2015-04-01 16:45:53 +00:00
Martin Diehl 693efcaa58 comment was confused by search and replace 2015-04-01 16:04:33 +00:00
Martin Diehl 63491bf268 for some reasons, ifort was complaining about mixed integer types.
also deleted #ifdef PETSC since we do not support the solver without PETSC
2015-03-30 20:49:17 +00:00
Pratheek Shanthraj c7418db9bd FEM : hierarchical ordering of h5 output data and more meaningful visualisation for multiple grains/crystallites/phases 2015-03-30 09:45:10 +00:00
Martin Diehl e0f1132a17 gfortran complaints: equal comparison of reals and external (MPI) functions 2015-03-29 12:54:13 +00:00
Martin Diehl 9e824d029e renamed to new names: spectralsolver and spectralfilter 2015-03-28 21:14:11 +00:00
Martin Diehl 1b2e4c0239 removed old solver 2015-03-28 17:41:40 +00:00
Martin Diehl 8a30021202 cleaned up core modul 2015-03-27 14:10:18 +00:00
Pratheek Shanthraj c8554a5641 corrected alignment of restart output byte offset. restart test working now 2015-03-27 07:17:24 +00:00
Pratheek Shanthraj 1ee81e74ea more patches to get the tests running again:
- fixed increment counting in postResults to current output scheme

- corrected node coordinates calculation

- some work on restart (restart test still not running)
2015-03-26 21:19:28 +00:00
Pratheek Shanthraj 86283359f1 PetscFinalize should be the last call before exiting 2015-03-26 13:43:18 +00:00
Pratheek Shanthraj a9cf8e8051 some patches to compile damask core module without fftw-mpi 2015-03-26 13:41:19 +00:00
Pratheek Shanthraj 2417877bdb dropped non petsc basic spectral solver 2015-03-25 16:08:41 +00:00
Pratheek Shanthraj d44fce4a76 Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me 2015-03-25 16:06:19 +00:00
Pratheek Shanthraj 37a7364a3e multi physics output now working for all solvers 2015-03-25 16:02:30 +00:00
Pratheek Shanthraj 2d4e7fb8e7 calculate dFi/dS (9x9 matrix inverse) for the analytic jacobian only when needed. should improve performance 2015-03-20 07:41:11 +00:00
Pratheek Shanthraj 2e44a846af need to keep track of temperature again since we can no longer get this from F_T 2015-03-18 18:03:18 +00:00
Pratheek Shanthraj c5a8543f17 updated petsc 2015-03-18 17:18:43 +00:00
Martin Diehl 4a0c1c1717 comments 2015-03-15 15:30:14 +00:00
Martin Diehl 12cb5e56fe fixed interpretation 'spectralOut' header. Seem to have been buggy all the time, but manifested itself only after keyword change. Also added ':' after keywords to be more robust 2015-03-15 15:29:38 +00:00
Martin Diehl ac3d82493e changed output keywords form resolution to grid and from dimension to size 2015-03-13 14:45:33 +00:00
Martin Diehl e4d70e971d cp is commercially pure ;) 2015-03-13 13:49:09 +00:00
Martin Diehl d8debc3396 updated reference geom files (used old names), updated numerics.config to reflect use of analytic jaco 2015-03-13 06:47:26 +00:00
Martin Diehl aa125bac5b aliases for grid and size not supported any more 2015-03-12 22:30:24 +00:00
Martin Diehl 00cba25a44 improved update of ip coordinates for spectral solver.
do not need to create a new fftw plan all the time, using data already defined for the convolution
2015-03-12 22:28:33 +00:00
Martin Diehl 263f997cf2 seems to be an old output 2015-03-12 08:51:56 +00:00
Martin Diehl b338f1aeb8 added some outputs 2015-03-09 23:10:36 +00:00
Su Leen Wong 6b0a2ed5a4 Changed initialization of projection matrix for transformation based on the number of slip systems 2015-03-09 12:05:38 +00:00
Martin Diehl a5ee02ff01 added some aliases for output to finally converge from 'shear_rate' to 'shearrate' and from 'accumulatedshear' to 'accumulated_shear' 2015-03-09 11:29:44 +00:00
Franz Roters 2ab90c28a5 replaced example ODF by reasonable data in fully correct format 2015-03-06 15:56:58 +00:00
Pratheek Shanthraj 35f9e91e73 some corrections to the damage models 2015-03-06 13:12:32 +00:00
Pratheek Shanthraj bbb5ff6ae9 changes related to intermediate configuration kinematics:
- switched Fi and Li from state variables to crystallite variables

- Lp and Li are now work conjugate with the corresponding mandel stresses defined in their respective configuration

- T, Lp and Li need to return tangents wrt Fi arising from the convection of the material frame due to Fi

- Updated analytic jacobian to take into account tangents wrt Fi

- Updated Lp and Li residual jacobians to take into account tangents wrt Fi
2015-03-06 13:09:00 +00:00
Yannick Naunheim 4fd7e0165e more meaningful name 2015-03-05 11:07:55 +00:00
Su Leen Wong 36f14a2d1d Removed nr = ns condition 2015-02-26 14:16:29 +00:00
Su Leen Wong 2baac6fc10 Fix for 3881 because gfortran was complaining 2015-02-24 10:15:07 +00:00
Pratheek Shanthraj 3c2a95efd8 if matrix inversion error encountered during analytic tangent calculation, give a warning and return elastic tangent 2015-01-30 14:00:32 +00:00
Pratheek Shanthraj f1f4f939a3 analytic jacobian default for all solvers 2015-01-29 13:59:49 +00:00
Pratheek Shanthraj dd8458a775 updated analytic jacobian calculation to correctly take into account intermediate configuration Fi. improved convergence of Li loop in stress integration 2015-01-29 13:58:25 +00:00
Pratheek Shanthraj 0b59519a2a updated damage models:
- coupling to plasticity handled within damage module instead of plasticity module
- anisotropic models more stable
2015-01-29 13:56:09 +00:00
Martin Diehl 12c586afb9 removed unused variables 2015-01-23 13:44:16 +00:00
Martin Diehl 1f666dffc9 fixed bug: output size was not set 2015-01-23 13:08:25 +00:00
Martin Diehl e9c3137a3e disabled a faulty ntrans loop (matmul was complaining) 2015-01-23 11:23:35 +00:00
Martin Diehl 43095b1f1b added threshold_stress_slip as output (got lost once upon a time) 2015-01-22 09:34:11 +00:00
Martin Diehl 1c2dd37e28 deleted unused variables in dislo UCLA 2015-01-21 15:14:00 +00:00
David Mercier e57cb2b4f4 Wrong number of slip families in BCC material config files 2015-01-20 17:22:52 +00:00
Luv Sharma 19dc945065 Inroduced Fd in kinematic decomposition. 2015-01-19 13:39:47 +00:00
Pratheek Shanthraj 54dc4943de adaptive penalty to bound Cv in (0,1) 2015-01-16 18:31:14 +00:00