fn_maxwell_l2N.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 FN_MAXWELL_L2N_HPP
#define FN_MAXWELL_L2N_HPP
 
#include "typedefs.hpp"
#include "cl_Matrix.hpp"
 
#include "cl_FEM_Calculator.hpp"
 
namespace belfem
{
    namespace fem
    {
        namespace maxwell
        {
            inline const Matrix< real > &
            l2N( Calculator * aCalc, const uint aIndex )
            {
                const Vector< real > & n = aCalc->Nvec( aIndex ) ;
                Matrix< real > & N = aCalc->matrix("L2N");
 
                uint nnodes = n.length();
                uint ndim = N.n_rows() ;
                for( uint j=0; j<ndim; ++j )
                {
                    uint c = j ;
                    for( uint i=0; i<nnodes; ++i )
                    {
                        N( j, c ) = n( i );
                        c+=ndim;
                    }
                }
 
                return N;
            }
        }
    }
}
#endif //FN_MAXWELL_L2N_HPP