<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">Dear Barry,<br>
<br>
Many thanks for your reply.<br>
<br>
I checked with the valgrind and did not obtain any outputs (<code
class="filename">massif.out.<pid></code>) because the job
was killed before it reached the end.<br>
<br>
Then I switch to a smaller case, it works well and one of the
output is attached (I did not find any useful information in it).
The output with the option -mat_coloring_view is followed, which
shows that the number of colors is 65. Any ideas for this?<br>
<br>
MatColoring Object: 480 MPI processes<br>
type: sl<br>
Weight type: RANDOM<br>
Distance 2, Max. Colors 65535<br>
Number of colors 65<br>
Number of total columns 1637350<br>
<br>
Best regards,<br>
Rongliang<br>
<br>
On 02/17/2016 01:13 AM, Barry Smith wrote:<br>
</div>
<blockquote
cite="mid:684F7147-4F58-42E3-9C9E-2C2297383276@mcs.anl.gov"
type="cite">
<pre wrap="">
How many colors are needed?
You need to produce a breakdown of where all the memory is being used. For example valgrind with the <a class="moz-txt-link-freetext" href="http://valgrind.org/docs/manual/ms-manual.html">http://valgrind.org/docs/manual/ms-manual.html</a>
Barry
</pre>
<blockquote type="cite">
<pre wrap="">On Feb 16, 2016, at 6:45 AM, Rongliang Chen <a class="moz-txt-link-rfc2396E" href="mailto:rongliang.chan@gmail.com"><rongliang.chan@gmail.com></a> wrote:
Dear all,
I am using the DMPlex to solve a PDE on a unstructured mesh and I use the SNESComputeJacobianDefaultColor to compute the Jacobian matrix.
My code works well for small problems (such as problem with 3.3x10^5 cells using 120 cores) but when I increase the number of cells (2.6x10^6 cells using 1920 cores), I go out of memory in the function MatColoringApply. It shows that one of the cores uses over 11G memory. I think this is unreasonable. Do you have any suggestions for debugging this problem?
Best regards,
Rongliang
</pre>
</blockquote>
<pre wrap="">
</pre>
</blockquote>
<br>
</body>
</html>