n.nl
From: Satish Balay
Sent: Friday, September 25, 2020 9:02 AM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] something wrong with digest?
I see one of the mailman settings is:
How big in Kb should a digest be before it gets sent out? 0 implies no maximum
Today I got more than 20 petsc-user Digest emails in my inbox (Vol 141, Issue
78 to 101), many containing only one message and being sent within a few
minutes of each other (10:39, 10:41 10:44, ...). Is this really how Digest is
supposed to work?
Chris
dr. ir. Christiaan Klaij | Senior
propeller<https://www.marin.nl/news/new-publication-on-noise-measurements-on-a-cavitating-propeller>
From: Zhang, Hong
Sent: Friday, September 11, 2020 4:57 PM
To: Klaij, Christiaan; petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] How to activate the modified Gram
Make me feel ancient. Would be nice to have the fortran binding though...
Chris
> --
>
> Message: 1
> Date: Thu, 10 Sep 2020 19:41:30 -0600
> From: Zhuo Chen
> To: "Zhang, Hong"
> Cc: "petsc-users@mcs.anl.gov"
> Subject: Re:
euvring-vessel-trials-held-on-North-Sea.htm>
From: Matthew Knepley
Sent: Thursday, April 11, 2019 2:16 PM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] Error with KSPSetUp and MatNest
On Thu, Apr 11, 2019 at 7:51 AM Klaij, Christiaan v
Matt,
As a happy MATNEST user, I got a bit worried when you wrote "we
should remove MatCreateNest() completely". This happened last
time I tried to use AIJ instead:
https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2016-August/029973.html
Has this problem been fixed in the meantime?
Chris
>
18 4:50 PM
To: Klaij, Christiaan; petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] PetscPartitioner is missing for fortran
Thanks, Chris, for your kind reply. Yes, I believe the mailing list is the most
convenient way for Petsc's to communicate with each other.
Also, I am wondering if you may ha
Hi Jiaoyan,
I've been using the fortran interface since 2004. During these
years I only found a handful of things missing. After an email to
this list the developers are usually quite fast in fixing the
problem.
Chris
Date: Mon, 26 Nov 2018 16:41:29 +
From: Jiaoyan Li
To:
r 19, 2018 10:26 PM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] fortran INTENT with petsc object gives segfault
after upgrade from 3.8.4 to 3.10.2
> On Oct 19, 2018, at 9:37 AM, Klaij, Christiaan wrote:
>
> As far as I (mis)understand fortran, this is a da
news:
http://www.marin.nl/web/News/News-items/Seminar-Scheepsbrandstof-en-de-mondiale-zwavelnorm-2020.htm
From: Smith, Barry F.
Sent: Friday, October 19, 2018 2:32 PM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] fortran IN
I've recently upgraded from petsc-3.8.4 to petsc-3.10.2 and was
surprised to find a number of segmentation faults in my test
cases. These turned out to be related to user-defined MatMult and
PCApply for shell matrices. For example:
SUBROUTINE systemMatMult(aa_system,xx_system,rr_system,ierr)
ws/News-items/MARIN-at-WindDays-2018-June-13-14-Rotterdam.htm>
From: Matthew Knepley <knep...@gmail.com>
Sent: Wednesday, January 18, 2017 4:13 PM
To: Klaij, Christiaan
Cc: Lawrence Mitchell; petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] monitoring the
With petsc-3.7.5, I had F90 code like this:
CALL MatSeqAIJSetPreallocation(aa_symmetric,PETSC_NULL_INTEGER,d_nnz,ierr);
CHKERRQ(ierr)
CALL
MatMPIAIJSetPreallocation(aa_symmetric,PETSC_NULL_INTEGER,d_nnz,PETSC_NULL_INTEGER,o_nnz,ierr);
CHKERRQ(ierr)
which worked fine. Now, with
<ba...@mcs.anl.gov>
Sent: Thursday, April 05, 2018 3:31 PM
To: Klaij, Christiaan
Cc: Gaetan Kenway; petsc-users
Subject: Re: [petsc-users] --download-hdf5 and then compiling cgns
Glad it works..
Perhaps its is simple fix to add in --download-cgns option to petsc configure..
[I see its a direct depen
november naar de open dag in
Wageningen!<http://www.marin.nl/web/News/News-items/Kom-zaterdag-10-november-naar-de-open-dag-in-Wageningen.htm>
From: Gaetan Kenway <gaet...@gmail.com>
Sent: Thursday, March 29, 2018 6:43 PM
To: petsc-users
Cc: Klaij, C
Satish,
I've build petsc-3.8.3 with --download-hdf5 and this works
fine. Now, I'm trying to compile cgns-3.3.1 and have it use the
hdf5 library from petsc. Therefore, in the cmake configure I have
set this option
-D HDF5_DIR=/path/to/petsc/hdf5
but cmake gives this message:
CMake Warning at
odel-for-area-risk-assessment-of-shipsource-oil-spills-in-Canadian-waters-1.htm
From: Satish Balay <ba...@mcs.anl.gov>
Sent: Thursday, March 22, 2018 4:03 PM
To: Klaij, Christiaan
Cc: petsc-users
Subject: Re: [petsc-users] problem configuring 3.8.3
anl.gov>
Sent: Thursday, March 22, 2018 3:14 PM
To: Klaij, Christiaan
Cc: petsc-users
Subject: Re: [petsc-users] problem configuring 3.8.3 with intel mpi
Sure - they have so many ways of tweaking things - so one can't things
will be the same across multiple users.
Most users dont' set these env
, March 22, 2018 2:54 PM
To: Klaij, Christiaan
Cc: petsc-users
Subject: Re: [petsc-users] problem configuring 3.8.3 with intel mpi
On Thu, 22 Mar 2018, Klaij, Christiaan wrote:
> OK, configure works with the option --with-cc=/path/to/mpicc,
> thanks!
>
> I'm not sure you are right abo
018-Shanghai.htm
________
From: Klaij, Christiaan
Sent: Thursday, March 22, 2018 2:43 PM
To: petsc-users
Subject: Re: [petsc-users] problem configuring 3.8.3 with intel mpi
OK, configure works with the option --with-cc=/path/to/mpicc,
thanks!
I'm not sure you are right about checkin
18 2:18 PM
To: petsc-users
Cc: Klaij, Christiaan
Subject: Re: [petsc-users] problem configuring 3.8.3 with intel mpi
Added this change to balay/mpi-dir-check-warning/maint and merged to next.
Satish
On Thu, 22 Mar 2018, Satish Balay wrote:
> The relevant change that is causing the difference
p://www.marin.nl/web/News/News-items/Office-tool-aNySIM-integrated-into-realtime-onboard-software-1.htm>
________
From: Matthew Knepley <knep...@gmail.com>
Sent: Thursday, March 22, 2018 1:27 PM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] problem con
GROW-partners-innovate-together-in-offshore-wind-industry.htm
From: Satish Balay <ba...@mcs.anl.gov>
Sent: Thursday, December 07, 2017 6:07 PM
To: Klaij, Christiaan
Cc: petsc-users
Subject: Re: [petsc-users] segfault after recent scientific linux upgrade
t-Marintec-China-Shanghai-December-58-1.htm>
From: Kong, Fande <fande.k...@inl.gov>
Sent: Thursday, December 07, 2017 4:26 PM
To: Klaij, Christiaan
Cc: petsc-users
Subject: Re: [petsc-users] segfault after recent scientific linux upgrade
On Thu, Dec 7, 2017 at 8:15 AM, Klaij, C
.marin.nl
MARIN news:
http://www.marin.nl/web/News/News-items/Simulator-facility-in-Houston-as-bridge-between-engineering-and-operations.htm
________
From: Klaij, Christiaan
Sent: Thursday, December 07, 2017 12:02 PM
To: petsc-users
Cc: Fande Kong
Subject:
ghai-December-58-1.htm
From: Satish Balay <ba...@mcs.anl.gov>
Sent: Wednesday, December 06, 2017 6:05 PM
To: Klaij, Christiaan
Cc: Fande Kong; petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] segfault after recent scientific linux upgrade
petsc 3.7 - an
geningen>
MARIN news: Seminar ‘Blauwe toekomst: versnellen van innovaties door
samenwerken<http://www.marin.nl/web/News/News-items/Seminar-Blauwe-toekomst-versnellen-van-innovaties-door-samenwerken.htm>
From: Fande Kong <fdkong...@gmail.com>
Sent: Tuesday
I'm running production software with petsc-3.7.5 and, among
others, superlu_dist 5.1.3 on scientific linux 7.4.
After a recent update of SL7.4, notably of the kernel and glibc,
we found that superlu is somehow broken. Below's a backtrace of a
serial example. Is this a known issue? Could you
.htm
From: Jed Brown <j...@jedbrown.org>
Sent: Thursday, August 24, 2017 6:01 PM
To: Klaij, Christiaan; Matthew Knepley; Barry Smith
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] Petsc ILU PC Change between 3.6.4 and 3.7.x?
"Klaij, Christiaan" <c.kl...@ma
dition number of A. If instead of A you
use P A, which we assume has a lower condition number, then
this bound is improved.
Thanks,
Matt
Barry
> On Aug 22, 2017, at 11:27 PM, Klaij, Christiaan
> <c.kl...@marin.nl<mailto:c.kl...@marin.nl>> wrote:
>
> Barry,
>
etsc-users] Petsc ILU PC Change between 3.6.4 and
> 3.7.x?
> Message-ID: <14d4c158-95b1-96ab-90b7-b81d3239b...@imperial.ac.uk>
> Content-Type: text/plain; charset=utf-8
>
>
>
> On 22/08/17 14:49, Klaij, Christiaan wrote:
> > CALL
> > PetscOption
tems/BlueWeek-October-911-Rostock.htm
From: Barry Smith <bsm...@mcs.anl.gov>
Sent: Tuesday, August 22, 2017 6:25 PM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] Petsc ILU PC Change between 3.6.4 and 3.7.x?
> On Aug 2
We also faced this problem in our code. So I've added:
CALL
PetscOptionsSetValue(PETSC_NULL_OBJECT,"-sub_pc_factor_shift_type","nonzero",ierr)
since there seems to be no setter function for this (correct me
if I'm wrong). Then everythings fine again.
Out of curiosity, what was the reason to
er.com/MARIN_nieuws> [Facebook]
<https://www.facebook.com/marin.wageningen>
MARIN news: Meet us again at the OTC
2017<http://www.marin.nl/web/News/News-items/Meet-us-again-at-the-OTC-2017.htm>
From: Matthew Knepley <knep...@gmail.com>
Sent: Tues
o A11, it does
reach the Schur complement and I do get convergence, see
attachment "out2".
Chris
From: Barry Smith <bsm...@mcs.anl.gov>
Sent: Monday, March 27, 2017 6:35 PM
To: Klaij, Christiaan
Cc: Lawrence Mitchell; Matthew Knepley; petsc-u
-Rotterdam.htm
From: Barry Smith <bsm...@mcs.anl.gov>
Sent: Saturday, March 25, 2017 1:29 AM
To: Klaij, Christiaan
Cc: Lawrence Mitchell; Matthew Knepley; petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] left and right preconditioning with a constant null
space
>
| Senior Researcher | Research & Development
MARIN | T +31 317 49 33 44 | mailto:c.kl...@marin.nl | http://www.marin.nl
MARIN news:
http://www.marin.nl/web/News/News-items/Meet-us-again-at-the-OTC-2017.htm
____
From: Klaij, Christiaan
Sent: Friday, March
http://www.marin.nl
MARIN news:
http://www.marin.nl/web/News/News-items/Comfort-and-Safety-at-Sea-March-29-Rotterdam.htm
From: Klaij, Christiaan
Sent: Friday, March 24, 2017 9:05 AM
To: Lawrence Mitchell; Matthew Knepley
Cc: petsc-users@mcs.anl.gov
Subject: Re: [p
list-Human-Performance.htm
From: Lawrence Mitchell <lawrence.mitch...@imperial.ac.uk>
Sent: Thursday, March 23, 2017 11:57 AM
To: Klaij, Christiaan; Matthew Knepley
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] left and right preco
nghai for the first time, March
28<http://www.marin.nl/web/News/News-items/Software-seminar-in-Shanghai-for-the-first-time-March-28.htm>
From: Matthew Knepley <knep...@gmail.com>
Sent: Wednesday, March 22, 2017 4:47 PM
To: Klaij, Christiaan
Cc: petsc-users@m
I'm solving the Navier-Stokes equations using PCFieldSplit type
Schur and Selfp. This particular case has only Neumann conditions
for the pressure field.
With left preconditioning and no nullspace, I see that the KSP
solver for S does not converge (attachment "left_nonullsp") in
either norm.
Barry,
I came across the same problem and decided to use KSPSetNormType
instead of KSPSetPCSide. Do I understand correctly that CG with
KSP_NORM_UNPRECONDITIONED would be as efficient as with
KSP_NORM_PRECONDITIONED? Since PC_RIGHT is not supported, I was
under the impression that the former
7 10:59 AM
To: petsc-users@mcs.anl.gov
Cc: bsm...@mcs.anl.gov; Klaij, Christiaan
Subject: Re: [petsc-users] monitoring the convergence of fieldsplit 0 and 1
On 18/01/17 08:40, Klaij, Christiaan wrote:
> Barry,
>
> I've managed to replicate the problem with 3.7.4
> snes/examples/tutorials/
MARIN | T +31 317 49 33 44 | mailto:c.kl...@marin.nl | http://www.marin.nl
MARIN news:
http://www.marin.nl/web/News/News-items/Few-places-left-for-Offshore-and-Ship-hydrodynamics-courses.htm
________
From: Klaij, Christiaan
Sent: Tuesday, January 17, 2017 8:
dsplit_1_inner_pc_type","jacobi",ierr);
CHKERRQ(ierr)
dr. ir. Christiaan Klaij | CFD Researcher | Research & Development
MARIN | T +31 317 49 33 44 | mailto:c.kl...@marin.nl | http://www.marin.nl
MARIN news:
http://www.marin.nl/web/News/News-items/Verification-and-validati
testing-continues.htm
From: Barry Smith <bsm...@mcs.anl.gov>
Sent: Friday, January 13, 2017 7:51 PM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] monitoring the convergence of fieldsplit 0 and 1
Yes, I would have exp
n.nl
MARIN news:
http://www.marin.nl/web/News/News-items/MARIN-wishes-you-a-challenging-inspiring-2017.htm
From: Barry Smith <bsm...@mcs.anl.gov>
Sent: Saturday, April 11, 2015 12:27 AM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject:
v>
Sent: Thursday, January 05, 2017 7:02 PM
To: Matthew Knepley
Cc: Klaij, Christiaan; petsc-users
Subject: Re: [petsc-users] problems after glibc upgrade to 2.17-157
On Thu, 5 Jan 2017, Matthew Knepley wrote:
> On Thu, Jan 5, 2017 at 2:37 AM, Klaij, Christiaan <c.kl...@marin.nl>
From: Satish Balay <ba...@mcs.anl.gov>
Sent: Wednesday, January 04, 2017 7:24 PM
To: petsc-users
Cc: Klaij, Christiaan
Subject: Re: [petsc-users] problems after glibc upgrade to 2.17-157
On Wed, 4 Jan 2017, Satish Balay wrote:
> So I guess your best bet is static libraries..
Or upgrad
From: Klaij, Christiaan
Sent: Wednesday, January 04, 2017 3:53 PM
To: Matthew Knepley
Cc: petsc-users; Satish Balay
Subject: Re: [petsc-users] problems after glibc upgrade to 2.17-157
So how would I do that? Does LIBS= accept spaces in the string?
Something like this perhaps:
LIBS="
rediction.htm>
From: Matthew Knepley <knep...@gmail.com>
Sent: Wednesday, January 04, 2017 3:13 PM
To: Klaij, Christiaan
Cc: petsc-users; Satish Balay
Subject: Re: [petsc-users] problems after glibc upgrade to 2.17-157
On Wed, Jan 4, 2017 at 7:37 AM, Klaij, Christ
<http://www.marin.nl/web/News/News-items/Modelling-natural-transition-on-hydrofoils-for-application-in-underwater-gliders-1.htm>
From: Matthew Knepley <knep...@gmail.com>
Sent: Wednesday, January 04, 2017 1:40 PM
To: Klaij, Christiaan
Cc: petsc-users; Sa
317 49 33 44 | mailto:c.kl...@marin.nl | http://www.marin.nl
MARIN news:
http://www.marin.nl/web/News/News-items/Software-seminar-in-Shanghai-for-the-first-time-March-28.htm
________
From: Klaij, Christiaan
Sent: Wednesday, January 04, 2017 9:26 AM
To: Matthew Knepley;
thew Knepley <knep...@gmail.com>
Sent: Tuesday, January 03, 2017 4:36 PM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] problems after glibc upgrade to 2.17-157
On Tue, Jan 3, 2017 at 9:04 AM, Klaij, Christiaan
<c.kl...@marin.nl<mailto:c.kl...@marin.nl>>
317 49 33 44 | mailto:c.kl...@marin.nl | http://www.marin.nl
MARIN news:
http://www.marin.nl/web/News/News-items/SSSRIMARIN-seminar-November-2-Shanghai.htm
From: Jed Brown <j...@jedbrown.org>
Sent: Monday, October 31, 2016 3:26 PM
To: Klaij, Christi
Satish,
I've noticed that SuperLU depends on metis and parmetis and that
PETSc downloads the versions 5.1.0-p3 and 4.0.3-p3. These are
different from the Karypis latest stable versions (without the
-p3). Do I really need these -p3 versions?
If so, after configure, compilation and installation by
eren-15-november.htm
From: Satish Balay <ba...@mcs.anl.gov>
Sent: Wednesday, October 26, 2016 5:32 PM
To: petsc-users
Cc: Klaij, Christiaan
Subject: Re: [petsc-users] petsc 3.7.4 with superlu_dist install problem
One additional note - there is this op
Here is a small complaint about the error message "unable to
download" that is given when using
--download-PACKAGENAME=/PATH/TO/package.tar.gz with the wrong
tarball. For example, for my previous install, I was using
petsc-3.5.3 with:
--download-ml=/path/to/ml-6.2-win.tar.gz
Using the same file
Jed,
Just a reminder, you haven't responded to this thread, notably
Matt's question below whether you can fix the L2G mapping.
Chris
> > From: Matthew Knepley <knep...@gmail.com>
> > Sent: Tuesday, August 02, 2016 12:28 AM
> > To: Klaij, Christiaan
> > Cc: pets
no
MatNestSetPreallocation, should one get the submats and
preallocate those if type=nest?
Chris
> Date: Tue, 2 Aug 2016 08:49:36 -0500
> From: Matthew Knepley <knep...@gmail.com>
> To: "Klaij, Christiaan" <c.kl...@marin.nl>
> Cc: "petsc-users@mcs.anl.gov" <petsc
.@gmail.com>
> Sent: Tuesday, August 02, 2016 12:28 AM
> To: Klaij, Christiaan
> Cc: petsc-users@mcs.anl.gov; Jed Brown
> Subject: Re: [petsc-users] block matrix without MatCreateNest
>
> On Mon, Aug 1, 2016 at 9:36 AM, Klaij, Christiaan <c.kl...@marin.nl> wrote:
>
>
___
From: Matthew Knepley <knep...@gmail.com>
Sent: Monday, August 01, 2016 4:15 PM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] block matrix without MatCreateNest
On Mon, Aug 1, 2016 at 8:59 AM, Klaij, Christiaan
<c.kl...@marin.nl<mailto:
Matt,
Thanks for your suggestions. Here's the outcome:
1) without the "if type=nest" protection, I get a "Cannot
locate function MatNestSetSubMats_C" error when using
type mpiaij, see below.
2) with the isg in a proper array, I get the same "Invalid
Pointer to Object" error, see below.
Chris
Matt, Barry
Thanks for your replies! I've added a call to MatNestSetSubMats()
but something's still wrong, see below.
Chris
$ cat mattry.F90
program mattry
use petscksp
implicit none
#include
PetscInt :: n=4 ! setting 4 cells per process
PetscErrorCode :: ierr
PetscInt
Anyone?
(my guess is an if-statement, something like "if type nest then
setup nest"...)
> Date: Thu, 28 Jul 2016 08:38:54 +0000
> From: "Klaij, Christiaan" <c.kl...@marin.nl>
> To: "petsc-users@mcs.anl.gov" <petsc-users@mcs.anl.gov>
> Sub
I'm trying to understand how to assemble a block matrix in a
format-independent manner, so that I can switch between types
mpiaij and matnest.
The manual states that the key to format-independent assembly is
to use MatGetLocalSubMatrix. So, in the code below, I'm using
this to assemble a 3-by-3
petsc-3.5.3 after patch
balay@asterix /home/balay/download-pine/x
$ ./ex1f
0
0
balay@asterix /home/balay/download-pine/x
Satish
$ On Thu, 30 Apr 2015, Klaij, Christiaan wrote:
Satish,
Thanks for the files! I've copied them to the correct location
: Friday, April 24, 2015 10:54 PM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] PETSC_NULL_OBJECT gets corrupt after call to
MatNestGetISs in fortran
Sorry for dropping the ball on this issue. I pushed the following fix to maint
branch.
https://bitbucket.org/petsc
a
PETSC_NULL_OBJECT argument.
Satish, could you please see if you can add a custom fortran stub for
this function in maint?
Thanks
Barry
On Feb 12, 2015, at 3:02 AM, Klaij, Christiaan C.Klaij at marin.nl wrote:
Using petsc-3.5.3, I noticed that PETSC_NULL_OBJECT gets corrupt after
calling
Barry, Matt
Thanks, the total iteration count is already very useful. Looking forward to
the custom monitors as well.
Chris
From: Barry Smith bsm...@mcs.anl.gov
Sent: Saturday, April 11, 2015 12:27 AM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Using petsc-3.5.3, I noticed that PETSC_NULL_OBJECT gets corrupt after calling
MatNestGetISs in fortran. Here's a small example:
$ cat fieldsplittry2.F90
program fieldsplittry2
use petscksp
implicit none
#include finclude/petsckspdef.h
PetscErrorCode :: ierr
PetscInt ::
Hi Fabian,
After reading your thread, I think we have the same
discretization. Some thoughts:
- The stabilization by pressure-weighted interpolation should
lead to a diagonally dominant system that can be solved by
algebraic multigrid (done for example in CFX), see
Hi Dave,
My understanding is that stabilization by pressure weighted
interpolation in FVM is in essence a fourth order pressure
derivative with proper scaling. But in FVM it is usually
implemented in defect correction form: only part of the
stabilization is in the A11 block of the coupled matrix,
In the program below, I'm using PCFieldSplitGetSubKSP to get the
sub KSP's of a Schur fieldsplit preconditioner. I'm setting
fieldsplit_0 to BICG+ILU and fieldsplit_1 to CG+ICC. Running
$ ./fieldsplittry -ksp_view
shows that this works as expected (full output below). Now, I
would like to change
To: Klaij, Christiaan; petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] fieldsplit_0_ monitor in combination with selfp
On Fri, Sep 5, 2014 at 1:34 AM, Klaij, Christiaan
c.kl...@marin.nlmailto:c.kl...@marin.nl wrote:
Matt,
I think the problem is somehow related to -pc_fieldsplit_schur_precondition
knep...@gmail.com
Sent: Wednesday, September 03, 2014 5:02 PM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] PCFieldSplitSetSchurPre in fortran
On Wed, Sep 3, 2014 at 9:46 AM, Klaij, Christiaan
c.kl...@marin.nlmailto:c.kl...@marin.nl wrote:
print
I'm playing with the selfp option in fieldsplit using
snes/examples/tutorials/ex70.c. For example:
mpiexec -n 2 ./ex70 -nx 4 -ny 6 \
-ksp_type fgmres \
-pc_type fieldsplit \
-pc_fieldsplit_type schur \
-pc_fieldsplit_schur_fact_type lower \
-pc_fieldsplit_schur_precondition selfp \
(1,0) : prefix=a10_, type=mpiaij, rows=24, cols=48
(1,1) : prefix=fieldsplit_1_, type=mpiaij, rows=24, cols=24
From: Matthew Knepley knep...@gmail.com
Sent: Thursday, September 04, 2014 2:20 PM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] fieldsplit_0_
I'm using PCFieldSplitGetSubKSP to access and modify the various
KSP's that occur in a Schur preconditioner. Now,
PCFieldSplitGetSubKSP returns two KSP's: one corresponding to A00
and one corresponding to the Schur complement.
However, in the LDU factorization, ksp(A00) occurs 4 times: in
the
immediately and delete your copy from your system.
From: Matthew Knepley knep...@gmail.com
Sent: Tuesday, September 02, 2014 5:42 PM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] PCFieldSplitSetSchurPre in fortran
On Tue, Sep 2, 2014
...
Chris
From: Matthew Knepley knep...@gmail.com
Sent: Wednesday, September 03, 2014 2:12 PM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] PCFieldSplitSetSchurPre in fortran
On Wed, Sep 3, 2014 at 7:00 AM, Klaij, Christiaan
c.kl
, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] PCFieldSplitSetSchurPre in fortran
On Wed, Sep 3, 2014 at 9:23 AM, Klaij, Christiaan
c.kl...@marin.nlmailto:c.kl...@marin.nl wrote:
Matt,
Thanks, after applying the fix to my petsc-3.5.1 install, the
small Fortran program works
: Wednesday, September 03, 2014 4:44 PM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] PCFieldSplitSetSchurPre in fortran
On Wed, Sep 3, 2014 at 9:43 AM, Klaij, Christiaan
c.kl...@marin.nlmailto:c.kl...@marin.nl wrote:
I'm sorry, how do I do that?
print it
Matt
Chris
-912.htm
This e-mail may be confidential, privileged and/or protected by copyright. If
you are not the intended recipient, you should return it to the sender
immediately and delete your copy from your system.
From: Klaij, Christiaan
Sent: Friday, August 29
it to the sender
immediately and delete your copy from your system.
From: Matthew Knepley knep...@gmail.com
Sent: Friday, August 29, 2014 4:14 PM
To: Klaij, Christiaan
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] PCFieldSplitSetSchurPre in fortran
MatNestGetISs works as expected in fortran with petsc-3.5.1. Thanks!
Chris
From: Satish Balay ba...@mcs.anl.gov
Sent: Monday, July 07, 2014 6:37 PM
To: Klaij, Christiaan
Cc: petsc-users
Subject: Re: [petsc-users] MatNestGetISs in fortran
My change
Just for future reference, below are the correct comments that
should be at the top of ex70.
(Jed, Matt, I sent these to the maintainer list a long time
ago ([petsc-maint #155997]), but you never put them in, it would
avoid a lot of confusion)
/* Poiseuille flow problem.
Satish,
Thanks for your reply! I've been on holiday for some time, do you still want me
to test this or has it been fixed with the 3.5 release?
Chris
From: Satish Balay ba...@mcs.anl.gov
Sent: Monday, June 16, 2014 5:41 PM
To: Klaij, Christiaan
Cc
317 49 32 45, I www.marin.nl
From: Klaij, Christiaan
Sent: Monday, June 02, 2014 9:54 AM
To: petsc-users@mcs.anl.gov
Subject: RE: MatNestGetISs in fortran
Just a reminder. Could you please add fortran support for MatNestGetISs
www.marin.nl
From: Klaij, Christiaan
Sent: Tuesday, May 27, 2014 3:47 PM
To: petsc-users@mcs.anl.gov
Subject: MatNestGetISs in fortran
I'm trying to use MatNestGetISs in a fortran program but it seems to be missing
from the fortran include file (PETSc 3.4).
I'm trying to use MatNestGetISs in a fortran program but it seems to be missing
from the fortran include file (PETSc 3.4).
dr. ir. Christiaan Klaij
CFD Researcher
Research Development
E mailto:c.kl...@marin.nl
T +31 317 49 33 44
MARIN
2, Haagsteeg, P.O. Box 28, 6700 AA Wageningen, The
I'm having problems using mat_view in petsc 3.4.3 in combination
with a prefix. For example in ../snes/examples/tutorials/ex70:
mpiexec -n 2 ./ex70 -nx 16 -ny 24 -ksp_type fgmres -pc_type fieldsplit
-pc_fieldsplit_type schur -pc_fieldsplit_schur_fact_type lower -user_ksp
-a00_mat_view
Something's wrong with
http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatSetOption.html
In Firefox 24.3.0, the option description is continued at the bottom of the
page, after the examples.
dr. ir. Christiaan Klaij
CFD Researcher
Research Development
E mailto:c.kl...@marin.nl
If I understand correctly, the first step to running SNES ex62.c is to create
the header file, however:
$ $PETSC_DIR/bin/pythonscripts/PetscGenerateFEMQuadrature.py dim order dim 1
laplacian dim order 1 1 gradient $PETSC_DIR/src/snes/examples/tutorials/ex62.h
Traceback (most recent call last):
I can't do that, it's a matshell.
You can use MatComputeExplicitOperator or MatFDColoring to get the entries.
Also, make sure the MATSHELL is really, truly a linear operator. You'll get
very confusing results if you accidentally leave some nonlinearity inside
that function.
Thanks for
What happens if you use -pc_type lu ?
Barry,
I can't do that, it's a matshell.
Chris
On Sep 21, 2012, at 3:29 AM, Klaij, Christiaan C.Klaij at marin.nl wrote:
When I use zero initial guess, GMRES with left PC gives a huge
jump in true resisdual between iteration 0 and 1 and GMRES
I'm solving a system with GMRES using the same preconditioner either on
the left or on the right.
For left preconditioning I get two orders of reduction for the
preconditioned residual in 20 its:
Notice here that your preconditioner is far from one. It manages to
blow
When I use zero initial guess, GMRES with left PC gives a huge
jump in true resisdual between iteration 0 and 1 and GMRES with
right PC is stuck, the solution remains zero, as mentioned before.
When I use the Knoll trick, both issues are gone (!) and I do get
similar results for left and right
I'm solving a system with GMRES using the same preconditioner either on
the left or on the right.
For left preconditioning I get two orders of reduction for the
preconditioned residual in 20 its:
Notice here that your preconditioner is far from one. It manages to blow
up the true
Is the preconditioner singular? How about the matrix? Is
the preconditioner a linear operator? Did you try KSPFGMRES
on it?
The matrix comes from incompressible Navier-Stokes and the
preconditioner is of Schur complement type (SIMPLE). Both are
shells. FGMRES (and GCR) give the same
1 - 100 of 167 matches
Mail list logo