Dear petsc-users,
I am a beginner in petsc and I have some question with the python interface. I
am trying to solve problem with mumps (or other direct sparse solver)
I have written the following piece of code
ksp = PETSc.KSP()
ksp.create(PETSc.COMM_WORLD)
ksp.setOperators(A)
Benoit :
mumps options can be called from PETSc C-code as
ierr = PetscOptionsInsertString(-mat_mumps_icntl_13 1 -mat_mumps_icntl_24
1 -mat_mumps_cntl_3 1e-12);CHKERRQ(ierr);
or see petsc/src/ksp/ksp/examples/tutorials/ex52.c
I do not know how to translate these calls into python API though.
On 15 April 2015 at 18:15, Hong hzh...@mcs.anl.gov wrote:
Benoit :
mumps options can be called from PETSc C-code as
ierr = PetscOptionsInsertString(-mat_mumps_icntl_13 1 -mat_mumps_icntl_24 1
-mat_mumps_cntl_3 1e-12);CHKERRQ(ierr);
or see petsc/src/ksp/ksp/examples/tutorials/ex52.c
I do
Hong and Emil,
You will need to answer this.
On Apr 15, 2015, at 6:24 AM, Dave Makhija makhi...@colorado.edu wrote:
Barry,
Is the method exact if I have a time dependent mass matrix (dF/dU_t changes
each time step)? I'm not sure it's possible if TSSetIJacobian expects a
On 4/15/15 12:37 PM, Barry Smith wrote:
Hong and Emil,
You will need to answer this.
On Apr 15, 2015, at 6:24 AM, Dave Makhija makhi...@colorado.edu wrote:
Barry,
Is the method exact if I have a time dependent mass matrix (dF/dU_t changes
each time step)? I'm not sure it's
On 3/24/15 5:31 AM, Lisandro Dalcin wrote:
Emil, is there any chance you can try my code with your problem? I
really need some feedback to push this to PETSc, otherwise
Hi Lisandro - we checked ts_alpha_adapt and we tested it on a small
system (mildly stiff van der Pol ODE). I enclosed a
The more adaptors the better
On Apr 15, 2015, at 8:51 PM, Emil Constantinescu emcon...@mcs.anl.gov wrote:
On 3/24/15 5:31 AM, Lisandro Dalcin wrote:
Emil, is there any chance you can try my code with your problem? I
really need some feedback to push this to PETSc, otherwise
Hi
Barry,
Is the method exact if I have a time dependent mass matrix (dF/dU_t changes
each time step)? I'm not sure it's possible if TSSetIJacobian expects a
function which computes the combination of dF/dU and a*dF/dU_t. But,
perhaps there is something I'm missing.
Otherwise, it looks very