Hello Jose.
Some months ago, I wrote to the petsc-users mailing list asking if it would
be possible to use an iterative solver in slepc which converges eigenvectors
in order of maximum projection on a given basis vector. Back then you told
me you might look into it in some time and let me know.
On Mon, Jul 16, 2012 at 11:44 PM, Matthew Knepley wrote:
> On Mon, Jul 16, 2012 at 3:47 PM, Thomas Hisch wrote:
>>
>> Hello list,
>>
>> is it possible to store just the nonzero pattern of a large (sparse)
>> matrix in a file? If I write the whole matrix to disk
>> (PetscViewerBinaryOpen/MatView)
--- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120717/d6df5ba7/attachment.html>
Great! Thank you, that was exactly what I was looking for! Could not
find it in the documentation, so I asked here.
Cheers,
-David
On Tue, Jul 17, 2012 at 4:50 PM, Satish Balay wrote:
> If PETSc is installed with hypre - then the corresponding petscconf.h
> will have a flag PETSC_HAVE_HYPRE defi
On 07/17/2012 04:04 PM, David Nies wrote:
> Hello all!
>
> Is there a way to find out if a given PETSc installation has been
> configured and built with '--with-hypre=1'? The background is that I
> want to write an M4 test for our local build system to see if we
> should support PETSc's HYPRE capab
Hello all!
Is there a way to find out if a given PETSc installation has been
configured and built with '--with-hypre=1'? The background is that I
want to write an M4 test for our local build system to see if we
should support PETSc's HYPRE capabilities. Of couse, I could test if
there is a 'libHYP
;> I modified your code slightly. All these options can be provided at runtime:
>>> '-ksp_type cg -pc_type bjacobi -sub_pc_type icc'
>>>
>>> Hong
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> On 07/13/2012 12:14 PM, Hong Zhang wrote:
>>>> Michele :
>>>>
>>>> I need to use the ICC factorization as preconditioner, but I noticed that
>>>> no parallel version is supported.
>>>> Is that correct?
>>>> Correct.
>>>>
>>>> If so, is there a work around, like building the preconditioner "by hand"
>>>> by using PETSc functions?
>>>> You may try block jacobi with icc in the blocks '-ksp_type cg -pc_type
>>>> bjacobi -sub_pc_type icc'
>>>>
>>>> Hong
>>>>
>>>
>>>
>>
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120717/31b6707d/attachment.html>
>> Please update to petsc-3.3. petsc-3.1 is too old.
> I would do that but the version installed on the platform (Intrepid at
> ALCF) I am working on is 3.1-p2.
Satish,
Please fix this.
Thanks
Barry
On Jul 17, 2012, at 1:36 PM, Michele Rosso wrote:
>
> On 07/17/2012
pc/impls/mg/mg.c
> >> [0]PETSC ERROR: PCSetUp_ML() line 668 in src/ksp/pc/impls/ml/ml.c
> >> [0]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c
> >> [0]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c
> >> [0]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c
> >>
> >> Regards
> >> Stijn
> >>
> >> Stijn A.M. Vantieghem
> >> Earth and Planetary Magnetism
> >> Institute for Geophysics
> >> ETH Z?rich
> >> Sonneggstrasse 5 - CH 8092 Z?rich
> >> tel: +41 44 632 39 90
> >> e-mail: stijn.vantieghem at erdw.ethz.ch
> >>
> >>
> >>
> >>
> >
>
>
>
>
>
> --
> 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
--
Ir. B. Sanderse
Centrum Wiskunde en Informatica
Science Park 123
1098 XG Amsterdam
t: +31 20 592 4161
e: sanderse at cwi.nl
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120717/10b73a22/attachment-0001.html>
On Jul 17, 2012, at 4:31 AM, TAY wee-beng wrote:
> Hi,
>
> I have 3 vectors of different dimension, create using 2 DM (da_dof1,da_dof2),
> with dof = 1 and 2.
>
> They are declared as:
>
> Vec u,v,duv
>
> call DMCreateLocalVector(da_dof1,u,ierr)
>
> call DMCreateLocalVector(da_dof1,v,ierr)
rus-free.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120717/75e52bcb/attachment.html>
-- next part --
A non-text attachment was scrubbed...
Name: not available
Type:
torization as preconditioner, but I noticed that
>> no parallel version is supported.
>> Is that correct?
>>
> Correct.
>
>
>> If so, is there a work around, like building the preconditioner "by
>> hand" by using PETSc functions?
>>
> Y
o which their
experiments lead.
-- Norbert Wiener
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120717/53dcf180/attachment.html>
If so, is there a work around, like building the
>> preconditioner "by hand" by using PETSc functions?
>>
>> You may try block jacobi with icc in the blocks '-ksp_type cg
>> -pc_type bjacobi -sub_pc_type icc'
>>
>> Hong
>>
>>
>
>
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120717/077b6027/attachment.html>
Hi,
I have 3 vectors of different dimension, create using 2 DM
(da_dof1,da_dof2), with dof = 1 and 2.
They are declared as:
Vec u,v,duv
call DMCreateLocalVector(da_dof1,u,ierr)
call DMCreateLocalVector(da_dof1,v,ierr)
call DMCreateLocalVector(da_dof2,duv,ierr)
How can I add duv to u and v t
bjacobi -sub_pc_type icc'
>
> Hong
>
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120717/2d7cd0fd/attachment.html>
Hello all,
I am trying to use ML to solve a Poisson equation with Neumann BC (singular
matrix) and found the thread below to avoid a zero pivot error. I used the
option that Barry suggests, -mg_coarse_pc_factor_shift_type nonzero, which
works, but only when I run on a single processor. For two
om your system and notify the sender immediately by
> return message.
>
> E-mail communication cannot be guaranteed to be timely secure, error or
> virus-free.
>
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any resul
If PETSc is installed with hypre - then the corresponding petscconf.h
will have a flag PETSC_HAVE_HYPRE defined.
[You could use this flag directly in your code to enable-disable
corresponding functionality.]
Satish
On Tue, 17 Jul 2012, David Nies wrote:
> Hello all!
>
> Is there a way to find
an that be done?
--
Yours sincerely,
TAY wee-beng
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120717/be63a2ed/attachment.html>
On Mon, Jul 16, 2012 at 11:35 PM, Satish Balay wrote:
> On Mon, 16 Jul 2012, Thomas Hisch wrote:
>
>> Is there any benefit to the user if he enables boost support in petsc
>> with the --with-boost configure option ? What is the difference
>> between a petsc build with enabled boost support and a b
; same error.
>
-mg_coarse_redundant_pc_factor_shift_type nonzero
Or use -mg_coarse_pc_type svd and forget about it (provided the coarsest
grid is not too big).
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments
0]PETSC ERROR: MatLUFactorNumeric_SeqAIJ() line 574 in
>> src/mat/impls/aij/seq/aijfact.c
>> >> [0]PETSC ERROR: MatLUFactorNumeric() line 2587 in
>> src/mat/interface/matrix.c
>> >> [0]PETSC ERROR: PCSetUp_LU() line 158 in
>> src/ksp/pc/impls/factor/lu/lu.c
>> >> [0]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c
>> >> [0]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c
>> >> [0]PETSC ERROR: PCSetUp_MG() line 602 in src/ksp/pc/impls/mg/mg.c
>> >> [0]PETSC ERROR: PCSetUp_ML() line 668 in src/ksp/pc/impls/ml/ml.c
>> >> [0]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c
>> >> [0]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c
>> >> [0]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c
>> >>
>> >> Regards
>> >> Stijn
>> >>
>> >> Stijn A.M. Vantieghem
>> >> Earth and Planetary Magnetism
>> >> Institute for Geophysics
>> >> ETH Z?rich
>> >> Sonneggstrasse 5 - CH 8092 Z?rich
>> >> tel: +41 44 632 39 90
>> >> e-mail: stijn.vantieghem at erdw.ethz.ch
>> >>
>> >>
>> >>
>> >>
>> >
>>
>>
>>
>
>
> --
> 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
>
>
> --
> Ir. B. Sanderse
>
> Centrum Wiskunde en Informatica
> Science Park 123
> 1098 XG Amsterdam
>
> t: +31 20 592 4161
> e: sanderse at cwi.nl
>
>
--
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/20120717/d06de09e/attachment-0001.html>
tely 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/20120717/9eda6375/attachment.html>
impls/ml/ml.c
> >> [0]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c
> >> [0]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c
> >> [0]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c
> >>
> >> Regards
> >> Stijn
> >>
> >> Stijn A.M. Vantieghem
> >> Earth and Planetary Magnetism
> >> Institute for Geophysics
> >> ETH Z?rich
> >> Sonneggstrasse 5 - CH 8092 Z?rich
> >> tel: +41 44 632 39 90
> >> e-mail: stijn.vantieghem at erdw.ethz.ch
> >>
> >>
> >>
> >>
> >
>
>
>
--
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/20120717/f5f7dc94/attachment.html>
imenters 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/20120717/7b730b64/attachment.html>
bled does not make any difference for me ?
>
Yes
Matt
--
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/20120717/243e94d4/attachment-0001.html>
27 matches
Mail list logo
