[petsc-users] a question about PetscSectionCreate

[Dave May]

On Sun, 28 May 2017 at 08:31, leejearl <leejearl at 126.com> wrote:

> Hi, PETSc developer:
>
>      I need to create a PetscSection with a struct. The struct is
> defined as follow,
>
>     typedef struct
>     {
>        PetscReal x;
>        PetscInt id;
>     } testStruct;
>
>     When I run the program, I got a wrong output as follow,
>
>     Vec Object: 1 MPI processes
>    type: seq
> 2.
> 4.94066e-324
> 2.
> 4.94066e-324
> 2.
> 4.94066e-324
> 2.
> 4.94066e-324
> 2.
> 4.94066e-324
> 2.
> 4.94066e-324
> 2.
> 4.94066e-324
> 2.
> 4.94066e-324
>
> But when I defined the struct as
>
>     typedef struct
>     {
>        PetscReal x;
>        PetscReal id;
>     } testStruct;
>
> The output is ok. It seems that  there is some wrong with the memories
> when I define the "id" as a PetscInt type.


Yep.


>
> I can not find out the reasons, and any one can help me with it?


The Vec object can only store quantities of type PetscScalar. It cannot
store PetscInt's and it definitely cannot represent a mixture of
PetscReal's and PetscInt's.


Thanks,
 Dave

The
> source file "test.c" is attached.
>
>
> Thanks,
>
> leejearl
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170528/c1bc4534/attachment.html>