Yes but in the documentation it should be clear that if the global
indices are not contiguous in each process, even if the memory is
locally contiguous in each process, then there will be communication to
build the matrix and the vector.
Cheers,
Enrico
On 19/10/2023 20:41, Matthew Knepley
Sorry but I don't want another partition, Petsc internally is changing
the partition. I would like to have the same partition that I have in
the application. Is the example not clear?
On 19/10/2023 19:43, Matthew Knepley wrote:
On Thu, Oct 19, 2023 at 1:00 PM Enrico <mailto:degreg...
) = ', xx_v(i)
end do
CALL VecRestoreArrayReadF90(b, xx_v, ierr)
deallocate(glb_index, array)
CALL VecDestroy(b,ierr)
CALL PetscFinalize(ierr)
end program main
On 19/10/2023 17:36, Matthew Knepley wrote:
On Thu, Oct 19, 2023 at 11:33 AM Enrico <mailto:degreg...@dkrz.de>&
have {2,3,4} on process 0 also in Petsc. Is it more clear now?
On 19/10/2023 17:36, Matthew Knepley wrote:
On Thu, Oct 19, 2023 at 11:33 AM Enrico <mailto:degreg...@dkrz.de>> wrote:
The layout is not poor, just the global indices are not contiguous,this
has nothing to do with
and there is no communication when I fill or
fetch the vectors and the matrix.
Cheers,
Enrico
On 19/10/2023 17:21, Matthew Knepley wrote:
On Thu, Oct 19, 2023 at 10:51 AM Enrico <mailto:degreg...@dkrz.de>> wrote:
In the application the storage is contiguous but the global indexing
etch the values
of the application ordering?
Cheers,
Enrico
On 19/10/2023 15:25, Matthew Knepley wrote:
On Thu, Oct 19, 2023 at 8:57 AM Enrico <mailto:degreg...@dkrz.de>> wrote:
Maybe I wasn't clear enough. I would like to completely get rid of
Petsc
orderi
19, 2023 at 6:51 AM Enrico <mailto:degreg...@dkrz.de>> wrote:
Hello,
if I create an application ordering using AOCreateBasic, should I
provide the same array for const PetscInt myapp[] and const PetscInt
mypetsc[] in order to get the same ordering of the application
with
defined in PETSC as in my application, how can I use it to create the
actual vector and matrix?
Thanks in advance for the help.
Cheers,
Enrico
On 18/10/2023 13:39, Matthew Knepley wrote:
On Wed, Oct 18, 2023 at 5:55 AM Enrico <mailto:degreg...@dkrz.de>> wrote:
Hello,
I'm
not a big issue, but after solving the linear
system I get the solution vector x and I want to get the values in the
correct processes. Is there a way to get vector values from other
processes or to get a mapping so that I can do it myself?
Cheers,
Enrico Degregori
lls.h:
CC-1043 craycc: ERROR File = /usr/include/bits/cmathcalls.h, Line = 55
_Complex can only be used with floating-point types.
__MATHCALL (cacos, (_Mdouble_complex_ __z));
^
Am I doing something wrong?
Regards,
Enrico Degregori
th log and
debugging set to 0 so this should be fine.
Do you know about a branch that support this kind of issue or a way to
solve it?
Enrico
11 matches
Mail list logo