<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=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<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:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-ligatures:none;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 2.0cm 70.85pt;}
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" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal">Hello,<br>
<span style="font-size:12.0pt;color:black">i am using petsc in a single cpu setup where I have preassembled crs matrices that I wrap via PetSC’s<o:p></o:p></span></p>
<p class="MsoNormal">MatCreateSeqAIJWithArrays Functionality.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Now I manipulate the values of these matrices (wohtout changing the sparsity) without using petsc,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">When I now want to solve again I have to call <o:p></o:p></p>
<p class="MsoNormal">PetscObjectStateIncrease((PetscObject)petsc_A);<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">So that oetsc actually solves again (otherwise thinking nothing hs changed ,<o:p></o:p></p>
<p class="MsoNormal">This means I have to include the private header<o:p></o:p></p>
<p class="MsoNormal">#include <petsc/private/petscimpl.h><o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Which makes a seamingless implementation of petsc into a cmake process more complicate (This guy has to be stated explicitly in the cmake process at the moment)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I would like to resolve that by “going” around the private header,<o:p></o:p></p>
<p class="MsoNormal">My first intuition was to increase the state by hand<o:p></o:p></p>
<p class="MsoNormal">((PetscObject)petsc_A_aux[the_sys])->state++;<br>
<br>
<o:p></o:p></p>
<p class="MsoNormal">This is the definition of petscstateincrease in the header. This throws me an
<o:p></o:p></p>
<p class="MsoNormal">error: invalid use of incomplete type ‘struct _p_PetscObject’<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">compilation error.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Is there any elegeant way around this?<o:p></o:p></p>
<p class="MsoNormal">This is the first time I use the petsc mailing list so apologies for any beginners mistake I did in formatting or anything else.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Best regards<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Franz Pichler<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>