TairaColoniusSolver.h

Go to the documentation of this file.

#pragma once

#include "NSWithBody.h"


template <typename memoryType>
class TairaColoniusSolver : public NSWithBody<memoryType>
{
private:
    cusp::coo_matrix<int, real, memoryType> 
      E,   
      ET;  
    
    // generate the transposed gradient matrix and the interpolation matrix
    virtual void generateQT();

    // update the interpolation matrix with current location of body points
    void updateQT();

    // generate the right hand-side of the Poisson system
    virtual void generateBC2();
    
    // update the location of the bodies and re-generate appropriate matrices
    virtual void updateSolverState();
    
    // calculate forces acting on each immersed body
    virtual void calculateForce();
    
public:
    // constructor -- copy the database and information about the computational grid
    TairaColoniusSolver(parameterDB *pDB=NULL, domain *dInfo=NULL);
    
    // initiliaze the solvers
    virtual void initialise();
    
    // calculate and write forces acting on each immersed body at current time
    virtual void writeData();
    
    virtual std::string name()
    {
        return "Taira & Colonius";
    }
}; // TairaColoniusSolver