On 05/23/2017 07:24 PM, Jack wrote:
Originally, I suspected that when solving two linear systems simultaneously by
two threads would reduce time. But now it seems that this idea increases the
complexity of communication between MPI communicators and coding, and also not
surely to decrease compu
On 05/24/2017 06:02 AM, Juan Carlos Araujo Cabarcas wrote:
I managed to do it by looking at: *PETScWrappers::MatrixBase*, but my solution
is extremely ugly!
|
doublemax_coup_dofs =dof_handler.max_couplings_between_dofs();
PETScWrappers::SparseMatrixlocal_M;
local_M.reinit (par.dofs,par.dofs,
I managed to do it by looking at: *PETScWrappers::MatrixBase*, but my
solution is extremely ugly!
double max_coup_dofs = dof_handler.max_couplings_between_dofs();
PETScWrappers::SparseMatrix local_M;
local_M.reinit (par.dofs,par.dofs, max_coup_dofs);
MatSetOption( local_M, MAT_NEW_NONZ
Dear all,
I have a matrix written by using PETSc in a Binary file, and I am trying
work with it in my dealii environment. For this, I read a system matrix
from PETSc:
Mat M_read;
PetscViewerBinaryOpen(PETSC_COMM_WORLD,file.c_str (),FILE_MODE_READ,&fd);
MatCreate(PETSC_COMM_WORLD,
Thanks, Denis and Wolfgang for the quick reply.
I will give it a try and report back whatever issue I encounter.
El miércoles, 24 de mayo de 2017, 7:26:42 (UTC+2), Denis Davydov escribió:
>
>
> On Tuesday, May 23, 2017 at 5:22:44 PM UTC+2, Juan Carlos Araujo Cabarcas
> wrote:
>>
>> Dear all,
>>
