<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Hello,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">We are using PETSC in our CFD code, and noticed that using “MatCreateMPIAIJWithSplitArrays” is almost 60% faster for large problem size (i.e DOF > 725M, using GAMG each time-step only takes 5sec, compared to 8.3 sec when assembling the
matrix one row at a time using matsetvaluesblocked() as recommended). <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The problem is that the memory usage goes up after each call to MatCreateMPIAIJWithSplitArrays to update the matrix values. As MatCreateMPIAIJWithSplitArrays is not supposed to copy the values, do we need to call it each time to update
the values? We tried to just update the values of the diagonal and off-diagonal part of the arrays passed to “MatCreateMPIAIJWithSplitArrays”, (the sparsity structure is fixed) but it looks like that the values are not updated, what is the proper way to update
the values of the matrix created by MatCreateMPIAIJWithSplitArrays?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thank you<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" width="800" style="border-collapse:collapse">
<tbody>
<tr>
<td width="800" colspan="3" style="width:600.0pt;background:white;padding:.75pt .75pt .75pt .75pt">
<p class="MsoNormal" style="page-break-after:avoid;text-autospace:none"><b><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#00201F">Hassan Raiesi,
<o:p></o:p></span></b></p>
</td>
</tr>
<tr>
<td width="800" colspan="3" style="width:600.0pt;background:white;padding:.75pt .75pt .75pt .75pt">
</td>
</tr>
<tr>
<td width="800" colspan="3" style="width:600.0pt;background:white;padding:.75pt .75pt .75pt .75pt">
<p class="MsoNormal" style="page-break-after:avoid;text-autospace:none"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#00201F">Advanced Aerodynamics Department<o:p></o:p></span></p>
</td>
</tr>
<tr>
<td width="800" colspan="3" style="width:600.0pt;background:white;padding:.75pt .75pt .75pt .75pt">
<p class="MsoNormal" style="page-break-after:avoid;text-autospace:none"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#00201F">Bombardier Aerospace<o:p></o:p></span></p>
<p class="MsoNormal" style="page-break-after:avoid;text-autospace:none"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#00201F"><o:p> </o:p></span></p>
</td>
</tr>
<tr>
<td width="800" colspan="3" style="width:600.0pt;background:white;padding:.75pt .75pt .75pt .75pt">
<p class="MsoNormal" style="page-break-after:avoid;text-autospace:none"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#00201F">hassan.raiesi@aero.bombardier.com<o:p></o:p></span></p>
</td>
</tr>
<tr>
<td width="800" colspan="3" style="width:600.0pt;background:white;padding:.75pt .75pt .75pt .75pt">
<p class="MsoNormal" style="page-break-after:avoid;text-autospace:none"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#00201F"><o:p> </o:p></span></p>
</td>
</tr>
<tr>
<td width="800" colspan="3" style="width:600.0pt;background:white;padding:.75pt .75pt .75pt .75pt">
<p class="MsoNormal" style="page-break-after:avoid;text-autospace:none"><i><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#00201F">2351 boul. Alfred-Nobel (BAN1)<o:p></o:p></span></i></p>
</td>
</tr>
<tr>
<td width="800" colspan="3" style="width:600.0pt;background:white;padding:.75pt .75pt .75pt .75pt">
<p class="MsoNormal" style="page-break-after:avoid;text-autospace:none"><i><span lang="FR-CA" style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#00201F">Ville Saint-Laurent, Québec, H4S 2A9<o:p></o:p></span></i></p>
</td>
</tr>
<tr>
<td width="24" style="width:.25in;background:white;padding:.75pt .75pt .75pt .75pt">
<p class="MsoNormal" style="page-break-after:avoid;text-autospace:none"><i><span lang="FR-CA" style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#00201F"><o:p> </o:p></span></i></p>
</td>
<td width="576" valign="top" style="width:6.0in;padding:.75pt .75pt .75pt .75pt">
<p class="MsoNormal" style="page-break-after:avoid;text-autospace:none"><i><span lang="FR-CA" style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#00201F"><o:p> </o:p></span></i></p>
</td>
<td style="border:none;padding:0in 0in 0in 0in" width="200">
<p class="MsoNormal"> </p>
</td>
</tr>
<tr>
<td width="24" style="width:.25in;background:white;padding:.75pt .75pt .75pt .75pt">
<p class="MsoNormal" style="page-break-after:avoid;text-autospace:none"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#00201F">Tél.<o:p></o:p></span></p>
</td>
<td width="576" style="width:6.0in;background:white;padding:.75pt .75pt .75pt .75pt">
<p class="MsoNormal" style="page-break-after:avoid;text-autospace:none"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#00201F"> 514-855-5001 # 62204<o:p></o:p></span></p>
</td>
<td style="border:none;padding:0in 0in 0in 0in" width="200">
<p class="MsoNormal"> </p>
</td>
</tr>
<tr>
<td width="600" colspan="2" style="width:6.25in;background:white;padding:.75pt .75pt .75pt .75pt">
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#00201F"><o:p> </o:p></span></p>
</td>
<td style="border:none;padding:0in 0in 0in 0in" width="200">
<p class="MsoNormal"> </p>
</td>
</tr>
<tr>
<td width="600" colspan="2" style="width:6.25in;background:white;padding:.75pt .75pt .75pt .75pt">
<p class="MsoNormal" style="page-break-after:avoid;text-autospace:none"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#00201F"><img width="162" height="34" id="Picture_x0020_1" src="cid:image001.png@01D1D2C3.1D169550"><o:p></o:p></span></p>
</td>
<td style="border:none;padding:0in 0in 0in 0in" width="200">
<p class="MsoNormal"> </p>
</td>
</tr>
<tr>
<td width="600" colspan="2" style="width:6.25in;background:white;padding:.75pt .75pt .75pt .75pt">
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b><span style="font-size:8.5pt;font-family:"Arial",sans-serif;color:#555555"><o:p> </o:p></span></b></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b><span style="font-size:8.5pt;font-family:"Arial",sans-serif;color:#555555">CONFIDENTIALITY NOTICE</span></b><span style="font-size:8.5pt;font-family:"Arial",sans-serif;color:#666666">
- This communication may contain privileged or confidential information.<br>
If you are not the intended recipient or received this communication by error, please notify the sender<br>
and delete the message without copying, forwarding and/or disclosing it.<o:p></o:p></span></p>
</td>
<td style="border:none;padding:0in 0in 0in 0in" width="200">
<p class="MsoNormal"> </p>
</td>
</tr>
</tbody>
</table>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:8.0pt;font-family:"Arial",sans-serif;color:#5F5F5F"><o:p> </o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>