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.

Reply via email to