On Wed, Aug 8, 2012 at 9:12 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</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">
<div class="im"><br>
On Aug 8, 2012, at 9:08 PM, Satish Balay <<a href="mailto:balay@mcs.anl.gov">balay@mcs.anl.gov</a>> wrote:<br>
<br>
> Barry - you forgot to commit/push the interface files in<br>
> src/ksp/pc/impls/fieldsplit/ftn-custom.<br>
<br>
</div>    Thanks. I forgot to hg add the new files :-(<br>
<div class="im"><br>
><br>
> <a href="http://petsc.cs.iit.edu/petsc/petsc-dev/rev/40f9515d8b01" target="_blank">http://petsc.cs.iit.edu/petsc/petsc-dev/rev/40f9515d8b01</a><br>
><br>
> Perhaps this interface addition can go into petsc-3.3?<br>
<br>
</div>    No<br>
<div class="im"><br>
><br>
> That patch [with a minor fix] can be easily applied to petsc-3.1 to<br>
> get the link working. [old iffy functionality is a different matter :)]<br>
<br>
</div>   Users are free to shot themselves in the foot if they want to so Colin certainly could back-port this but I recommend upgrading.</blockquote><div><br></div><div>Colin, I also recommend upgrading since with the new stuff, you should be able to do everything you want</div>
<div>to do from the command line without messing around pulling out KSP objects.</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">
<span class="HOEnZb"><font color="#888888"><br>
   Barry<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
><br>
> Satish<br>
><br>
> On Wed, 8 Aug 2012, Barry Smith wrote:<br>
><br>
>><br>
>>  Well the field split stuff in 3.1 is rather iffy anyways.   You really need to take the couple hours and do the upgrade before trying to add new features to the fortran code you are using.  The time to upgrade will be much less than time spent in work arounds.<br>

>><br>
>>    Barry<br>
>><br>
>> On Aug 8, 2012, at 8:52 PM, Colin McAuliffe <<a href="mailto:cjm2176@columbia.edu">cjm2176@columbia.edu</a>> wrote:<br>
>><br>
>>> Is there an alternative to this routine that will allow me to define hierarchical field splits in 3.1p8? It will be a fairly significant undertaking to upgrade the fortran code I am using to be compatible with newer versions of petsc.<br>

>>><br>
>>> Colin<br>
>>><br>
>>> Quoting Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>>:<br>
>>><br>
>>>><br>
>>>>  Colin,<br>
>>>><br>
>>>>   I'm sorry we never had a FORTRAN interface for this routine.  You'll need to switch to petsc-dev  <a href="http://www.mcs.anl.gov/petsc/developers/index.html" target="_blank">http://www.mcs.anl.gov/petsc/developers/index.html</a> to use that  function.<br>

>>>><br>
>>>>   Barry<br>
>>>><br>
>>>><br>
>>>> On Aug 8, 2012, at 7:54 PM, Colin McAuliffe <<a href="mailto:cjm2176@columbia.edu">cjm2176@columbia.edu</a>> wrote:<br>
>>>><br>
>>>>> Thanks I will look into DM. Is PCFieldSplitGetSubKSP in petsc  version 3.1p8 callable from fortran? I am getting the following  error when I compile the fortran code, and I cant figure out if  this is an error in my code or what.<br>

>>>>><br>
>>>>> Undefined symbols for architecture x86_64:<br>
>>>>> "_pcfieldsplitgetsubksp_", referenced from:<br>
>>>>>    _umacr6_ in umacr6.o<br>
>>>>> ld: symbol(s) not found for architecture x86_64<br>
>>>>><br>
>>>>><br>
>>>>> thanks<br>
>>>>> Colin<br>
>>>>><br>
>>>>><br>
>>>>> Quoting Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>>:<br>
>>>>><br>
>>>>>> On Wed, Aug 8, 2012 at 8:48 AM, Colin McAuliffe  <<a href="mailto:cjm2176@columbia.edu">cjm2176@columbia.edu</a>>wrote:<br>
>>>>>><br>
>>>>>>> Hello all,<br>
>>>>>>><br>
>>>>>>> When using PCFieldSplitSetIS to define splits within splits, should the IS<br>
>>>>>>> contain indices in the original matrix or in the next highest split level?<br>
>>>>>>><br>
>>>>>><br>
>>>>>> The next highest split level. You might consider using a DM if you have<br>
>>>>>> some structure (like discretization on a mesh).<br>
>>>>>> This should handle all the index-wrangling for you. It is new, but intended<br>
>>>>>> for just this kind of thing.<br>
>>>>>><br>
>>>>>> Also, my matrix is in aij format but one of the sub fields has a block<br>
>>>>>>> diagonal structure. Is it still possible use block jacobi on this field?<br>
>>>>>>><br>
>>>>>><br>
>>>>>> Yes.<br>
>>>>>><br>
>>>>>> Matt<br>
>>>>>><br>
>>>>>><br>
>>>>>>> Thanks and all the best<br>
>>>>>>> Colin<br>
>>>>>>><br>
>>>>>>><br>
>>>>>>><br>
>>>>>>><br>
>>>>>>> --<br>
>>>>>>> Colin McAuliffe<br>
>>>>>>> PhD Candidate<br>
>>>>>>> Columbia University<br>
>>>>>>> Department of Civil Engineering and Engineering Mechanics<br>
>>>>>>><br>
>>>>>><br>
>>>>>><br>
>>>>>><br>
>>>>>> --<br>
>>>>>> What most experimenters take for granted before they begin their<br>
>>>>>> experiments is infinitely more interesting than any results to which their<br>
>>>>>> experiments lead.<br>
>>>>>> -- Norbert Wiener<br>
>>>>>><br>
>>>>><br>
>>>>><br>
>>>>><br>
>>>>> --<br>
>>>>> Colin McAuliffe<br>
>>>>> PhD Candidate<br>
>>>>> Columbia University<br>
>>>>> Department of Civil Engineering and Engineering Mechanics<br>
>>>><br>
>>>><br>
>>>><br>
>>><br>
>>><br>
>>><br>
>>> --<br>
>>> Colin McAuliffe<br>
>>> PhD Candidate<br>
>>> Columbia University<br>
>>> Department of Civil Engineering and Engineering Mechanics<br>
>><br>
>><br>
><br>
<br>
</div></div></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>