My Project
2d/ppmatrix.hh
Go to the documentation of this file.
1 /* -*- mia-c++ -*-
2  *
3  * This file is part of MIA - a toolbox for medical image analysis
4  * Copyright (c) Leipzig, Madrid 1999-2017 Gert Wollny
5  *
6  * MIA is free software; you can redistribute it and/or modify
7  * it under the terms of the GNU General Public License as published by
8  * the Free Software Foundation; either version 3 of the License, or
9  * (at your option) any later version.
10  *
11  * This program is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  * GNU General Public License for more details.
15  *
16  * You should have received a copy of the GNU General Public License
17  * along with MIA; if not, see <http://www.gnu.org/licenses/>.
18  *
19  */
20 
21 #ifndef mia_2d_ppmatrix_hh
22 #define mia_2d_ppmatrix_hh
23 
24 #include <mia/2d/interpolator.hh>
25 #include <mia/2d/datafield.hh>
26 #include <mia/core/vector.hh>
27 
29 
41 {
42 public:
53  C2DPPDivcurlMatrix(const C2DBounds& size, const C2DFVector& range, const CSplineKernel& kernel,
54  double wd, double wr);
55 
57 
63  double operator * (const C2DFVectorfield& coefficients) const;
64 
70  double operator * (const T2DDatafield<C2DDVector>& coefficients) const;
71 
72 
80  double evaluate(const C2DFVectorfield& coefficients, CDoubleVector& gradient) const;
81 
89  double evaluate(const T2DDatafield<C2DDVector>& coefficients, CDoubleVector& gradient) const;
90 
94  const C2DBounds& get_size() const;
95 
99  void reset(const C2DBounds& size, const C2DFVector& range, const CSplineKernel& kernel,
100  double wd, double wr);
101 private:
102  class C2DPPDivcurlMatrixImpl *impl;
103 
104 };
106 
107 #endif
a 2D field of floating point single accuracy 2D vectors
Precomputed matrix for the DivCurl regularization.
Definition: 2d/ppmatrix.hh:41
const C2DBounds & get_size() const
C2DPPDivcurlMatrix(const C2DBounds &size, const C2DFVector &range, const CSplineKernel &kernel, double wd, double wr)
double evaluate(const C2DFVectorfield &coefficients, CDoubleVector &gradient) const
void reset(const C2DBounds &size, const C2DFVector &range, const CSplineKernel &kernel, double wd, double wr)
double evaluate(const T2DDatafield< C2DDVector > &coefficients, CDoubleVector &gradient) const
Base class for all spline based interpolation kernels.
Definition: splinekernel.hh:46
A class to hold data on a regular 2D grid.
Definition: 2d/datafield.hh:90
A wrapper around the c-array to provide an STL like interface for iterators.
Definition: core/vector.hh:78
#define EXPORT_2D
Definition: defines2d.hh:37
#define NS_MIA_BEGIN
conveniance define to start the mia namespace
Definition: defines.hh:33
#define NS_MIA_END
conveniance define to end the mia namespace
Definition: defines.hh:36
Matrix EXPORT_GSL operator*(const Matrix &lhs, const Matrix &rhs)