---
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/47f9e6c7/attachment.html>
Barry,
Are you able to reproduce the error?
I've updated petsc-dev it is still there unfortunately.
On 06.05.2012 20:00, Barry Smith wrote:
>Alexander,
>
> Sorry, never mind. I found the code in the earlier email.,
>
> barry
>
> On May 6, 2012, at 12:49 PM, Barry Smith wrote:
>
>> Al
mf,&dmc); // r=0,c=1
> >
> >
> > I would like a way to identify dmc as being on the "same level" as dm0.
>
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/dee975c2/attachment.html>
t ------
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/95b17000/attachment.html>
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-dev/attachments/20120506/75542bd0/attachment.html>
some similar ugly name.
>
PCMGSetGridSequencingCoarseLevel() ?
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/a5ffabb7/attachment.html>
You present a case where the one level method should be a direct solve. I
don't think this disproves my point that there is NO correct general abstract
or practical answer to what it should be.
Perhaps this should be a property of PCMG, if in the one level case it is a
direct solve or an
I think it belongs in the DM (or the matrix if it was put in the matrix).
When you select a subpart of the DM that new DM should contain the near null
space for the subpart. The DM associated with the field defining IS is where
all properties should go, not the IS itself. Of course, we reall
ing 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-dev/attachments/20120506/bc2840ad/attachment.html>
t; I would like a way to identify dmc as being on the "same level" as dm0.
>>
>>
>
--
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-dev/attachments/20120506/7c80c3bc/attachment.html>
near
null space is not a property of an IS at all.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/c063c54c/attachment.html>
attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/8d88c854/attachment.html>
hat to do for both of these separately, even
though the problem is of the same size in both cases, so it's usually
desirable for the same method to be used in both cases.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/ddb326a2/attachment.html>
is with 100 iterations.
Without valgrind the error is always there.
--
Regards,
Alexander
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/a7b674d7/attachment.html>
ot;same level" as dm0.
>
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/698faadf/attachment.html>
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-dev/attachments/20120506/ba21b6b1/attachment.html>
0287==by 0x4509F8: main (solveTest.c:66)
>
> Full log is attached.
>
> Important.
> If I comment this line:
> KSPComputeExtremeSingularValues(ksp, &maxx, &minx);
>
> It works.
>
> --
> Regards,
> Alexander
--
Regards,
Alexander
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/192a9cce/attachment.html>
Alexander,
I cannot reproduce this on my mac with 3 different blas/lapack.
Could you please run the case below but with --download-f-blas-lapack
(you forgot the -f last time)? Send us the valgrind results. This will tell use
the exact line number in dlasq3() that is triggering th
0287==by 0x6AA839: KSPDestroy (itfunc.c:780)
==20287==by 0x4509F8: main (solveTest.c:66)
Full log is attached.
Important.
If I comment this line:
KSPComputeExtremeSingularValues(ksp, &maxx, &minx);
It works.
--
Regards,
Alexander
-- next part --
An HTML
Fine
On May 6, 2012, at 2:44 PM, Jed Brown wrote:
> Should the refinement level be copied over by DMCoarsen (and the coarsen
> level be copied by DMRefine)?
>
> It's useful for diagnostics to be able to define a universal level. If I use
> PCMG and -snes_grid_sequence, there is effectively
t; Thanks.
>
> --
> Regards,
> Alexander
>
>
>
>
> --
> 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
--
Regards,
Alexander
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/b76d0547/attachment.html>
-- next part --
An embedded and charset-unspecified text was scrubbed...
Name: test.c
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/b76d0547/attachment.c>
pinit.
>
> I did not configure for a threaded version of PETSc.
>
> --
> Stefano
>
>
>
>
> --
> Stefano
>
>
>
>
> --
> Stefano
>
>
>
> --
> Stefano
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/5e8c4b85/attachment.html>
(dm0,&dmf); // r=1,c=0
DMCoarsen(dmf,&dmc); // r=0,c=1
I would like a way to identify dmc as being on the "same level" as dm0.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/b6e5c3fd/attachment.html>
Hello,
I use KSP and random rhs to compute largest singular value:
! create solver and set options for singular value estimation
call KSPCreate(MPI_COMM_WORLD,ksp,ierr);CHKERRQ(ierr)
call KSPSetType(ksp,KSPGMRES,ierr);CHKERRQ(ierr)
call
KSPSetTolerances(ksp,solvertol,PETSC_DEFAULT_DO
omm.c
>>> [0]PETSC ERROR: PetscRegisterFinalizeAll() line 313 in
>>> /work/adz/zampini/MyWorkingCopyOfPetsc/petsc-dev/src/sys/objects/destroy.c
>>> [0]PETSC ERROR: PetscFinalize() line 1080 in
>>> /work/adz/zampini/MyWorkingCopyOfPetsc/petsc-dev/src/sys/objects/pinit.
>>>
>>> I did not configure for a threaded version of PETSc.
>>>
>>> --
>>> Stefano
>>>
>>>
>>>
>>
>>
>> --
>> Stefano
>>
>>
>>
>
>
> --
> Stefano
>
--
Stefano
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/d44e46f2/attachment.html>
Alexander,
Sorry, never mind. I found the code in the earlier email.,
barry
On May 6, 2012, at 12:49 PM, Barry Smith wrote:
>
> Alexander,
>
> No attached test program. (only valgrind log).
>
> Barry
>
> On May 6, 2012, at 8:21 AM, Alexander Grayver wrote:
>
>> On 06.05.2012 15
Alexander,
No attached test program. (only valgrind log).
Barry
On May 6, 2012, at 8:21 AM, Alexander Grayver wrote:
> On 06.05.2012 15:07, Matthew Knepley wrote:
>>> Hello,
>>>
>>> I use KSP and random rhs to compute largest singular value:
>>>
>>> 1) Is this the whole program? If no
I see. So if the Mat is not known when PCSetFromOptions() is called then
PCSetFromOptions() generates a bad guess for the PC to use.
Presumably this has always been true.
Presumably the PCSetFromOptions() has access to the DM? Can it inquire of
the DM what PC to default to if the Mat
PetscErrorCode PCSetUp(PC pc)
{
PetscErrorCode ierr;
const char *def;
PetscFunctionBegin;
PetscValidHeaderSpecific(pc,PC_CLASSID,1);
if (!pc->mat)
SETERRQ(((PetscObject)pc)->comm,PETSC_ERR_ARG_WRONGSTATE,"Matrix must be set
first");
How can it be in there without a Mat? The Mat
On May 6, 2012, at 12:19 PM, Jed Brown wrote:
> On Sun, May 6, 2012 at 11:07 AM, Barry Smith wrote:
> Is your specific questions is: should we replace all MPI_Init() calls in
> PETSc with MPI_Init_thread(...,MPI_THREAD_FUNNELED,...)?
>
> Can't do this because of MPI-1 compatibility, but yeah
On May 6, 2012, at 11:09 AM, Jed Brown wrote:
> When threading is used, PETSc should be calling MPI_Init_thread() instead of
> MPI_Init(). For our current use case, MPI_THREAD_FUNNELED is the natural mode
> to request. As opposed to MPI_THREAD_SINGLE, this causes the MPI
> implementation to ma
L attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/89ec94fb/attachment.html>
scErrorCode (*f)(Mat,MatReuse,Mat*);
ierr =
PetscObjectQueryFunction((PetscObject)pc->pmat,"MatGetDiagonalBlock_C",(void
(**)(void))&f);CHKERRQ(ierr);
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/69dfe419/attachment.html>
the way up
here?
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/fab774b8/attachment.html>
pipermail/petsc-dev/attachments/20120506/e6501668/attachment.html>
actly the same thing in both cases.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/2deac2a9/attachment.html>
:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/1bbfb89e/attachment.html>
> Important.
> If I comment this line:
> KSPComputeExtremeSingularValues(ksp, &maxx, &minx);
>
> It works.
>
> --
>
> Regards,
> Alexander
>
>
>
> --
> Regards,
> Alexander
>
>
--
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-dev/attachments/20120506/69670179/attachment.html>
> KSPReset,FP=7fffd4765840
>> KSPDestroy, FP=7fffd47659a0
>> kspdestroy_, FP=7fffd47659d0
>>
>>
>> Any ideas?
>>
>> Thanks.
>>
>> --
>> Regards,
>> Alexander
>>
>>
>
>
> --
> 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
>
>
>
> --
> Regards,
> Alexander
>
>
--
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-dev/attachments/20120506/090a61ec/attachment.html>
0
> KSPReset,FP=7fffd4765840
> KSPDestroy, FP=7fffd47659a0
> kspdestroy_, FP=7fffd47659d0
>
>
> Any ideas?
>
> Thanks.
>
> --
> Regards,
> Alexander
>
>
--
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-dev/attachments/20120506/1f6fa4da/attachment.html>
40 matches
Mail list logo
