BEM_Poisson_Boltzmann.h

Go to the documentation of this file.

// 
//  This file is part of the FFEA simulation package
//  
//  Copyright (c) by the Theory and Development FFEA teams,
//  as they appear in the README.md file. 
// 
//  FFEA is free software: you can redistribute it and/or modify
//  it under the terms of the GNU General Public License as published by
//  the Free Software Foundation, either version 3 of the License, or
//  (at your option) any later version.
// 
//  FFEA is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.
// 
//  You should have received a copy of the GNU General Public License
//  along with FFEA.  If not, see <http://www.gnu.org/licenses/>.
// 
//  To help us fund FFEA development, we humbly ask that you cite 
//  the research papers on the package.
//

/*
 *      BEM_Poisson_Boltzmann.h
 *  Author: Robin Richardson, University of Leeds
 *  Email: pyrar@leeds.ac.uk
 */

#ifndef BEM_POISSON_BOLTZMANN_H_INCLUDED
#define BEM_POISSON_BOLTZMANN_H_INCLUDED

#include <math.h>

#include "FFEA_return_codes.h"
#include "NearestNeighbourLinkedListCube.h"
#include "SparseMatrixUnknownPattern.h"

#include "GaussianQuadrature_tri.h"
#include "GaussianQuadrature_1d.h"

class BEM_Poisson_Boltzmann : GaussianQuadrature_tri, GaussianQuadrature_1d {
public:
    BEM_Poisson_Boltzmann();
    ~BEM_Poisson_Boltzmann();
    int init(NearestNeighbourLinkedListCube *lookup);
    void set_kappa(scalar kappa);
    void build_BEM_matrices();
    void perform_integrals_for_lookup_cell_self(LinkedListNode<Face> *l_i, vector3 gqp[4]);
    void perform_integrals_for_lookup_cell_relative(LinkedListNode<Face> *l_i, vector3 gqp[4], int dx, int dy, int dz);
    void print_matrices();
    SparseMatrixUnknownPattern * get_C();
    SparseMatrixUnknownPattern * get_D();

private:

    NearestNeighbourLinkedListCube *lookup;

    int num_faces;


    SparseMatrixUnknownPattern *mat_C, *mat_D;

    scalar kappa;

    scalar u_4pi(scalar r);

    scalar grad_u_4pi(scalar r, scalar r2);

    /*
                    scalar screened_R_theta(scalar r_perp_mag, scalar half_theta_max, scalar theta_bar, scalar xi);
     */

    void gauss_quadrature_4_point(vector3 gqp[4], vector3 *p, scalar *int_u, scalar *int_du, Face *f);

    scalar self_term(vector3 *n0, vector3 *n1, vector3 *n2, int precision);

    scalar f_1d(scalar r);

    scalar f_3d(vector3 *p, vector3 *q);
};

#endif