On Thu, Jun 30, 2011 at 5:11 PM, Barry Smith wrote:
>
> On Jun 29, 2011, at 10:35 PM, John Fettig wrote:
>
>> On Wed, Jun 29, 2011 at 11:21 PM, Barry Smith wrote:
>>> ?If you pass in SAME_PRECONDITIONER to the KSPSetOperators() it is simply
>>> impossible that it will rebuild the preconditioner.
Aldo,
This is a bug. The allocated number was not multipled by the block size
squared. Because of the this the compute unneeded number was garbage. We fixed
this a while ago in petsc-dev (yes it should have been fixed as a patch for 3.1
but that didn't happen.)
Barry
On Jun 30, 2
On Jun 29, 2011, at 10:35 PM, John Fettig wrote:
> On Wed, Jun 29, 2011 at 11:21 PM, Barry Smith wrote:
>> If you pass in SAME_PRECONDITIONER to the KSPSetOperators() it is simply
>> impossible that it will rebuild the preconditioner. There must be another
>> explanation as to what is going o
fests when I try to create
>>>>>>> the
>>>>>>> larger matrix from the individual smaller ones. When the matrix is
>>>>>>> constructed I'm trying to use MatSetValues and having each node pass in
>>>>>>> it&#x
the program interpreting that as a 12x12 instead of as a
>>>>>> 6x6 with the blocks in a grid.
>>>>>>
>>>>>> My question is then: is it possible to fill a matrix as a grid of
>>>>>> blocks, or can I only fill it in groups of rows or columns? Also, am I
>>>>>> approaching this problem the correct way, or are there more efficient
>>>>>> ways
>>>>>> of building this matrix with the ultimate goal of inverting it?
>>>>>>
>>>>>> I have included my copy of a modified example if it helps. I do
>>>>>> apologize if this is answered somewhere in the documentation, I have been
>>>>>> unable to find a solution.
>>>>>>
>>>>>> Respectfully,
>>>>>> Adam
>>>>>>
>>>>>
>>>>>
>>>>>
>>>>> --
>>>>> What most experimenters take for granted before they begin their
>>>>> experiments is infinitely more interesting than any results to which their
>>>>> experiments lead.
>>>>> -- Norbert Wiener
>>>>>
>>>>
>>>>
>>>
>>>
>>> --
>>> What most experimenters take for granted before they begin their
>>> experiments is infinitely more interesting than any results to which their
>>> experiments lead.
>>> -- Norbert Wiener
>>>
>>
>>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110630/4064bee5/attachment.htm>
ger matrix from the individual smaller ones. When the matrix is
>>>>>>> constructed I'm trying to use MatSetValues and having each node pass in
>>>>>>> it's
>>>>>>> 3x3 block. I end up with an error that the sum
Dear all,
I have a doubt concerning the interpretation of the output from -vec_view.
I am using ILU(k) on a BAIJ matrix; the run is uni-processor.
> KSP Object:
> type: gmres
> GMRES: restart=30, using Classical (unmodified) Gram-Schmidt
> Orthogonalization with no iterative refinement
>
On Wed, Jun 29, 2011 at 11:35 PM, John Fettig wrote:
>> ? Is this only happing with PCML and not with other preconditioners like
>> PCILU?
>
> I have only tried PCML. ?Let me try PCILU and get back to you.
PCILU doesn't work at all on the second timestep in either version.
The solution hasn't me
