On Wed, Dec 7, 2011 at 2:54 AM, Parker, Andrew (UK Filton) <span dir="ltr"><<a href="mailto:Andrew.Parker2@baesystems.com">Andrew.Parker2@baesystems.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Legends, nice one!<br>
<br>
Got to that point in the souce code yesterday, saw that 7 was missing<br>
from the switch statement, and got worried.<br>
<br>
Is there any symmetry between the code in block size 5->6->7->8? I've<br>
asked that without looking, so could be stupid question. The reason I<br>
ask is that if there is this is a classic use of something like<br>
boost-preprocessor macros or a bit of templating, as in the compiler<br>
would expand it all out given a pattern at compile time? Just<br>
interest....<br></blockquote><div><br></div><div>This is the kind of thing you could do with templates. However, that would mean</div><div><br></div><div> 1) It only works in C++, and the logic is convoluted</div><div><br>
</div><div> 2) You can't see the generated code</div><div><br></div><div> 3) Its impossible to debug (and usually optimize)</div><div><br></div><div>We are looking at other approaches.</div><div><br></div><div> Thanks,</div>
<div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Cheers again,<br>
Andy<br>
<br>
-----Original Message-----<br>
From: <a href="mailto:petsc-users-bounces@mcs.anl.gov">petsc-users-bounces@mcs.anl.gov</a><br>
[mailto:<a href="mailto:petsc-users-bounces@mcs.anl.gov">petsc-users-bounces@mcs.anl.gov</a>] On Behalf Of Satish Balay<br>
Sent: 06 December 2011 19:31<br>
To: PETSc users list<br>
Subject: Re: [petsc-users] MATBDIAG in > 3.1<br>
<br>
*** WARNING ***<br>
<br>
This message has originated outside your organisation,<br>
either from an external partner or the Global Internet.<br>
Keep this in mind if you answer this message.<br>
<br>
<br>
On Tue, 6 Dec 2011, Matthew Knepley wrote:<br>
<br>
> On Tue, Dec 6, 2011 at 7:43 AM, Parker, Andrew (UK Filton) <<br>
> <a href="mailto:Andrew.Parker2@baesystems.com">Andrew.Parker2@baesystems.com</a>> wrote:<br>
><br>
> > So having got petsc 3.2 install and all its deps, I now get<br>
> > this....****<br>
> ><br>
> > ** **<br>
> ><br>
> > [0]PETSC ERROR: PCSetUp_PBJacobi() line 207 in<br>
> > src/ksp/pc/impls/pbjacobi/pbjacobi.c not supported for block size<br>
> > 7****<br>
> ><br>
> > ** **<br>
> ><br>
> > ** **<br>
> ><br>
> > Any thoughts?<br>
> ><br>
><br>
> Sorry about that. There was confusion on our end. I have added the<br>
> code you need to 3.2. It will go out with the next patch release.<br>
> Alternatively, you can retrieve the source tarball from the repository<br>
<br>
> site or use the repository directly.<br>
<br>
<br>
I'm attaching this change as a patchfile - which can be applied to the<br>
current<br>
petsc-3.2 sources you have.<br>
<br>
cd petsc-3.2<br>
patch -Np1 < bpjacobi.patch<br>
<br>
Satish<br>
<br>
********************************************************************<br>
This email and any attachments are confidential to the intended<br>
recipient and may also be privileged. If you are not the intended<br>
recipient please delete it from your system and notify the sender.<br>
You should not copy it or use it for any purpose nor disclose or<br>
distribute its contents to any other person.<br>
********************************************************************<br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>