vectorBCFunction.cc

Go to the documentation of this file.

/*
 * vectorBCFunction.cc
 *
 *  Created on: Feb 22, 2017
 *      Author: stephendewitt
 */

#include <deal.II/base/function.h>
#include <deal.II/lac/parallel_vector.h>
#include "../../include/vectorBCFunction.h"

template <int dim>
vectorBCFunction<dim>::vectorBCFunction(std::vector<double> input_values) : dealii::Function<dim>(dim), BC_values (input_values) {}

template <int dim>
void vectorBCFunction<dim>::vector_value(const dealii::Point<dim> &p, dealii::Vector<double> &values) const {

    for (unsigned int i=0; i<dim; i++) {
        values(i) = BC_values[i];
    }
}

template <int dim>
void vectorBCFunction<dim>::vector_value_list (const std::vector<dealii::Point<dim> > &points, std::vector<dealii::Vector<double> > &value_list) const{
    const unsigned int n_points = points.size();
    for (unsigned int p=0; p<n_points; ++p)
        vectorBCFunction<dim>::vector_value(points[p],value_list[p]);
}


template class vectorBCFunction<1>;
template class vectorBCFunction<2>;
template class vectorBCFunction<3>;