<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
Hi Stefano,</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
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.<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted0">
{'PETSC_DIR': '/home/pesun/.emopt', 'PETSC_ARCH': ''}</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted0">
<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted0">
Best regards,</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted0">
Peng Sun<br>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Stefano Zampini <stefano.zampini@gmail.com><br>
<b>Sent:</b> Thursday, October 13, 2022 1:57 PM<br>
<b>To:</b> Matthew Knepley <knepley@gmail.com><br>
<b>Cc:</b> Peng Sun <psun@outlook.com>; petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov><br>
<b>Subject:</b> Re: [petsc-users] Issue with single precision complex numbers in petsc4py</font>
<div> </div>
</div>
<div class="" style="word-wrap:break-word; line-break:after-white-space">Matt
<div class=""><br class="">
</div>
<div class="">Yes, petsc4py does the right thing. This is probably. Picking up the wrong PETSc arch.</div>
<div class=""><br class="">
</div>
<div class="">Peng, can you please run this?</div>
<div class="">
<div class=""><br class="">
</div>
<div class="">import petsc4py<br class="">
<div>petsc4py.init()</div>
<div>print(petsc4py.get_config())</div>
<div><br class="">
<blockquote type="cite" class="">
<div class="">On Oct 13, 2022, at 11:23 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com" class="">knepley@gmail.com</a>> wrote:</div>
<br class="x_Apple-interchange-newline">
<div class="">
<div dir="ltr" class="" style="font-family:Helvetica; font-size:12px; font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
Lisandro,
<div class=""><br class="">
</div>
<div class="">PETSc is compiled for single. Does petsc4py respect this, or does it always use double for getArray() and friends?</div>
<div class=""><br class="">
</div>
<div class=""> Thanks,</div>
<div class=""><br class="">
</div>
<div class=""> Matt</div>
</div>
<br class="" style="font-family:Helvetica; font-size:12px; font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
<div class="x_gmail_quote" style="font-family:Helvetica; font-size:12px; font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
<div dir="ltr" class="x_gmail_attr">On Thu, Oct 13, 2022 at 11:42 AM Peng Sun <<a href="mailto:psun@outlook.com" class="">psun@outlook.com</a>> wrote:<br class="">
</div>
<blockquote class="x_gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left-width:1px; border-left-style:solid; border-left-color:rgb(204,204,204); padding-left:1ex">
<div class="x_msg-9108617864171965949">
<div dir="ltr" class="">
<div class="" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt">Hi Matt,</div>
<div class="" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt"><br class="">
</div>
<div class="" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt">Sure, please see the attached configure.log file. Thanks!</div>
<div class="" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt"><br class="">
</div>
<div class="" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt">Best regards,</div>
<div class="" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt">Peng Sun<br class="">
</div>
<div class="" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt"><br class="">
</div>
<div class="" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt"><br class="">
</div>
<div id="x_m_-7847157880467528269appendonsend" class=""></div>
<hr class="" style="display:inline-block; width:965.015625px">
<div id="x_m_-7847157880467528269divRplyFwdMsg" dir="ltr" class=""><font face="Calibri, sans-serif" class="" style="font-size:11pt"><b class="">From:</b><span class="x_Apple-converted-space"> </span>Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank" class="">knepley@gmail.com</a>><br class="">
<b class="">Sent:</b><span class="x_Apple-converted-space"> </span>Thursday, October 13, 2022 6:34 AM<br class="">
<b class="">To:</b><span class="x_Apple-converted-space"> </span>Peng Sun <<a href="mailto:psun@outlook.com" target="_blank" class="">psun@outlook.com</a>><br class="">
<b class="">Cc:</b><span class="x_Apple-converted-space"> </span><a href="mailto:petsc-users@mcs.anl.gov" target="_blank" class="">petsc-users@mcs.anl.gov</a><span class="x_Apple-converted-space"> </span><<a href="mailto:petsc-users@mcs.anl.gov" target="_blank" class="">petsc-users@mcs.anl.gov</a>><br class="">
<b class="">Subject:</b><span class="x_Apple-converted-space"> </span>Re: [petsc-users] Issue with single precision complex numbers in petsc4py</font>
<div class=""> </div>
</div>
<div class="">
<div dir="ltr" class="">First send configure.log so we can see the setup.
<div class=""><br class="">
</div>
<div class=""> <span class="x_Apple-converted-space"> </span>Thanks,</div>
<div class=""><br class="">
</div>
<div class=""> <span class="x_Apple-converted-space"> </span>Matt</div>
</div>
<br class="">
<div class="">
<div dir="ltr" class="">On Thu, Oct 13, 2022 at 12:53 AM Peng Sun <<a href="mailto:psun@outlook.com" target="_blank" class="">psun@outlook.com</a>> wrote:<br class="">
</div>
<blockquote class="" style="margin:0px 0px 0px 0.8ex; border-left-width:1px; border-left-style:solid; border-left-color:rgb(204,204,204); padding-left:1ex">
<div class="">
<div dir="ltr" class="">
<div class="" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt">
<p class="" style="margin-right:0in; margin-left:0in; font-size:11pt; font-family:Calibri,sans-serif">
Dear PETSc community,<u class=""></u> <u class=""></u></p>
<p class="" style="margin-right:0in; margin-left:0in; font-size:11pt; font-family:Calibri,sans-serif">
<br class="">
</p>
<p class="" style="margin-right:0in; margin-left:0in; font-size:11pt; font-family:Calibri,sans-serif">
I have a question regarding the single precision complex numbers of petsc4py.<span class=""> <span class="x_Apple-converted-space"> </span></span>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. </p>
<p class="" style="margin-right:0in; margin-left:0in; font-size:11pt; font-family:Calibri,sans-serif">
<br class="">
</p>
<p class="" style="margin-right:0in; margin-left:0in; font-size:11pt; font-family:Calibri,sans-serif">
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!</p>
<p class="" style="margin-right:0in; margin-left:0in; font-size:11pt; font-family:Calibri,sans-serif">
<br class="">
</p>
<div class="" style="margin:0px 0in; font-size:11pt; font-family:Calibri,sans-serif">
<pre class="" style="margin-top:0px; width:auto; max-height:600px; overflow:auto; text-align:left"><code class="" style="margin:0px"><span class="" style="margin:0px; font-size:13px">import</span> petsc4py
<span class="" style="margin:0px; font-size:13px">import</span> sys
petsc4py.init(sys.argv)
<span class="" style="margin:0px; font-size:13px">from</span> petsc4py <span class="" style="margin:0px; font-size:13px">import</span> PETSc
da=PETSc.DA().create(sizes=[<span class="" style="margin:0px; font-size:13px">2</span>,<span class="" style="margin:0px; font-size:13px">2</span>,<span class="" style="margin:0px; font-size:13px">2</span>],dof=<span class="" style="margin:0px; font-size:13px">1</span>,stencil_type=<span class="" style="margin:0px; font-size:13px">0</span>,stencil_width=<span class="" style="margin:0px; font-size:13px">1</span>,boundary_type=<span class="" style="margin:0px; font-size:13px">1</span>)
da_1 = da.createGlobalVec()
<span class="" style="margin:0px; font-size:13px">print</span>(petsc4py.PETSc.ComplexType)
<span class="" style="margin:0px; font-size:13px">print</span>(da_1.getArray().dtype)</code></pre>
<br class="">
</div>
<p class="" style="margin-right:0in; margin-left:0in; font-size:11pt; font-family:Calibri,sans-serif">
<br class="">
</p>
<p class="" style="margin-right:0in; margin-left:0in; font-size:11pt; font-family:Calibri,sans-serif">
<br class="">
</p>
<p class="" style="margin-right:0in; margin-left:0in; font-size:11pt; font-family:Calibri,sans-serif">
Best regards,</p>
<p class="" style="margin-right:0in; margin-left:0in; font-size:11pt; font-family:Calibri,sans-serif">
Peng Sun</p>
</div>
</div>
</div>
</blockquote>
</div>
<br clear="all" class="">
<div class=""><br class="">
</div>
--<span class="x_Apple-converted-space"> </span><br class="">
<div dir="ltr" class="">
<div dir="ltr" class="">
<div class="">
<div dir="ltr" class="">
<div class="">
<div dir="ltr" class="">
<div class="">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br class="">
-- Norbert Wiener</div>
<div class=""><br class="">
</div>
<div class=""><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank" class="">https://www.cse.buffalo.edu/~knepley/</a><br class="">
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br clear="all" class="" style="font-family:Helvetica; font-size:12px; font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
<div class="" style="font-family:Helvetica; font-size:12px; font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
<br class="">
</div>
<span class="" style="font-family:Helvetica; font-size:12px; font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none; float:none; display:inline!important">--<span class="x_Apple-converted-space"> </span></span><br class="" style="font-family:Helvetica; font-size:12px; font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
<div dir="ltr" class="x_gmail_signature" style="font-family:Helvetica; font-size:12px; font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
<div dir="ltr" class="">
<div class="">
<div dir="ltr" class="">
<div class="">
<div dir="ltr" class="">
<div class="">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br class="">
-- Norbert Wiener</div>
<div class=""><br class="">
</div>
<div class=""><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank" class="">https://www.cse.buffalo.edu/~knepley/</a></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</body>
</html>