fixup! Added: virtual method ASMbase::findPoints() searching for a list of spatial points. The default implementation invokes ASMbase::findPoint() in a loop over the specified points. It is overridden in ASMu2DLag to perform parametric inversion for bilinear quads and triangles.

Author: kmokstad

src/ASM/ASMu2DLag.C

@@ -39,6 +39,10 @@ extern "C" {
   //! \brief Evaluate the bilinear shape functions.
   void shape2_(const int& iop, const int& nenod, const double* xi,
                double* shp, const int& ipsw, const int& iwr, int& ierr);
+  //! \brief Opens a temporary file for Fortran print.
+  int openftnfile_(const char* fname, const ssize_t nchar);
+  //! \brief Closes a Fortran file.
+  void closeftnfile_(const int& iunit);
 }
 
 
@@ -769,19 +773,30 @@ bool ASMu2DLag::paramInvert (Vec3& X, PointParams& elm) const
   const double eps = 1.0e-15;
   const double tol = 1.0e-3;
 #if SP_DEBUG > 2
-  const int ipsw = SP_DEBUG;
+  int ipsw = SP_DEBUG;
 #else
-  const int ipsw = 0;
+  int ipsw = 0;
 #endif
-  const int iwr = 6;
+  // Open a temporary file for the Fortran output
+  std::string tmpfile = "/tmp/" + std::string(getenv("USER")) + ".ftn";
+  const int iwr = openftnfile_(tmpfile.c_str(),tmpfile.size());
+  if (iwr <= 0) ipsw = -1; // suppress all output
   int ierr = 0;
   if (nelnod == 3)
     t3_inv_(nsd,tol,eps,X.ptr(),Xn.ptr(),elm.u,ipsw,iwr,ierr);
   else if (nelnod == 4)
     q4_inv_(nsd,tol,eps,X.ptr(),Xn.ptr(),elm.u,ipsw,iwr,ierr);
   else
     ierr = -99;
-
+  if (iwr > 0 && ipsw >= 0)
+  {
+    // Copy fortran output to IFEM::cout
+    closeftnfile_(iwr);
+    std::string cline;
+    std::ifstream is(tmpfile);
+    while (std::getline(is,cline))
+      IFEM::cout << cline << std::endl;
+  }
   if (ierr != 0 && ierr != 2)
     return false;