fn_inv2.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_INV2_HPP
#define BELFEM_FN_INV2_HPP
 
#include "typedefs.hpp"
#include "assert.hpp"
#include "cl_Matrix.hpp"
 
 
namespace belfem
{
//------------------------------------------------------------------------------

    inline real
    inv2( const Matrix< real > & aA, Matrix< real > & aB )
    {
        // compute determinant
        real aDetA = aA( 0, 0 )*aA( 1,1 )
                   - aA( 1, 0 )*aA( 0,1 );
 
        aB( 0, 0 ) =   aA( 1, 1 );
        aB( 1, 0 ) = - aA( 1, 0 );
        aB( 0, 1 ) = - aA( 0, 1 );
        aB( 1, 1 ) =   aA( 0, 0 );
 
        // singularity test on the hadamard ratio |det| / prod( ||row_i|| ),
        // which is dimensionless: an absolute bound on the determinant is a
        // scale test, and rejects well conditioned but small matrices such as
        // any fine mesh jacobian. squared here to avoid the square roots
        BELFEM_ASSERT( aDetA * aDetA > BELFEM_EPSILON * BELFEM_EPSILON
                       * ( aA( 0, 0 ) * aA( 0, 0 ) + aA( 0, 1 ) * aA( 0, 1 ) )
                       * ( aA( 1, 0 ) * aA( 1, 0 ) + aA( 1, 1 ) * aA( 1, 1 ) ),
                       "Can't invert a singular matrix" );
 
        aB /= aDetA ;
 
        return aDetA ;
    }
 
//------------------------------------------------------------------------------
}
#endif //BELFEM_FN_INV2_HPP