Paras, You could try to use SuperLU_MT (see https://portal.nersc.gov/project/sparse/superlu/) but we don't have wrapper for it. The Krylov solver in deal.II are multithreaded but the preconditioners are not. What you can try is to use deal.II solvers and STRUMPACK (https://github.com/pghysels/STRUMPACK) for the preconditioner but again we don't have wrapper for it.
Best, Bruno On Friday, February 14, 2020 at 7:22:15 AM UTC-5, Paras Kumar wrote: > > 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/2a2ce382-894d-42ff-a9ea-d1c15612a309%40googlegroups.com.
