Dear all,
I am experiencing long times when computing the assembling and I would like
to ask if this is common or there is something wrong with my implementation.
My model is built in a similar way as step-29 and step-40 (using complex
values ad solving with a direct solver using distributed parallel
implementation).
Now I am running larger systems with 3.5million dof and the assembling took
16h, while the solver function took much less.
I can show the structure of my assembly_system() function to ask if there
is something that can be done in order to speed up the process:
*void Problem<dim>::assemble_system()*
*{*
* for (unsigned int i = 0; i < dofs_per_cell; ++i) {*
* for (unsigned int j = 0; j < dofs_per_cell; ++j)*
* {*
* for (unsigned int q_point = 0; q_point < n_q_points; ++q_point)*
* { *
* if (fe.system_to_component_index(i).first ==
fe.system_to_component_index(j).first)*
* {*
* cell_matrix(i, j) += ....*
* }*
* if (fe.system_to_component_index(i).first !=
fe.system_to_component_index(j).first)*
* {*
* cell_matrix(i, j) += ....*
* }*
* }*
// Boundaries
* if (fe.system_to_component_index(i).first ==
fe.system_to_component_index(j).first)*
* {*
* for (unsigned int face_no : GeometryInfo<dim>::face_indices())*
* if (cell->face(face_no)->at_boundary() &&
(cell->face(face_no)->boundary_id() == 0))*
* {*
* fe_face_values.reinit(cell, face_no);*
* for (unsigned int q_point = 0; q_point < n_face_q_points; ++q_point)*
* cell_matrix(i, j) += ....*
* }*
* }*
* if (fe.system_to_component_index(i).first !=
fe.system_to_component_index(j).first)*
* {*
* for (unsigned int face_no : GeometryInfo<dim>::face_indices())*
* {*
* if (cell->face(face_no)->at_boundary() &&
(cell->face(face_no)->boundary_id() == 0))*
* {*
* fe_face_values.reinit(cell, face_no);*
* for (unsigned int q_point = 0; q_point < n_face_q_points;++q_point)*
* cell_matrix(i, j) += ....*
* }*
* }*
* }*
* } *
*}*
Thank you very much.
Regards,
Hermes
--
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/299f4da5-cf2f-470f-8076-c759c092549fn%40googlegroups.com.
