Hong,
Thanks for trying out.
The extra printings are not properly guarded by the print level. I will
fix that. I will look into the hang problem soon.
Sherry
On Mon, Jul 27, 2015 at 7:50 PM, Hong wrote:
> Sherry,
>
> I can repeat hang using petsc/src/ksp/ksp/examples/tutorials/ex10.c:
> m
Sherry,
I can repeat hang using petsc/src/ksp/ksp/examples/tutorials/ex10.c:
mpiexec -n 4 ./ex10 -f0 /homes/hzhang/tmp/Amat_binary.m -rhs 0 -pc_type lu
-pc_factor_mat_solver_package superlu_dist -mat_superlu_dist_parsymbfact
...
.. Starting with 1 OpenMP threads
[0] .. BIG U size 1342464
[0] .. BI
Mark Adams writes:
> I have a 4x4 block matrix, in one big AIJ matrix. I use (global) ISs to
> address individual fields. I have a mass matrix for one field. I would
> like to add a scaled mass matrix to one field at a time.
>
> MatGetLocalSubMatrix looks like what I want but it takes local in
You need to explicitly put a zero entry at all diagonal locations that have
zero.
Barry
1) this is too keep the logic of the symbolic and numeric factorization easier
2) but we don't want to force all matrices to have locations on the diagonal so
we don't automatically insert zero values
Hi all,
I am solving a saddle point problem using LU, and the problem has some
zeros at diagonal. I got error messages as follows:
[0]PETSC ERROR: - Error Message
--
[0]PETSC ERROR: Object is in wrong state
[0]PETSC E
Mat's only have a LocalToGlobal mapping if it is specifically set (for
example DMCreateMatrix() usually sets one for you or you can use
MatSetLocalToGlobalMapping()) (read the manual page) otherwise the Mat simply
does not have the mapping. All this means is one cannot use MatSetValuesLocal()
"Klinvex, Alicia Marie" writes:
> ISLocalToGlobalMapping PETScRowOwnership, PETScColOwnership;
>
> ierr = MatGetLocalToGlobalMapping(A, &PETScRowOwnership,
> &PETScColOwnership); CHKERRV(ierr);
It appears that the local-to-global mapping has not been set. Note that
these mappings do not sign
Hello,
I am trying to determine how the rows and columns of my PETSc Mat are being
mapped to different MPI processes. I call MatGetLocalToGlobalMapping, but the
ISLocalToGlobalMapping obtained for both the rows and columns is null, and I
receive the following error message when I call ISLocal
On Sun, Jul 26, 2015 at 8:48 AM, mahir.ulker-kaust...@tyrens.se <
mahir.ulker-kaust...@tyrens.se> wrote:
> Thank you very much for the input. I will test this approach. I have
> some reading (about the numerical techniques that you have implemented) to
> do as well, so it might take some time for
Hi Hong,
No that is not the correct matrix. Note that I forgot to mention that it is
a complex matrix. I tried loading the matrix I sent you this morning with:
!...Load a Matrix in Binary Format
call
PetscViewerBinaryOpen(PETSC_COMM_WORLD,"Amat_binary.m",FILE_MODE_READ,viewer,ierr)
ca
Zin:
See
petsc/src/mat/examples/tests/ex130.c
petsc/src/ksp/ksp/examples/tutorials/ex52.c
Hong
Hi
>
> I would like to know how I can retrieve the lower triangular matrix
> (possibly permuted or preferably the inverse of L) and the upper triangular
> matrix (preferably the inverse of U) from the L
Hi
I would like to know how I can retrieve the lower triangular matrix
(possibly permuted or preferably the inverse of L) and the upper triangular
matrix (preferably the inverse of U) from the LU factorization done by an
external package (such as MUMPS).
I need them to build a preconditioner of a
12 matches
Mail list logo
