general_hypersingular_integral_operator.hpp

// Copyright (C) 2011-2012 by the BEM++ Authors
//
// Permission is hereby granted, free of charge, to any person obtaining a copy
// of this software and associated documentation files (the "Software"), to deal
// in the Software without restriction, including without limitation the rights
// to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
// copies of the Software, and to permit persons to whom the Software is
// furnished to do so, subject to the following conditions:
//
// The above copyright notice and this permission notice shall be included in
// all copies or substantial portions of the Software.
//
// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
// IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
// AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
// OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
// THE SOFTWARE.

#ifndef bempp_general_hypersingular_integral_operator_hpp
#define bempp_general_hypersingular_integral_operator_hpp

#include "hypersingular_integral_operator.hpp"

namespace Bempp
{

template <typename BasisFunctionType_, typename KernelType_, typename ResultType_>
class GeneralHypersingularIntegralOperator :
        public HypersingularIntegralOperator<
            BasisFunctionType_, KernelType_, ResultType_>
{
    typedef HypersingularIntegralOperator<
    BasisFunctionType_, KernelType_, ResultType_> Base;
public:
    typedef typename Base::BasisFunctionType BasisFunctionType;
    typedef typename Base::KernelType KernelType;
    typedef typename Base::ResultType ResultType;
    typedef typename Base::CoordinateType CoordinateType;
    typedef typename Base::CollectionOfShapesetTransformations
    CollectionOfShapesetTransformations;
    typedef typename Base::CollectionOfBasisTransformations
    CollectionOfBasisTransformations;
    typedef typename Base::CollectionOfKernels CollectionOfKernels;
    typedef typename Base::TestKernelTrialIntegral TestKernelTrialIntegral;

    template <typename KernelFunctor,
              typename TestTransformationsFunctor,
              typename TrialTransformationsFunctor,
              typename IntegrandFunctor,
              typename OffDiagonalKernelFunctor,
              typename OffDiagonalTestTransformationsFunctor,
              typename OffDiagonalTrialTransformationsFunctor,
              typename OffDiagonalIntegrandFunctor>
    GeneralHypersingularIntegralOperator(
            const shared_ptr<const Space<BasisFunctionType_> >& domain,
            const shared_ptr<const Space<BasisFunctionType_> >& range,
            const shared_ptr<const Space<BasisFunctionType_> >& dualToRange,
            const std::string& label,
            int symmetry,
            const KernelFunctor& kernelFunctor,
            const TestTransformationsFunctor& testTransformationsFunctor,
            const TrialTransformationsFunctor& trialTransformationsFunctor,
            const IntegrandFunctor& integrandFunctor,
            const OffDiagonalKernelFunctor& offDiagonalKernelFunctor,
            const OffDiagonalTestTransformationsFunctor& offDiagonalTestTransformationsFunctor,
            const OffDiagonalTrialTransformationsFunctor& offDiagonalTrialTransformationsFunctor,
            const OffDiagonalIntegrandFunctor& offDiagonalIntegrandFunctor);

    template <typename KernelFunctor,
              typename TestTransformationsFunctor,
              typename TrialTransformationsFunctor,
              typename OffDiagonalKernelFunctor,
              typename OffDiagonalTestTransformationsFunctor,
              typename OffDiagonalTrialTransformationsFunctor>
    GeneralHypersingularIntegralOperator(
            const shared_ptr<const Space<BasisFunctionType_> >& domain,
            const shared_ptr<const Space<BasisFunctionType_> >& range,
            const shared_ptr<const Space<BasisFunctionType_> >& dualToRange,
            const std::string& label,
            int symmetry,
            const KernelFunctor& kernelFunctor,
            const TestTransformationsFunctor& testTransformationsFunctor,
            const TrialTransformationsFunctor& trialTransformationsFunctor,
            const shared_ptr<Fiber::TestKernelTrialIntegral<
            BasisFunctionType_, KernelType_, ResultType_> >& integral,
            const OffDiagonalKernelFunctor& offDiagonalKernelFunctor,
            const OffDiagonalTestTransformationsFunctor& offDiagonalTestTransformationsFunctor,
            const OffDiagonalTrialTransformationsFunctor& offDiagonalTrialTransformationsFunctor,
            const shared_ptr<Fiber::TestKernelTrialIntegral<
            BasisFunctionType_, KernelType_, ResultType_> >& offDiagonalIntegral);

    virtual const CollectionOfKernels& kernels() const
    { return *m_kernels; }
    virtual const CollectionOfShapesetTransformations& testTransformations() const
    { return *m_testTransformations; }
    virtual const CollectionOfShapesetTransformations& trialTransformations() const
    { return *m_trialTransformations; }
    virtual const TestKernelTrialIntegral& integral() const
    { return *m_integral; }

    virtual const CollectionOfKernels& offDiagonalKernels() const
    { return *m_offDiagonalKernels; }
    virtual const CollectionOfShapesetTransformations& offDiagonalTestTransformations() const
    { return *m_offDiagonalTestTransformations; }
    virtual const CollectionOfShapesetTransformations& offDiagonalTrialTransformations() const
    { return *m_offDiagonalTrialTransformations; }
    virtual const TestKernelTrialIntegral& offDiagonalIntegral() const
    { return *m_offDiagonalIntegral; }

private:
    shared_ptr<CollectionOfKernels> m_kernels;
    shared_ptr<CollectionOfShapesetTransformations> m_testTransformations;
    shared_ptr<CollectionOfShapesetTransformations> m_trialTransformations;
    shared_ptr<TestKernelTrialIntegral> m_integral;
    shared_ptr<CollectionOfKernels> m_offDiagonalKernels;
    shared_ptr<CollectionOfShapesetTransformations> m_offDiagonalTestTransformations;
    shared_ptr<CollectionOfShapesetTransformations> m_offDiagonalTrialTransformations;
    shared_ptr<TestKernelTrialIntegral> m_offDiagonalIntegral;
};

} // namespace Bempp

#endif