On Oct 13, 2022, at 11:23 PM, Matthew Knepley <knepley@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@outlook.com> wrote:

Hi Matt,

Sure, please see the attached configure.log file.  Thanks!

Best regards,

Peng Sun

---

From: Matthew Knepley <knepley@gmail.com>
Sent: Thursday, October 13, 2022 6:34 AM
To: Peng Sun <psun@outlook.com>
Cc: petsc-users@mcs.anl.gov <petsc-users@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@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 petsc4py
import 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/

-- 

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/