I have a system of coupled heat equations, written for MPI using LinearAlgebraTrilinos::MPI::SparseMatrix system_matrix;
as the sparse system matrix. Now for certain values the program runs fine, but for other values (which I am interested in) the RAM consumption goes up rapidly after several steps in time and refinement. Valgrind tells me that the biggest reason for that is the system matrix. This makes the program unusable workstations with 32 GByte RAM. Is that something I could expect, or is that a bug in my program? If it is something I could expect: What are possible solutions? I am solving the equations with the method shown in step-15 and step-33, i.e. I iterate over an inner loop for every time step until I am below a certain l2-value, and then go to the next time point. One possible approach I could think of is to set the threshold higher, but that would only delay the problem. Another approach would be to check if the l2-value did not change over the last five internal loops, and then continue anyway, but that would again increase the error. Is there a better way? Thanks! -- 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. For more options, visit https://groups.google.com/d/optout.
