src/EfieldIntegrand.h

Go to the documentation of this file.

//
// C++ Interface: EfieldIntegrand
//
// Description: 
//
//
// Author: Roman Geus <geus@maxwell>, (C) 2004
//
// Copyright: See COPYING file that comes with this distribution
//
//

#ifndef EfieldIntegrand_H
#define EfieldIntegrand_H

#include "vector3.h"
#include "vector.h"
#include "nedelecmesh.h"

template <typename Curve>
class EfieldDs {
public:
    EfieldDs(Curve curve,
             const NedelecMesh& nedelec_mesh,
             const colarray::Vector<double>& q)
        : curve_(curve),
          nedelec_mesh_(nedelec_mesh),
          q_(q)
    { }
    double operator()(const double t) {
        mesh::Vector3 field = nedelec_mesh_.eval(curve_(t), q_);
        return field.dot_product(curve_.deriv(t));
    }
private:
    Curve curve_;
    const NedelecMesh& nedelec_mesh_;
    const colarray::Vector<double>& q_;
};

class EfieldIntegrand {
public:
    EfieldIntegrand(const mesh::Vector3& start,
                    const mesh::Vector3& direction,
                    const NedelecMesh& nedelec_mesh,
                    const colarray::Vector<double>& q)
        : start_(start),
          direction_(direction),
          nedelec_mesh_(nedelec_mesh),
          q_(q)
    { }
    double operator()(const double x) {
        mesh::Vector3 p(start_ + x*direction_);
        mesh::Vector3 field = nedelec_mesh_.eval(p, q_);
        return field.dot_product(direction_);
    }
private:
    const mesh::Vector3& start_;
    const mesh::Vector3& direction_;
    const NedelecMesh& nedelec_mesh_;
    const colarray::Vector<double>& q_;
};

#endif

---