<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><blockquote type="cite" class=""><div class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><div class="">To me the best solution is to first split the files when possible (it is often not or just too painful) and then adding imports if needed. </div></div></blockquote><div class=""><br class=""></div>I’m not sure I agree. <span style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);" class="">Bear in mind my familiarity with fortran is very limited, but to me the “use” statement is similar to #include <header> in C or using namespace xxx in cpp.</span><div class=""><font color="#000000" class=""><span style="caret-color: rgb(0, 0, 0);" class=""><br class=""></span></font><div class="">Our fortran interface is not optimal code atm, it’s as if we put a "#include <petscvec.h>” before every Vec function in the C source (assuming include guards are not used). Sure you can fix this by splitting the vec src files in two, but at the end of the day that’s a bandaid solution. I think I can get the import stuff to work reasonably well (a lot of the infrastructure to capture types is already present in the current generatefortranstubs).<div class=""><div class=""><div class=""><br class=""><div class="">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div>Best regards,<br class=""><br class="">Jacob Faibussowitsch<br class="">(Jacob Fai - booss - oh - vitch)<br class="">Cell: (312) 694-3391</div></div>
</div>
<div><br class=""><blockquote type="cite" class=""><div class="">On Mar 16, 2021, at 12:56, Barry Smith <<a href="mailto:bsmith@petsc.dev" class="">bsmith@petsc.dev</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html; charset=utf-8" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div class=""><br class=""></div> Jacob,<div class=""><br class=""></div><div class=""> I split the mat definitions in the MR <a href="https://gitlab.com/petsc/petsc/-/merge_requests/3696" class="">https://gitlab.com/petsc/petsc/-/merge_requests/3696</a> and this reduced memory requirements enough to get builds through on some VM that failed previously ran out of memory (with gfortran). </div><div class=""><br class=""></div><div class=""> The petsc*mod.F90 files are all handwritten so it is ok to manually change them with imports if that helps the IBM compiler. To me the best solution is to first split the files when possible (it is often not or just too painful) and then adding imports if needed. </div><div class=""><br class=""></div><div class=""> Note also the MR <a href="https://gitlab.com/petsc/petsc/-/merge_requests/3715/diffs" class="">https://gitlab.com/petsc/petsc/-/merge_requests/3715/diffs</a> which may help a great deal with the IBM compiler or not.</div><div class=""><br class=""></div><div class=""> Thanks</div><div class=""><br class=""></div><div class=""> Barry</div><div class=""><br class=""></div><div class=""><br class=""><div class=""><br class=""><blockquote type="cite" class=""><div class="">On Mar 16, 2021, at 11:47 AM, Jacob Faibussowitsch <<a href="mailto:jacob.fai@gmail.com" class="">jacob.fai@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html; charset=utf-8" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Ok something I have gotten to work, but only doing things by hand in petscvecmod.F90:<div class=""><br class=""></div><div class="">diff --git a/src/vec/f90-mod/petscvecmod.F90 b/src/vec/f90-mod/petscvecmod.F90<br class="">index 0c447156b9..ef3e2e2e55 100644<br class="">--- a/src/vec/f90-mod/petscvecmod.F90<br class="">+++ b/src/vec/f90-mod/petscvecmod.F90<br class=""></div><div class=""> module petscisdef<br class=""> use petscisdefdummy<br class=""> interface operator(.ne.)</div><div class=""> function isnotequal(A,B)<br class="">- use petscisdefdummy</div><div class="">+ import tIS<br class=""> logical isnotequal <br class=""> type(tIS), intent(in) :: A,B<br class=""> end function isnotequal<br class=""><div class=""><br class=""></div><div class="">This works for everything wrapped in a module which contains an interface block. For dangling functions the following works:</div><div class=""><br class=""></div><div class=""> function isnotequal(A,B)</div><div class="">- use petscisdefdummy</div><div class="">+ use petscisdef, only: tIs<br class=""> logical isnotequal<br class=""> type(tIS), intent(in) :: A,B<br class=""> isnotequal = (A%v .ne. B%v)<br class=""> end function isnotequal</div><div class=""><br class=""></div><div class="">Do our fortran stub generators have any notion of types? Specifically types that originate from petsc?</div><div class=""><br class=""></div><div class=""><div class="">
<div dir="auto" style="caret-color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div class="">Best regards,<br class=""><br class="">Jacob Faibussowitsch<br class="">(Jacob Fai - booss - oh - vitch)<br class="">Cell: (312) 694-3391</div></div>
</div>
<div class=""><br class=""><blockquote type="cite" class=""><div class="">On Mar 16, 2021, at 11:11, Satish Balay <<a href="mailto:balay@mcs.anl.gov" class="">balay@mcs.anl.gov</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="">On Tue, 16 Mar 2021, Jacob Faibussowitsch wrote:<br class=""><br class=""><blockquote type="cite" class=""><blockquote type="cite" class="">My [partial] change is in branch balay/reorg-f90-for-xlf<br class=""></blockquote><br class="">Satish is this branch pushed? I’d like to send it to the ibm folks to see if it works for them as well.<br class=""></blockquote><br class="">Sorry - pushed now. From what I remember - it didn't work.<br class=""><br class="">Satish<br class=""><br class=""><br class=""><blockquote type="cite" class="">They also added this extra follow up:<br class=""><br class="">The change we did in the source files is to replace all the "use pet*" statements in all the Interface blocks with IMPORT statement.<br class=""><br class="">The nature of this workaround is to reduce the number of symbols that the compiler have to create, which exceeded the limitation and caused ICE.<br class=""><br class="">Every USE statement in an interface block opens up the module symbol file and reads all the symbols from it and creates an entity for each symbol in compiler. This is unnecessary when the module already has the same USE statement in the module scope. Instead, users can use IMPORT statement to make the module symbols accessible inside interface face blocks.<br class=""><br class="">With the change, the compiler would only read the module symbol file once and create one set of symbols where the old code reads the module symbol files as many times as the number of the USE statements in Interface blocks and create that many sets of duplicate symbols. Replacing those USE statements with IMPORT statements also shortens the compile time significantly.<br class=""><br class="">Best regards,<br class=""><br class="">Jacob Faibussowitsch<br class="">(Jacob Fai - booss - oh - vitch)<br class="">Cell: (312) 694-3391<br class=""><br class=""><blockquote type="cite" class="">On Mar 3, 2021, at 13:43, Satish Balay via petsc-dev <<a href="mailto:petsc-dev@mcs.anl.gov" class="">petsc-dev@mcs.anl.gov</a>> wrote:<br class=""><br class="">The only change I can get working (i.e avoid compile errors) is to split petscvecmod.F90 into 2 source files - but I don't know if this will help with xlf..<br class=""><br class="">My [partial] change is in branch balay/reorg-f90-for-xlf<br class=""><br class="">Satish<br class=""><br class="">On Wed, 3 Mar 2021, Satish Balay via petsc-dev wrote:<br class=""><br class=""><blockquote type="cite" class="">On Wed, 3 Mar 2021, Satish Balay via petsc-dev wrote:<br class=""><br class=""><blockquote type="cite" class="">Sure - once any change works locally [for gcc and xlf]<br class=""><br class="">When I try - I get a bunch of errors.. [yet to digest them.]<br class=""></blockquote><br class=""><blockquote type="cite" class=""><blockquote type="cite" class=""><blockquote type="cite" class=""><blockquote type="cite" class=""><blockquote type="cite" class="">Can you please give the following source code workaround a try?<br class="">Since there is already "use petscvecdefdummy" at the module scope, one workaround might be to remove the unnecessary "use petscvecdefdummy" in vecnotequal and vecequals <br class="">and all similar procedures.<br class=""><br class="">For example, the test case has:<br class=""> module petscvecdef<br class=""> use petscvecdefdummy<br class="">...<br class=""> function vecnotequal(A,B)<br class=""> use petscvecdefdummy<br class=""> logical vecnotequal<br class=""> type(tVec), intent(in) :: A,B<br class=""> vecnotequal = (A%v .ne. B%v)<br class=""> end function<br class=""></blockquote></blockquote></blockquote></blockquote></blockquote><br class=""><br class="">Ok - try this suggestion:<br class=""><br class="">diff --git a/src/vec/f90-mod/petscvecmod.F90 b/src/vec/f90-mod/petscvecmod.F90<br class="">index 0c447156b9..81968c7ca1 100644<br class="">--- a/src/vec/f90-mod/petscvecmod.F90<br class="">+++ b/src/vec/f90-mod/petscvecmod.F90<br class="">@@ -77,7 +77,6 @@<br class=""> use petscvecdefdummy<br class=""> interface operator(.ne.)<br class=""> function vecnotequal(A,B)<br class="">- use petscvecdefdummy<br class=""> logical vecnotequal<br class=""> type(tVec), intent(in) :: A,B<br class=""> end function<br class=""><br class=""><br class=""><blockquote type="cite" class=""><blockquote type="cite" class=""><blockquote type="cite" class=""><blockquote type="cite" class=""><blockquote type="cite" class=""><br class=""></blockquote></blockquote></blockquote></blockquote></blockquote> FC arch-linux-c-debug/obj/vec/f90-mod/petscvecmod.o<br class="">/home/balay/petsc/src/vec/f90-mod/petscvecmod.F90:81:22:<br class=""><br class=""> 81 | type(tVec), intent(in) :: A,B<br class=""> | 1<br class="">Error: Derived type ‘tvec’ at (1) is being used before it is defined<br class="">/home/balay/petsc/include/../src/vec/f90-mod/ftn-auto-interfaces/petscis.h90:2:10:<br class=""><br class=""> 2 | use petscvecdef<br class=""> | 1<br class="">Fatal Error: Cannot open module file ‘petscvecdef.mod’ for reading at (1): No such file or directory<br class=""><<<<<<<br class=""><br class="">Satish<br class=""></blockquote></blockquote><br class=""><br class=""></blockquote></div></div></blockquote></div><br class=""></div></div></div></div></blockquote></div><br class=""></div></div></div></blockquote></div><br class=""></div></div></div></div></div></body></html>