fn_rho_database_is_current.hpp

Go to the documentation of this file.

/*
 * BELFEM -- The Berkeley Lab Finite Element Framework
 * Copyright (c) 2026, The Regents of the University of California, through
 * Lawrence Berkeley National Laboratory (subject to receipt of any required
 * approvals from the U.S. Dept. of Energy).  All rights reserved.
 *
 * Developers: Christian Messe, Gregory Giard
 *
 * See the top-level LICENSE file for the complete license and disclaimer.
 */
 
#ifndef BELFEM_FN_RHO_DATABASE_IS_CURRENT_HPP
#define BELFEM_FN_RHO_DATABASE_IS_CURRENT_HPP
 
#include "typedefs.hpp"
#include "hdf5_tools.hpp"
 
namespace belfem
{
    inline bool
    rho_database_is_current( const string & aPath )
    {
#ifdef BELFEM_HDF5
        hid_t tFile = H5Fopen( aPath.c_str(), H5F_ACC_RDONLY, H5P_DEFAULT );
 
        if ( tFile < 0 )
        {
            return false;
        }
 
        // H5Lexists returns a tri-state: > 0 exists, 0 absent, < 0 error.
        // Only a definite hit counts as current -- a probe ERROR must land
        // on the rebuild side, not be swallowed by a bool conversion
        htri_t tExists = H5Lexists( tFile, "RRR", H5P_DEFAULT );
 
        H5Fclose( tFile );
 
        return tExists > 0;
#else
        // without HDF5 the database path cannot work at all; keep the
        // pre-existing behavior rather than forcing a rebuild here
        return true;
#endif
    }
}
#endif //BELFEM_FN_RHO_DATABASE_IS_CURRENT_HPP