Merge pull request #71 from robertodr/fix-fortran-api

Fix API for the surface function NumPy save/load functions

Author: bast

src/interface/Meddle.cpp

@@ -351,8 +351,7 @@ void pcmsolver_load_surface_function(pcmsolver_context_t * context,
 void pcm::Meddle::loadSurfaceFunction(const char * name) const {
   std::string functionName(name);
   hostWriter_("\nLoading surface function " + functionName + " from .npy file");
-  std::string fname = functionName + ".npy";
-  Eigen::VectorXd values = cnpy::custom::npy_load<double>(fname);
+  Eigen::VectorXd values = cnpy::custom::npy_load<double>(functionName + ".npy");
   if (values.size() != cavity_->size())
     PCMSOLVER_ERROR("The loaded " + functionName +
                     " surface function is bigger than the cavity!");

tests/C_host/C_host.c

@@ -157,8 +157,9 @@ int main() {
   test_surface_functions(
       output, grid_size, mep, asc_Ag, asc_B3g, asc_neq_B3g, areas);
 
-  pcmsolver_save_surface_function(pcm_context, mep_lbl);
+  pcmsolver_save_surface_functions(pcm_context);
   pcmsolver_save_surface_function(pcm_context, asc_lbl);
+  pcmsolver_load_surface_function(pcm_context, mep_lbl);
 
   pcmsolver_write_timings(pcm_context);
 

tests/Fortran_host/Fortran_host.f90

@@ -168,6 +168,9 @@ program pcm_fortran_host
       ! Surface functions
       call test_surface_functions(grid_size, mep, asc_Ag, asc_B3g, asc_neq_B3g, areas)
 
+      call pcmsolver_save_surface_function(pcm_context, mep_lbl)
+      call pcmsolver_load_surface_function(pcm_context, mep_lbl)
+
       call pcmsolver_write_timings(pcm_context)
 
       call pcmsolver_delete(pcm_context)