Dear deal.ii Community, I am working on finite deformation hyperelasticity problem which is essentially a nonlinear-vector-valued problem with displacement as the unknown at each support point(dim=2,3). *With regards to parallelism, we currently restrict ourselves to shared memory parallel (SMP) only.* The FE assembly process has been paralleized using the workstream function.
This question pertains to (possibly faster) linear solvers. Currently. For the current 2D problem, we use the SparseDirectUMFPACK solver. As can been seen in the time log below (ths example was computed on a 36 core machine with 64 threads), the linear solver consumes the most time. I also tried using the CG solver, but it was much much slower, probably due to the large condition number stemming out of highly refined mesh in regions of interest. Triangulation: Number of active cells: 2059646 Number of Degrees of Freedom: 4119454 +---------------------------------------------+------------+------------+ | Total wallclock time elapsed since start | 7.81e+03s | | | | | | | Section | no. calls | wall time | % of total | +---------------------------------+-----------+------------+------------+ | Assemble Linear System | 50 | 875s | 11% | | Linear Solver | 40 | 5.55e+03s | 71% | | Make Constraints | 50 | 255s | 3.3% | | PBC setup | 1 | 1.02s | 0% | | Set Boundary IDs | 1 | 0.109s | 0% | | Stress Output Computation | 11 | 464s | 5.9% | | Stress Output Overall | 12 | 1e+03s | 13% | | Stress Output Writing | 11 | 379s | 4.9% | | Stress output initialize | 1 | 158s | 2% | | System Setup | 1 | 260s | 3.3% | +---------------------------------+-----------+------------+------------+ While going through the documentation, I came across options for DMP linear solvers but couldn't find any with SMP. Could somenone please guide me to any such possibilities within the library or through interface to some external library. Best regards, Paras -- 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/b5508fa4-6a2c-4368-b727-52f5dd47f393%40googlegroups.com.