[petsc-users] Issue with single precision complex numbers in petsc4py

[Matthew Knepley]

On Thu, Oct 13, 2022 at 5:27 PM Peng Sun <psun at outlook.com> wrote:

> Hi Stefano,
>
> Sure, please see the following.  The PETSC_ARCH field is empty in the
> printout despite the fact that it was set to 'arch-linux-c-opt' in the
> shell.
>
> {'PETSC_DIR': '/home/pesun/.emopt', 'PETSC_ARCH': ''}
>

Can you show the whole output?

Also, did you remember to 'export' it so that it goes to subshells?

   Matt


> Best regards,
> Peng Sun
> ------------------------------
> *From:* Stefano Zampini <stefano.zampini at gmail.com>
> *Sent:* Thursday, October 13, 2022 1:57 PM
> *To:* Matthew Knepley <knepley at gmail.com>
> *Cc:* Peng Sun <psun at outlook.com>; petsc-users at mcs.anl.gov <
> petsc-users at mcs.anl.gov>
> *Subject:* Re: [petsc-users] Issue with single precision complex numbers
> in petsc4py
>
> Matt
>
> Yes, petsc4py does the right thing. This is probably. Picking up the wrong
> PETSc arch.
>
> Peng, can you please run this?
>
> import petsc4py
> petsc4py.init()
> print(petsc4py.get_config())
>
> On Oct 13, 2022, at 11:23 PM, Matthew Knepley <knepley at gmail.com> wrote:
>
> Lisandro,
>
> PETSc is compiled for single. Does petsc4py respect this, or does it
> always use double for getArray() and friends?
>
>   Thanks,
>
>      Matt
>
> On Thu, Oct 13, 2022 at 11:42 AM Peng Sun <psun at outlook.com> wrote:
>
> Hi Matt,
>
> Sure, please see the attached configure.log file.  Thanks!
>
> Best regards,
> Peng Sun
>
>
> ------------------------------
> *From:* Matthew Knepley <knepley at gmail.com>
> *Sent:* Thursday, October 13, 2022 6:34 AM
> *To:* Peng Sun <psun at outlook.com>
> *Cc:* petsc-users at mcs.anl.gov <petsc-users at mcs.anl.gov>
> *Subject:* Re: [petsc-users] Issue with single precision complex numbers
> in petsc4py
>
> First send configure.log so we can see the setup.
>
>   Thanks,
>
>       Matt
>
> On Thu, Oct 13, 2022 at 12:53 AM Peng Sun <psun at outlook.com> wrote:
>
> Dear PETSc community,
>
>
> I have a question regarding the single precision complex numbers of
> petsc4py.  I configured PETSc with the “--with-scalar-type=complex
> --with-precision=single" option before compiling, but all the DA structures
> I created with petsc4py had double precision.
>
>
> Here is a minimum test code on Python 3.8/PETSc 3.12/petsc4py 3.12: both
> print commands show data type of complex128.  Could anybody please help
> me?  Thanks!
>
>
> import petsc4pyimport sys
> petsc4py.init(sys.argv)from petsc4py import PETSc
>
> da=PETSc.DA().create(sizes=[2,2,2],dof=1,stencil_type=0,stencil_width=1,boundary_type=1)
> da_1 = da.createGlobalVec()print(petsc4py.PETSc.ComplexType)print(da_1.getArray().dtype)
>
>
>
>
> Best regards,
>
> Peng Sun
>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>
> https://www.cse.buffalo.edu/~knepley/
> <http://www.cse.buffalo.edu/~knepley/>
>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>
> https://www.cse.buffalo.edu/~knepley/
> <http://www.cse.buffalo.edu/~knepley/>
>
>
>

-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221013/83c2a932/attachment-0001.html>