Hello,
I have a time-dependent model that I solve using TSSolve.
And I am trying to adaptively change the step size (dt).
I found that there are some TSAdapt schemes already available.
I have tried TSADAPTBASIC and TSADAPTCFL.
The former runs without any problems, whereas the latter yields th
you.
Best,
Zakariae
From: Barry Smith
Sent: Wednesday, April 12, 2023 3:32:11 PM
To: Jorti, Zakariae
Cc: petsc-users@mcs.anl.gov
Subject: [EXTERNAL] Re: [petsc-users] Question about -memory_view
On Apr 12, 2023, at 4:08 PM, Jorti, Zakariae via petsc-users
Hello,
I am running some matrix computations on 64 nodes, using 640 MPI tasks.
And I wanted to check the memory usage with the -memory_view flag.
I get the following output:
Summary of Memory Usage in PETSc
Maximum (over computational time) process memory:total 3.1056e+11 max
5.9918e
Hello,
I have a code that handles a PETSc Vec on many procs and I would like to use
VecScatterCreateToZero to have all elements of this vector on a single proc,
call VecGetArrayRead on this proc to get the corresponding array to carry out
some diagnostics.
Unfortunately, what I noticed is tha
bly misleading
2) It would be nice to have a view that displayed the nested fieldsplit
preconditioners more clearly
On Jun 9, 2022, at 5:19 PM, Jorti, Zakariae via petsc-users
mailto:petsc-users@mcs.anl.gov>> wrote:
Hi,
I am solving non-linear problem that has 5 unknowns {ni, T, E,
Hi,
I am solving non-linear problem that has 5 unknowns {ni, T, E, B, V}, and for
the preconditioning part, I am using a FieldSplit preconditioner. At the last
fieldsplit/level, we are left with a {B,V} block that tried to precondition in
2 different ways:
a) SuperLU:
-fieldsplit_TEBV_fieldspli
ng you obtained with MatColoringCreate() etc. The same command line
arguments as above allow you to control the coloring algorithm used and to view
them etc.
Barry
On May 10, 2022, at 11:40 AM, Jorti, Zakariae via petsc-users
mailto:petsc-users@mcs.anl.gov>> wrote:
Hi,
I am s
Hi,
I am solving a non-linear problem and using a finite difference approximation
with coloring to compute the Jacobian matrix.
There are several coloring algorithms available in PETSc as indicated here:
https://petsc.org/release/docs/manualpages/Mat/MatColoring.html
And I was wondering how to
Hello,
I am using a TS to solve a differential algebraic equation (DAE).
I do not provide the Jacobian matrix but instead set the TS to use a finite
difference jacobian with coloring.
For debugging, I only solve this DAE on one time step.
Is there a way to output this finite difference jacobian i
use PETScSection. Sorry about the delay!
You can only use PCFieldSplitSetDetectSaddlePoint when your diagonal entries
being zero or non-zero defines the splits correctly.
Am 17.04.2021 um 21:09 schrieb Matthew Knepley
mailto:knep...@gmail.com>>:
On Fri, Apr 16, 2021 at 8:39 PM Jorti, Zaka
l entries
being zero or non-zero defines the splits correctly.
Am 17.04.2021 um 21:09 schrieb Matthew Knepley
mailto:knep...@gmail.com>>:
On Fri, Apr 16, 2021 at 8:39 PM Jorti, Zakariae via petsc-users
mailto:petsc-users@mcs.anl.gov>> wrote:
Hello,
I have a DMStag grid with one dof on
,
Matt
Thanks,
Qi
On Oct 15, 2021, at 3:07 PM, Jorti, Zakariae via petsc-users
mailto:petsc-users@mcs.anl.gov>> wrote:
Hello,
Does the Jacobian approximation using coloring and finite differencing of the
function evaluation work in DMStag?
Thank you.
Best regards,
Zakariae
--
W
Hello,
Does the Jacobian approximation using coloring and finite differencing of the
function evaluation work in DMStag?
Thank you.
Best regards,
Zakariae
sure it
does not happen.
Thanks.
Zakariae
From: Jed Brown
Sent: Monday, August 16, 2021 10:40:20 PM
To: Jorti, Zakariae; petsc-users@mcs.anl.gov
Cc: Tang, Xianzhu
Subject: [EXTERNAL] Re: [petsc-users] malloc error
"Jorti, Zakariae via petsc-users&quo
Hello,
I am using TSSolve to solve a linear problem.
In the FormIJacobian function that I provide to TSSetIJacobian, I first set the
coefficients of both J and Jpre matrices the same way (J and Jpre matrices are
equal in the first step). Then I call MatAXPY to prepare Jpre (Jpre := Jpre -
ano
Hello,
I have a matrix A = [A00 , A01 ; A10, A11].
I extract the submatrix A11 with MatGetSubMatrix.
I only know the global IS is1 and is2, so to get A11 I call:
MatGetSubMatrix(A,is2,is2,MAT_INITIAL_MATRIX,&A11);
I want to modify A11 and update the changes on the global matrix A but I could
rom: Matthew Knepley
Sent: Wednesday, July 7, 2021 12:11:10 PM
To: Jorti, Zakariae
Cc: petsc-users@mcs.anl.gov; Tang, Qi; Tang, Xianzhu
Subject: [EXTERNAL] Re: [petsc-users] Problem with PCFIELDSPLIT
On Wed, Jul 7, 2021 at 1:51 PM Jorti, Zakariae via petsc-users
mailto:petsc-users@mcs.anl.gov&
Hi,
I am trying to build a PCFIELDSPLIT preconditioner for a matrix
J = [A00 A01]
[A10 A11]
that has the following shape:
M_{user}^{-1} = [I -ksp(A00) A01] [ksp(A00) 0] [I
0]
[0I] [0
SaddlePoint when your diagonal entries
being zero or non-zero defines the splits correctly.
Am 17.04.2021 um 21:09 schrieb Matthew Knepley
mailto:knep...@gmail.com>>:
On Fri, Apr 16, 2021 at 8:39 PM Jorti, Zakariae via petsc-users
mailto:petsc-users@mcs.anl.gov>> wrote:
Hello,
I hav
will work well on
the Schur complement.
Barry
On Apr 20, 2021, at 5:26 PM, Matthew Knepley
mailto:knep...@gmail.com>> wrote:
On Tue, Apr 20, 2021 at 4:33 PM Jorti, Zakariae via petsc-users
mailto:petsc-users@mcs.anl.gov>> wrote:
Hello,
I have a block matrix A,
[A_{00} A_{
Hello,
I have a block matrix A,
[A_{00} A_{01}]
[A_{10} A_{11}] ,
and am trying to use PCFIELDSPLIT as a preconditioner on it. The preconditioner
has this shape:
[inv(A_{00}) 0 ]
[ 0-ksp(S)] ;
where S is the Schur Complement matrix:
S = A_{11} - A_{10} *
Hello,
I have a DMStag grid with one dof on each edge and face center.
I want to use a PCFieldSplit preconditioner on a Jacobian matrix that I assume
is already split but I am not sure how to determine the fields.
In the DMStag examples (ex2.c and ex3.c), the function
PCFieldSplitSetDetectSad
Hi Patrick,
Thanks for your responses.
As for the code, I was not granted permission to share it yet. So, I cannot
send it to you for the moment. I apologize for that.
I wanted to let you know that while I was testing my code, I discovered that
when the periodic boundary conditions are activ
Hi,
I implemented a PETSc code to solve Maxwell's equations for the magnetic and
electric fields (B and E) in a cylinder:
0 < r_min <= r <= r_max;with r_max > r_min
phi_min = 0 <= r <= phi_max = 2 π
z_min <= z =< z_max;with z_max > z_min.
I am using a PE
Hi,
I have a PETSc code which works in some cases circumstances whereas in others
it gives me the error message below.
The code works on my machine when I compile PETSc with debug flag on:
./configure PETSC_DIR=/Users/zjorti/software/petsc-3.14.5 PETSC_ARCH=macx
--with-fc=0 --with-mpi-dir=$HO
25 matches
Mail list logo
