Sorry, I meant: what is the reason for MatAssemblyBegin/End being run
for matrix A, in the case when A != B?
Best regards,
Torquil Sørensen
On 25/01/16 11:58, Torquil Macdonald Sørensen wrote:
> Hi!
>
> I have been looking at some of the PETSc examples where TSSetIJacobian,
> and
Hi!
Is it intentional that Petsc prints duplicates of the matrix-related
options in the following test program that creates two matrices?:
-
Mat A;
ierr = MatCreate(PETSC_COMM_SELF, A); CHKERRQ(ierr);
ierr = MatSetType(A, MATSEQAIJ); CHKERRQ(ierr);
Mat B;
ierr =
,PETSCSPRNG)
Barry
On Jun 30, 2015, at 5:42 PM, Torquil Macdonald Sørensen torq...@gmail.com
wrote:
Thanks. So what can I do to create a PetscRandom object that is not
affected by the SLEPc-functions that I'm calling?
Best regards,
Torquil Sørensen
On 30/06/15 23:43, Jose E. Roman wrote
Hi!
I'm experiencing some problems using PetscRandom and SVD from SLEPc. For
some reason, the order of two seemingly unrelated parts of the program
below affects the random values in a vector. The test program below, run
with mpiexec -n 2, will print a vector with two identical component
values.
values in each process.
Jose
El 30/6/2015, a las 23:05, Torquil Macdonald Sørensen torq...@gmail.com
escribió:
Hi!
I'm experiencing some problems using PetscRandom and SVD from SLEPc. For
some reason, the order of two seemingly unrelated parts of the program
below affects the random values
Hi!
Using Petsc 3.4.4, compiled to support complex numbers, I'm getting a
completely white figure when plotting a matrix with purely imaginary
elements, using MatView(mat, PETSC_VIEWER_DRAW_WORLD). So I cannot see the
nonzero structure from the plot.
For real matrices (still using the
On 01/04/14 16:41, Jed Brown wrote:
Torquil Macdonald Sørensen torq...@gmail.com writes:
Hi!
Using Petsc 3.4.4, compiled to support complex numbers, I'm getting a
completely white figure when plotting a matrix with purely imaginary
elements, using MatView(mat, PETSC_VIEWER_DRAW_WORLD). So I
On 19/12/13 00:50, Jed Brown wrote:
I'm using PCMG. My minimal testcase program has everything
hardcoded. I'm not providing any runtime options other than
-start_in_debugger, so everything can be seen from the code. I've
since made a couple of small changes that might influence line number
Hi!
In the documentation of MatCreateMPIAdj it says that the fifth argument
j should be sorted for each row:
http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatCreateMPIAdj.html
The same page links to an example:
How well is PETSC suited for solving elliptic PDEs on four-dimensional
domains? Looking at the list of commands, I see e.g. DMDACreate{1d,2d,3d}
with no 4d version. On the other hand, I see you have
VecGetArray{1d,2d,3d,4d}.
Is there e.g. a convenient way of using the PETSCs data structures for
On 30 January 2014 22:05, Jed Brown j...@jedbrown.org wrote:
I did not have an opportunity to test what happens if I don't set any
restriction, but I'm hoping it will then use P^T, as before.
Yes, it will.
Only setting the interpolation, and not the restriction, results in a
crash for me
On 30/01/14 17:45, Jed Brown wrote:
Boris Kaus k...@uni-mainz.de writes:
While doing multigrid it is possible to explicitly set restriction and
prolongation/interpolation operators in PETSC using PCMGSetRestriction
and PCMGSetInterpolation. In many cases, the restriction (R) and
On 30/01/14 22:05, Jed Brown wrote:
Torquil Macdonald Sørensen torq...@gmail.com writes:
This question was similar to my recent unanswered one:
http://lists.mcs.anl.gov/pipermail/petsc-users/2014-January/020382.html
That's why I Cc'd you.
Thanks, I didn't notice that. I am actually very
Hi!
I'm using todays Petsc git master on Linux and tried:
./prog -mg_coarse_ksp_plot_eigenvalues -draw_pause -1
but the plot window simply appears and dissappears at once. On the other
hand, if I do
./prog -mg_levels_1_est_ksp_gmres_krylov_monitor -draw_pause -1
it works as expected, i.e. I
The window appeared and dissapeared, so I just assumed that something was
plotted. But probably not, for the reason you mention.
Thanks!
Best regards
Torquil Sørensen
On 27 January 2014 16:57, Barry Smith bsm...@mcs.anl.gov wrote:
On Jan 27, 2014, at 5:59 AM, Torquil Macdonald Sørensen
On 27/01/14 23:27, Danyang Su wrote:
Does the following convergence monitoring make sense? I thought the
solver should stop before iteration 36, but it stops at iteration 70.
0 KSP preconditioned resid norm 4.208619583781e-002 true resid norm
2.227983625699e-005 ||r(i)||/||b||
Hi!
I'm experimenting with the PCMG preconditioner with two levels. I have
created a system matrix A, an interpolation matrix P and a restriction
matrix R which differs from P^T. In addition, I have included
PCMGSetGalerkin(pc, PETSC_TRUE).
When I extract the system matrix A, interpolation
PCSetUp_MG that called DMCoarsen, which then crashed. So I guess I'm
wondering how to tell PETSc that I'm not using any DM structure.
Best regards
Torquil Sørensen
On 14 December 2013 14:46, Torquil Macdonald Sørensen torq...@gmail.comwrote:
Hi!
I'm getting segfault problems when trying
On 18/12/13 21:19, Jed Brown wrote:
In short, I'm wondering if it is possible to use TSSolve without using
any DM structure:
When I was debugging my program I tried to analyse the difference between
the TSSolve-based time iteration program, and a version where I do the
time-stepping manually
Hi!
I'm getting segfault problems when trying to use TSStep or TSSolve with
TSBEULER and TSLINEAR, together with PCMG. So I have reduced my code to
a minimal (and mathematically trivial) test case which displays the same
error, and I'm hoping someone here will spot any errors that I've made
in my
Hi!
I'm having some problems getting PCMG to work when manually populating the
interpolation matrix, i.e. using MatSetValue to set its values, instead of
the use of DMCreateInterpolation, as I have seen in an example. I have to
use MatSetValue to populate my PETSc matrices, because I'm obtaining
, 2013, at 1:16 PM, Torquil Macdonald Sørensen torq...@gmail.com
wrote:
Hi!
I'm having some problems getting PCMG to work when manually populating the
interpolation matrix, i.e. using MatSetValue to set its values, instead of
the use of DMCreateInterpolation, as I have seen in an example. I
On 30/10/13 19:01, Xiaoye S. Li wrote:
The sequential SuperLU has ILUTP implementation, not in parallel
versions. PETSc already supports the option of using SuperLU.
In SuperLU distribution:
EXAMPLE/zitersol.c : an example to use GMRES with ILUTP
preconditioner (returned from driver
On 30/10/13 16:48, Matthew Knepley wrote:
On Wed, Oct 30, 2013 at 8:58 AM, Torquil Macdonald Sørensen
torq...@gmail.com mailto:torq...@gmail.com wrote:
Thanks Matthew!
The problem I'm working on is the Dirac equation, in various
number of dimensions. I'm going to take a stab
Macdonald Sørensen
torq...@gmail.com wrote:
Hi!
I have a linear problem Ax=b, with complex values, that is solved very
well using the ILUTP and GMRES implementations in GMM++. None of the other
preconditioners available in GMM++ would work (e.g. ILU and ILUT) for this
problem.
Then I
Hi!
I have a linear problem Ax=b, with complex values, that is solved very well
using the ILUTP and GMRES implementations in GMM++. None of the other
preconditioners available in GMM++ would work (e.g. ILU and ILUT) for this
problem.
Then I tried the same problem using ILU and GMRES in PETSc,
Hi!
I'm using PETSc 3.4.3, with only default runtime settings, to solve Ax=b,
for complex numbers (I have used --with-scalar-type=complex).
After reading sections 4.3.1, 4.3.2, 4.3.3 in the manual, I decided to
investigate the convergence in my program, since I'm having some problems
when the
On 25/10/13 18:59, Matthew Knepley wrote:
On Fri, Oct 25, 2013 at 11:33 AM, Torquil Macdonald Sørensen
torq...@gmail.com mailto:torq...@gmail.com wrote:
Hi!
I'm using PETSc 3.4.3, with only default runtime settings, to
solve Ax=b, for complex numbers (I have used
--with-scalar
/ex1.c
FAQ #6 mentions MUMPS, which I don't have the moment, so I'll have to look
into that.
- Torquil
On 18 October 2013 16:46, Jose E. Roman jro...@dsic.upv.es wrote:
El 18/10/2013, a las 16:42, Torquil Macdonald Sørensen escribió:
Hi!
I'm just starting to use SLEPc on a quad-core
you build PETSc --with-debugging=0 ?
You should send performance data obtained with -log_summary, otherwise we
cannot guess what happens.
Jose
El 18/10/2013, a las 17:21, Torquil Macdonald Sørensen escribió:
It is a standard eigenvalue problem, and I'm only using default algorithm
options, so
30 matches
Mail list logo