gsFlowTerms.hpp

Go to the documentation of this file.

 
#pragma once
#include <gsIncompressibleFlow/src/gsFlowTerms.h>
 
namespace gismo
{
 
template<class T>
gsVector<T> gsFlowTerm<T>::getCoeffGeoMapProduct(const gsMapData<T>& mapData)
{
    evalCoeff(mapData);
 
    index_t npts = mapData.points.cols();
    gsVector<T> result(npts);
 
    for (index_t k = 0; k < npts; k++)
    {
        index_t coeffID = (m_coeff.size() == 1) ? 0 : k;
        result(k) = m_coeff(coeffID) * mapData.measure(k);
    }
 
    return result;
}
 
 
template<class T>
void gsFlowTerm_ValVal<T>::assemble(const gsMapData<T>& mapData, const gsVector<T>& quWeights, const std::vector< gsMatrix<T> >& testFunData, const std::vector< gsMatrix<T> >& shapeFunData, gsMatrix<T>& localMat)
{ 
    gsVector<T> coeffMeasure = this->getCoeffGeoMapProduct(mapData);
 
    const gsMatrix<T>& testFunVals = testFunData[0];
    const gsMatrix<T>& shapeFunVals = shapeFunData[0];
 
    const index_t nQuPoints = quWeights.rows();
 
    for (index_t k = 0; k < nQuPoints; k++)
    {
        const T weight = quWeights(k) * coeffMeasure(k);
        localMat.noalias() += weight * (testFunVals.col(k) * shapeFunVals.col(k).transpose());
    }
}
 
// ===================================================================================================================
 
template<class T>
void gsFlowTerm_GradGrad<T>::assemble(const gsMapData<T>& mapData, const gsVector<T>& quWeights, const std::vector< gsMatrix<T> >& testFunData, const std::vector< gsMatrix<T> >& shapeFunData, gsMatrix<T>& localMat)
{ 
    gsVector<T> coeffMeasure = this->getCoeffGeoMapProduct(mapData);
 
    const gsMatrix<T>& testFunGrads = testFunData[1];
    const gsMatrix<T>& shapeFunGrads = shapeFunData[1];
 
    const index_t nQuPoints = quWeights.rows();
    gsMatrix<T> testFunPhysGrad, shapeFunPhysGrad;
 
    for (index_t k = 0; k < nQuPoints; k++)
    {
        const T weight = quWeights(k) * coeffMeasure(k);
 
        transformGradients(mapData, k, testFunGrads, testFunPhysGrad);
        transformGradients(mapData, k, shapeFunGrads, shapeFunPhysGrad);
 
        localMat.noalias() += weight * (testFunPhysGrad.transpose() * shapeFunPhysGrad);
    }
}
 
// ===================================================================================================================
 
template<class T>
void gsFlowTerm_rhs<T>::assemble(const gsMapData<T>& mapData, const gsVector<T>& quWeights, const std::vector< gsMatrix<T> >& testFunData, const std::vector< gsMatrix<T> >& shapeFunData, gsMatrix<T>& localMat)
{ 
    m_pRhsFun->eval_into(mapData.values[0], m_rhsVals);
 
    const index_t nQuPoints = quWeights.rows();
 
    for (index_t k = 0; k < nQuPoints; k++)
    {
        const T weight = quWeights(k) * mapData.measure(k);
 
        localMat.noalias() += weight * (testFunData[0].col(k) *  m_rhsVals.col(k).transpose());
    }
}
 
 
} // namespace gismo