gsINSTerms.hpp

Go to the documentation of this file.

 
#pragma once
#include <gsIncompressibleFlow/src/gsINSTerms.h>
 
namespace gismo
{
 
 
template<class T>
void gsINSTerm_PvalUdiv<T>::assemble(const gsMapData<T>& mapData, const gsVector<T>& quWeights, const std::vector< gsMatrix<T> >& testFunData, const std::vector< gsMatrix<T> >& shapeFunData, std::vector< gsMatrix<T> >& localMat)
{ 
    gsVector<T> coeffMeasure = this->getCoeffGeoMapProduct(mapData);
 
    const gsMatrix<T>& testFunGrads = testFunData[1];
    const gsMatrix<T>& shapeFunVals = shapeFunData[0];
 
    gsMatrix<T> testFunPhysGrad;
 
    const index_t nQuPoints = quWeights.rows();
 
    for (index_t k = 0; k < nQuPoints; k++)
    {
        const T weight = quWeights(k) * coeffMeasure(k);
 
        transformGradients(mapData, k, testFunGrads, testFunPhysGrad);
 
        for (size_t i = 0; i != localMat.size(); ++i)
            localMat[i].noalias() += weight * (shapeFunVals.col(k) * testFunPhysGrad.row(i)).transpose();
    }
}
 
// ===================================================================================================================
 
template<class T>
void gsINSTerm_UdivPval<T>::assemble(const gsMapData<T>& mapData, const gsVector<T>& quWeights, const std::vector< gsMatrix<T> >& testFunData, const std::vector< gsMatrix<T> >& shapeFunData, std::vector< gsMatrix<T> >& localMat)
{ 
    gsVector<T> coeffMeasure = this->getCoeffGeoMapProduct(mapData);
 
    const gsMatrix<T>& testFunVals = testFunData[0];
    const gsMatrix<T>& shapeFunGrads = shapeFunData[1];
 
    gsMatrix<T> shapeFunPhysGrad;
 
    const index_t nQuPoints = quWeights.rows();
 
    for (index_t k = 0; k < nQuPoints; k++)
    {
        const T weight = quWeights(k) * coeffMeasure(k);
 
        transformGradients(mapData, k, shapeFunGrads, shapeFunPhysGrad);
 
        for (size_t i = 0; i != localMat.size(); ++i)
            localMat[i].noalias() += weight * (shapeFunPhysGrad.row(i) * testFunVals(k) );
    }
}
 
// ===================================================================================================================
 
template<class T>
void gsINSTerm_UsolGradVal<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)
{ 
    this->computeCoeffSolU(mapData);
    gsVector<T> coeffMeasure = this->getCoeffGeoMapProduct(mapData);
 
    const gsMatrix<T>& testFunVals = testFunData[0];
    const gsMatrix<T>& shapeFunGrads = shapeFunData[1];
 
    gsMatrix<T> shapeFunPhysGrad;
 
    const index_t nQuPoints = quWeights.rows();
 
    for (index_t k = 0; k < nQuPoints; k++)
    {
        const T weight = quWeights(k) * coeffMeasure(k);
 
        transformGradients(mapData, k, shapeFunGrads, shapeFunPhysGrad);
 
        localMat.noalias() += weight * (testFunVals.col(k) * (this->m_solUVals.col(k).transpose() * shapeFunPhysGrad));
    }
}
 
 
} // namespace gismo