On 1/6/23 22:01, Wasim Niyaz Munshi ce21d400 wrote:

H_plus is being called in Assemble_damage and not assemble_elastic. It uses elastic solution, cell and gauss point to evaluate strain at a gauss point. Then some quantity is evaluated based on the strain.

Similarly I have another function damage_gauss which is being called in assemble_elastic that evaluates damage at a gauss point using the damage solution, cell and gauss point.

Wasim -- that is not helpful. You complain that the function takes a lot of time to run, but your description contains nothing that would anyone help understanding *why* that might be so. For all we know "then some quantity is evaluated based on the strain" could be an operation that solves a differential equation in itself.

Show us the function, or better even: Do some timings of that function yourself to figure out which line or code block within that function is expensive. Then we can think of improvements.

Best
 W.

--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 bange...@colostate.edu
                           www: http://www.math.colostate.edu/~bangerth/


--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/a94b3828-77e2-f3a9-5cb0-e1a4690a5d87%40colostate.edu.

Reply via email to