Hello,

I make a DG programme on libmesh, there is a problem that the first step
takes so long when I use
high order basis, the performance log is as follow:

------------------------------------------------------------------------------
| Matrix Assembly Performance: Alive time=683.179, Active time=682.44 |
------------------------------------------------------------------------------
| Event nCalls Total Avg Percent of |
| Time Time Active Time |
|------------------------------------------------------------------------------|
| |
| Assemble init 1 0.0001 0.000120 0.00 |
| bounray matrix 134 0.0041 0.000030 0.00 |
| elem init 1470 0.0757 0.000051 0.01 |
| interior matrix 4276 0.2926 0.000068 0.04 |
| mass matrix 5880 0.0607 0.000010 0.01 |
| matrix insertion 1470 173.4288 0.117979 25.41 |
| neighbor matrix insertion 4276 508.5776 0.118938 74.52 |
------------------------------------------------------------------------------
| Totals: 17507 682.4396 100.00 |
------------------------------------------------------------------------------

Why the matrix insertion spend more than 10 minutes even there is only
17460 dofs?
code like this:
perf_log.stop_event("interior matrix");

perf_log.start_event ("neighbor matrix insertion");

euler_system.matrix->add_matrix(Kn,dof_indices,neighbor_dof_indices);

perf_log.stop_event ("neighbor matrix insertion");

perf_log.start_event ("matrix insertion");

euler_system.matrix->add_matrix(Ke, dof_indices);
euler_system.rhs->add_vector(Fe, dof_indices);

perf_log.stop_event ("matrix insertion");



-------------------------------------------------------------------------
This SF.net email is sponsored by the 2008 JavaOne(SM) Conference 
Don't miss this year's exciting event. There's still time to save $100. 
Use priority code J8TL2D2. 
http://ad.doubleclick.net/clk;198757673;13503038;p?http://java.sun.com/javaone
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to