
<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>

</head>

<body dir="ltr">

<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

just in case it helps, here's one way I have to create the near nullspace:</div>

<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

<br>

</div>

<div class="elementToProof"><span style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);"><a href="https://github.com/seamplex/feenox/blob/main/src/pdes/mechanical/init.c#L468" id="LPlnk365748">https://github.com/seamplex/feenox/blob/main/src/pdes/mechanical/init.c#L468</a></span></div>

<div class="elementToProof"><br>

</div>

<div id="appendonsend"></div>

<div style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">

--</div>

<div style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">

jeremy</div>

<hr style="display: inline-block; width: 98%;">

<div id="divRplyFwdMsg" dir="ltr"><span style="font-family: Calibri, sans-serif; font-size: 11pt; color: rgb(0, 0, 0);"><b>From:</b> petsc-users <petsc-users-bounces@mcs.anl.gov> on behalf of Jordi Manyer Fuertes via petsc-users <petsc-users@mcs.anl.gov><br>

<b>Sent:</b> Tuesday, December 5, 2023 9:57 AM<br>

<b>To:</b> Matthew Knepley <knepley@gmail.com>; bsmith@petsc.dev <bsmith@petsc.dev><br>

<b>Cc:</b> petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov><br>

<b>Subject:</b> Re: [petsc-users] Help for MatNullSpaceCreateRigidBody</span>

<div> </div>

</div>

<div style="padding: 2pt; border-width: 1pt; border-style: solid; border-color: rgb(156, 101, 0);">

<p style="line-height: 12pt; background-color: rgb(255, 235, 156);"><span style="font-size: 10pt; color: red;"><b>[External Sender]</b></span></p>

</div>

<p>Thanks for the prompt response. Both answers look like what I'm doing.</p>

<p>After playing a bit more with solver, I managed to make it run in parallel with different boundary conditions (full dirichlet bcs, vs mixed newmann + dirichlet). This raises two questions:</p>

<p>- How relevant are boundary conditions (eliminating dirichlet rows/cols vs weak newmann bcs) to the solver? Should I modify something when changing boundary conditions?</p>

<p>- Also, the solver did well with the old bcs when run in a single processor (but not in parallel). This seems odd, since parallel and serial behavior should be consistent (or not?). Could it be fault of the PCGAMG? I believe the default local solver is ILU,

 shoud I be changing it to LU or something else for these kind of problems?</p>

<p>Thank you both again,</p>

<p>Jordi</p>

<p><br>

</p>

</body>

</html>