This has happened to me if I forget to call MatAssemblyBegin/End. The entry
will still print out in the MatView(), but it won't compute.
Might not be your problem, but worth a shot.
Original message
Date: Wed, 31 Aug 2011 18:15:14 -0600
From: petsc-users-bounces at mcs.anl.gov (on
Dear all,
I used MatCreate() to create a non-square matrix B, and set values to it
in parallel. Now i want to get the tranpose of it to solve Ax=B', where A
is SeqDense.
I define a SeqDense C=B' and use MatTranspose() to get it, then apply
MatMatSolve(fact,C,x), but i get the following error:
the whole stack trace.
It looks like C=B' is not being partitioned correctly.
Did you take my earlier suggestions about how to partition this problem?
-- next part --
An HTML attachment was scrubbed...
URL:
http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110901
Thank you for your reply~
I tried to use
PetscSplitOwnership(PETSC_COMM_WORLD,n,N);
MatCreateSeqDense(PETSC_COMM_SELF,M,n,PETSC_NULL,B);
But it seems not working. This command is make B store in N/n processors
with n columns in each processor? I try to calculate the N/n parts in
parallel, but
://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110901/de2b5e5b/attachment.htm
Could be random memory corruption. perhaps you can valgrind the code
on a small problem [on smaller set of nodes] - with --download-mpich
build of petsc.
Satish
On Thu, 1 Sep 2011, Rongliang Chen wrote:
Hi Jed and Andrew,
Thanks for your reply.
It's a little strange that this problem
Dear Hong,
Yes, I am using this SAME_NONZERO_PATTERN to optimize the code,hehe!
Thanks very much!
Best,
Rebecca
On Sep 1, 2011, at 7:35 AM, Hong Zhang wrote:
Rebecca:
Thanks for your kind reply. So if I set up the following:
ierr = DAGetMatrix(DMMGGetDA(dmmg), MATAIJ,
Thank you for your patience.
I still have some confusions. When computing B before MPI_Scan(), could i
compute the values in parallel? After using MPI_Scan(), does that mean the
columns of B will be gathered in one processor?
Here are main steps i took based on your suggestion,
attachment was scrubbed...
URL:
http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110901/ee8df973/attachment.htm
Hello!
I am getting an error with a different number, but the message is the same:
59]PETSC ERROR:
[59]PETSC ERROR: MatFactorNumeric_MUMPS() line 517 in
src/mat/impls/aij/mpi/mumps/mumps.c
[10]PETSC ERROR:
Gabriele:
This time,
you get INFO(1)=-1, An error occurred on processor INFO(2)=39,
according to mumps manual. Please ask mumps developer about
INFO(1)=-1.
My guess is that your matrix might be numerical singular, causing
numerical factorization fail.
Hong
I am getting an error with a
Thank you very much.
With the command:
for(j=0; jn; j++)
{
for(i=0; iM; i++)
{
MatSetValues(B, 1, i, 1, j, Value, INSERT_VALUES);
}
}
This only defines the columns from 0 to n-1? How exactly MPI_Scan()
getting used here?
And when i try to print out the result with
][]
are ghost grids right next to grid of u[i = 16][]).
Thank you so much.
best,
Alan
-- next part --
An HTML attachment was scrubbed...
URL:
http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110901/adbeecac/attachment.htm
Rebecca:
Thanks for your kind reply. So if I set up the following:
? ? ? ?ierr = DAGetMatrix(DMMGGetDA(dmmg), MATAIJ, jacobian);CHKERRQ(ierr);
? ? ? ?KSP ? ? ksp;
? ? ? ?ierr = SNESGetKSP(DMMGGetSNES(dmmg),ksp);CHKERRQ(ierr);
? ? ? ?ierr =
14 matches
Mail list logo