Thanks for your explanation so much.
On Tue, Jul 11, 2023 at 12:41 PM Barry Smith wrote:
>
>
> On Jul 11, 2023, at 12:16 PM, neil liu wrote:
>
> I think I understood that.
>
> 1) In order to save memory, use a larger real system
>
>
>This will save a little memory, but requires a slight cha
> On Jul 11, 2023, at 12:16 PM, neil liu wrote:
>
> I think I understood that.
>
> 1) In order to save memory, use a larger real system
This will save a little memory, but requires a slight change to your
model/code for your problem
> 2) vec will store the real part with imaginary part e
I think I understood that.
1) In order to save memory, use a larger real system
2) vec will store the real part with imaginary part equal to zero in the
context of --scalar =complex. But this will not save memory.
Right?
Thanks,
Xiaodong
On Tue, Jul 11, 2023 at 12:03 PM neil liu wrote:
> Do
Do you mean that,
even if I configure with --scalar=complex, when I use vecGetarray to get
the coords, the coords array is still stored in a real manner?
Thanks,
Xiaodong
On Tue, Jul 11, 2023 at 11:54 AM Matthew Knepley wrote:
> On Tue, Jul 11, 2023 at 11:48 AM neil liu wrote:
>
>> Thanks fo
On Tue, Jul 11, 2023 at 11:48 AM neil liu wrote:
> Thanks for your reply.
> I am using DMPLEX to read GMSH. And it seems I have to use
> DMGetCoordinatesLocal();
> VecGetArray(veccords, &coords);//it seems coords has to be defined by
> complex either.
>
> 2) Store the coordinates in coords with i
Thanks for your reply.
I am using DMPLEX to read GMSH. And it seems I have to use
DMGetCoordinatesLocal();
VecGetArray(veccords, &coords);//it seems coords has to be defined by
complex either.
2) Store the coordinates in coords with imaginary part equal to zero.
I got the coordinates from above fu
Not really. PetscScalar (and entries of vectors etc) are either real or
complex, they cannot be mixed.
What parts of PETSc are you using, DM?
You can
1) store the coordinates in arrays and not use Vec for them. (of if you are not
using DM)
2) sore them in Vec but with imaginary part
Dear Petsc developers,
My computation involves both real and complex numbers.
For example, my coordinate is real, but my matrix is complex.
I use VecGetArray to get the coordinates. But it seems coordinates have to
be defined in a complex manner to be consistent with complex matrix.
Is there som
