<div dir="ltr"><div><span class="gmail-adx">I just tried to change PetscScalar to the specific datatype PetscComplex in my Fortran code, but the compiler raised the following error:</span></div><div><span class="gmail-adx"></span>PetscComplex   my_coeff</div><div>Error: Unclassifiable statement at (1)</div><div><br></div><div>It seems that the compiler does not know this datatype.<br></div><div><br></div><div>My Petsc was configured with ./configure   --with-scalar-type=complex --with-precision=double <br></div><div><br></div><div>Before the modification, the code was complied successfully. <br></div><div>Could someone help?</div><div><br></div><div>Thanks,</div><div>Youjun<br></div><br><div class="gmail_quote"><div dir="ltr">On Sat, Aug 4, 2018 at 12:28 PM Smith, Barry F. <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
<br>
> On Aug 4, 2018, at 1:11 PM, Youjun Hu <<a href="mailto:youjunhu@gmail.com" target="_blank">youjunhu@gmail.com</a>> wrote:<br>
> <br>
> Thank Barry for the clarification.<br>
> In summary, application code developers should choose the specific datatypes (PetscReal and PetscComplex) rather than the more general datatype PetscScalar<br>
> because application codes usually have fixed datatype and the flexibility provided by PetscScalar is for PETSc library developers and usually not for application developers.<br>
> If an application developer uses PetscScalar to declare a variable that is intended by the developer to be complex and is assigned a complex value later in the code and then used in, e.g.,  MatSetValues, but it happens that petsc is configured with --with-scalar-type=real<br>
> then, it is not possible for this error to be detected at compiling time. If the developer specify PetscComplex, then this inconsistency can be seen by compilers and hopefully an error message can be raised during compiling.<br>
> <br>
> Conclusion: avoid using PetscScalar in application codes unless it is absolutely needed (which is not common). The example application codes provided in tutor folds of petsc often use PetscScalar, which is a little misleading and should better be revised to use concrete datatypes.<br>
<br>
    The reason the examples use PetscScalar is so that all the examples can be used to test both the real and complex versions of the libraries. If they used the specific types they could only be testing one or the other.<br>
<br>
    Note there is no harm in writing applications with PetscScalar it is just usually not needed.<br>
<br>
   Barry<br>
<br>
> Correct me if anybody finds this conclusion is not justified.<br>
> <br>
> Cheers,<br>
> Youjun<br>
> <br>
> <br>
> <br>
> On Sat, Aug 4, 2018 at 11:03 AM Smith, Barry F. <<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>> wrote:<br>
> <br>
>   For a particular configuration (./configure --with-scalar-type= real or complex) the PETSc objects, such as vectors and matrices either contain real values or complex values. In the same running code it is not possible to have some vectors with real numbers and some vectors with complex numbers.<br>
> <br>
>   Depending on how you configure PETSc ./configure --with-scalar-type= real or complex the variable name PetscScalar becomes either PetscReal or PetscComplex. PETSc code is written with PetscScalar so we can reconfigure the same source code and have it work with real numbers or complex.<br>
> <br>
>    For users who just work with real numbers you can simply use PetscReal as your variable type. If working with just with complex numbers you can use PetscComplex. If you want to sometimes configure with real numbers and sometimes configure with complex (which is not common) you use PetscScalar. <br>
> <br>
>    Barry<br>
> <br>
> <br>
> > On Aug 4, 2018, at 11:39 AM, Youjun Hu <<a href="mailto:youjunhu@gmail.com" target="_blank">youjunhu@gmail.com</a>> wrote:<br>
> > <br>
> > Hi all,<br>
> > <br>
> > I am just wondering why PETSc introduces PetscScalar datatype in addition to the PetscReal and PetscComplex.<br>
> > Could someone help clarify the rationale behind this design?<br>
> > I am primarily a Fortran guy. The additional datatype PetscScalar seems to only confuse me rather than provide me any flexibility.<br>
> > Cheers,<br>
> > Youjun Hu<br>
> <br>
<br>
</blockquote></div></div>