<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Feb 20, 2018 at 9:40 PM, TAY wee-beng <span dir="ltr"><<a href="mailto:zonexo@gmail.com" target="_blank">zonexo@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi,<br>
<br>
Indeed, replacing tvec with t_vec solves the problem. Now I'm trying to debug step by step. I got into problem when calling:<br>
<br>
call DMCreateGlobalVector(da_u,u_gl<wbr>obal,ierr)<br>
<br>
The error is:<br>
<br>
[0]PETSC ERROR: --------------------- Error Message ----------------------------<br>
------------------------------<wbr>----<br>
[0]PETSC ERROR: Null argument, when expecting valid pointer<br>
[0]PETSC ERROR: Null Object: Parameter # 2<br>
[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" rel="noreferrer" target="_blank">http://www.mcs.anl.gov/petsc/d<wbr>ocumentation/faq.html</a> for trou<br>
ble shooting.<br>
[0]PETSC ERROR: Petsc Release Version 3.8.3, Dec, 09, 2017<br>
[0]PETSC ERROR: C:\Obj_tmp\ibm3d_IIB_mpi\Debug<wbr>\ibm3d_IIB_mpi.exe on a petsc-3.8.<br>
3_win64_msmpi_vs2008 named 1C3YYY1-PC by tsltaywb Wed Feb 21 10:20:20 2018<br>
[0]PETSC ERROR: Configure options --with-cc="win32fe icl" --with-fc="win32fe ifo<br>
rt" --with-cxx="win32fe icl" --download-fblaslapack --with-mpi-include="[/cygdri<br>
ve/c/Program Files (x86)/Microsoft SDKs/MPI/Include,/cygdrive/c/P<wbr>rogram Files (x....<br>
<br>
But all I changed is from:<br>
<br>
module global_data<br>
#include "petsc/finclude/petsc.h"<br>
use petsc<br>
use kdtree2_module<br>
implicit none<br>
save<br>
!grid variables<br>
<br>
integer :: size_x,s....<br>
<br>
...<br>
<br>
to<br>
<br>
module global_data<br>
use kdtree2_module<br>
implicit none<br>
save<br>
#include "petsc/finclude/petsc.h90"<br>
!grid variables<br>
integer :: size_x,s...<br>
<br>
...<br>
<br>
da_u, u_global were declared thru:<br>
<br>
DM da_u,da_v,...<br>
DM da_cu_types ...<br>
Vec u_local,u_global,v_local...<br>
<br>
So what could be the problem?<br></blockquote><div><br></div><div>If you are using the latest release (or master) then</div><div><br></div><div> <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Sys/UsingFortran.html">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Sys/UsingFortran.html</a></div><div><br></div><div>you include petsc.h not petsc.h90, and you need to 'use petsc' as well, so the first form</div><div>looks correct, not the second.</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
Thank you very much.<br>
<br>
Yours sincerely,<br>
<br>
==============================<wbr>==================<br>
TAY Wee-Beng (Zheng Weiming) 郑伟明<br>
Personal research webpage: <a href="http://tayweebeng.wixsite.com/website" rel="noreferrer" target="_blank">http://tayweebeng.wixsite.com/<wbr>website</a><br>
Youtube research showcase: <a href="https://www.youtube.com/channel/UC72ZHtvQNMpNs2uRTSToiLA" rel="noreferrer" target="_blank">https://www.youtube.com/channe<wbr>l/UC72ZHtvQNMpNs2uRTSToiLA</a><br>
linkedin: <a href="http://www.linkedin.com/in/tay-weebeng" rel="noreferrer" target="_blank">www.linkedin.com/in/tay-weeben<wbr>g</a><br>
==============================<wbr>==================<br>
<br>
On 20/2/2018 10:46 PM, Jose E. Roman wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Probably the first error is produced by using a variable (mpi_comm) with the same name as an MPI type.<br>
<br>
The second error I guess is due to variable tvec, since a Fortran type tVec is now being defined in src/vec/f90-mod/petscvec.h<br>
<br>
Jose<br>
<br>
<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
El 20 feb 2018, a las 15:35, Smith, Barry F. <<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>> escribió:<br>
<br>
<br>
Please run a clean compile of everything and cut and paste all the output. This will make it much easier to debug than trying to understand your snippets of what is going wrong.<br>
<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
On Feb 20, 2018, at 1:56 AM, TAY Wee Beng <<a href="mailto:tsltaywb@nus.edu.sg" target="_blank">tsltaywb@nus.edu.sg</a>> wrote:<br>
<br>
Hi,<br>
<br>
I was previously using PETSc 3.7.6 on different clusters with both Intel<br>
Fortran and GNU Fortran. After upgrading, I met some problems when<br>
trying to compile:<br>
<br>
On Intel Fortran:<br>
<br>
Previously, I was using:<br>
<br>
#include "petsc/finclude/petsc.h90"<br>
<br>
in *.F90 when requires the use of PETSc<br>
<br>
I read in the change log that h90 is no longer there and so I replaced<br>
with #include "petsc/finclude/petsc.h"<br>
<br>
It worked. But I also have some *.F90 which do not use PETSc. However,<br>
they use some modules which uses PETSc.<br>
<br>
Now I can't compile them. The error is :<br>
<br>
math_routine.f90(3): error #7002: Error in opening the compiled module<br>
file. Check INCLUDE paths. [PETSC]<br>
use mpi_subroutines<br>
<br>
mpi_subroutines is a module which uses PETSc, and it compiled w/o problem.<br>
<br>
The solution is that I have to compile e.g. math_routine.F90 as if they<br>
use PETSc, by including PETSc include and lib files.<br>
<br>
May I know why this is so? It was not necessary before.<br>
<br>
Anyway, it managed to compile until it reached hypre.F90.<br>
<br>
Previously, due to some bugs, I have to compile hypre with the -r8<br>
option. Also, I have to use:<br>
<br>
integer(8) mpi_comm<br>
<br>
mpi_comm = MPI_COMM_WORLD<br>
<br>
to make my codes work with HYPRE.<br>
<br>
But now, compiling gives the error:<br>
<br>
hypre.F90(11): error #6401: The attributes of this name conflict with<br>
those made accessible by a USE statement. [MPI_COMM]<br>
integer(8) mpi_comm<br>
------------------------------<wbr>--------^<br>
hypre.F90(84): error #6478: A type-name must not be used as a<br>
variable. [MPI_COMM]<br>
mpi_comm = MPI_COMM_WORLD<br>
----^<br>
hypre.F90(84): error #6303: The assignment operation or the binary<br>
expression operation is invalid for the data types of the two<br>
operands. [1140850688]<br>
mpi_comm = MPI_COMM_WORLD<br>
---------------^<br>
hypre.F90(100): error #6478: A type-name must not be used as a<br>
variable. [MPI_COMM]<br>
call HYPRE_StructGridCreate(mpi_com<wbr>m, 3, grid_hypre, ierr)<br>
...<br>
<br>
What's actually happening? Why can't I compile now?<br>
<br>
On GNU gfortran:<br>
<br>
I tried to use similar tactics as above here. However, when compiling<br>
math_routine.F90, I got the error:<br>
<br>
math_routine.F90:1333:21:<br>
<br>
call subb(orig,vert1,tvec)<br>
1<br>
Error: Invalid procedure argument at (1)<br>
math_routine.F90:1339:18:<br>
<br>
qvec = cross_pdt2(tvec,edge1)<br>
1<br>
Error: Invalid procedure argument at (1)<br>
math_routine.F90:1345:21:<br>
<br>
uu = dot_product(tvec,pvec)<br>
1<br>
Error: ‘vector_a’ argument of ‘dot_product’ intrinsic at (1) must be<br>
numeric or LOGICAL<br>
math_routine.F90:1371:21:<br>
<br>
uu = dot_product(tvec,pvec)<br>
<br>
These errors were not present before. My variables are mostly vectors:<br>
<br>
real(8), intent(in) ::<br>
orig(3),infinity(3),vert1(3),v<wbr>ert2(3),vert3(3),normal(3)<br>
<br>
real(8) :: uu,vv,dir(3)<br>
<br>
real(8) :: edge1(3),edge2(3),tvec(3),pvec<wbr>(3),qvec(3),det,inv_det,epsilo<wbr>n,d,t<br>
<br>
I wonder what happened?<br>
<br>
Please advice.<br>
<br>
<br>
--<br>
Thank you very much.<br>
<br>
Yours sincerely,<br>
<br>
==============================<wbr>==================<br>
TAY Wee-Beng 郑伟明<br>
Research Scientist<br>
Experimental AeroScience Group<br>
Temasek Laboratories<br>
National University of Singapore<br>
T-Lab Building<br>
5A, Engineering Drive 1, #02-02<br>
Singapore 117411<br>
Phone: <a href="tel:%2B65%2065167330" value="+6565167330" target="_blank">+65 65167330</a><br>
E-mail: <a href="mailto:tsltaywb@nus.edu.sg" target="_blank">tsltaywb@nus.edu.sg</a><br>
<a href="http://www.temasek-labs.nus.edu.sg/program/program_aeroexperimental_tsltaywb.php" rel="noreferrer" target="_blank">http://www.temasek-labs.nus.ed<wbr>u.sg/program/program_aeroexper<wbr>imental_tsltaywb.php</a><br>
Personal research webpage: <a href="http://tayweebeng.wixsite.com/website" rel="noreferrer" target="_blank">http://tayweebeng.wixsite.com/<wbr>website</a><br>
Youtube research showcase: <a href="https://www.youtube.com/channel/UC72ZHtvQNMpNs2uRTSToiLA" rel="noreferrer" target="_blank">https://www.youtube.com/channe<wbr>l/UC72ZHtvQNMpNs2uRTSToiLA</a><br>
linkedin: <a href="http://www.linkedin.com/in/tay-weebeng" rel="noreferrer" target="_blank">www.linkedin.com/in/tay-weeben<wbr>g</a><br>
==============================<wbr>==================<br>
<br>
<br>
______________________________<wbr>__<br>
<br>
Important: This email is confidential and may be privileged. If you are not the intended recipient, please delete it and notify us immediately; you should not copy or use it for any purpose, nor disclose its contents to any other person. Thank you.<br>
</blockquote></blockquote></blockquote>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>