laplacian2.h

Go to the documentation of this file.

// ========================================================================= //
// This file is part of MyraMath, copyright (c) 2014-2019 by Ryan A Chilton  //
// and distributed by MyraCore, LLC. See LICENSE.txt for license terms.      //
// ========================================================================= //

#ifndef MYRAMATH_SPARSE_LAPLACIAN2_H
#define MYRAMATH_SPARSE_LAPLACIAN2_H

#include <myramath/utility/detail/LIBPUBLIC.h>
#include <myramath/utility/Number.h>

#include <utility>
#include <vector>

namespace myra {

// Forward declarations.
class Pattern;
class Permutation;
template<class Number> class SparseMatrix;

class LIBPUBLIC Natural2D
  {

  public:

    Natural2D(int in_I, int in_J);

    std::vector<int> size() const;
    
    int operator() (int i, int j) const;

    std::pair<Natural2D, Natural2D> i_split() const;

    std::pair<Natural2D, Natural2D> j_split() const;
    
  private:

    // Full constructor used by splitting methods.
    Natural2D(int in_I, int in_i_begin, int in_i_end, int in_J, int in_j_begin, int in_j_end);

    // Private data - total size (I/J) and begin/end specification of subwindow.
    int I, i_begin, i_end;
    int J, j_begin, j_end;
  };

LIBPUBLIC Pattern stencil2(int I, int J);

template<class Number> SparseMatrix<Number> laplacian2 (int I, int J);
template<> LIBPUBLIC SparseMatrix<NumberS> laplacian2<NumberS>(int I, int J);
template<> LIBPUBLIC SparseMatrix<NumberD> laplacian2<NumberD>(int I, int J);
template<> LIBPUBLIC SparseMatrix<NumberC> laplacian2<NumberC>(int I, int J);
template<> LIBPUBLIC SparseMatrix<NumberZ> laplacian2<NumberZ>(int I, int J);

LIBPUBLIC Permutation bisect2 (int I, int J);

LIBPUBLIC Pattern stencil2_unsymmetric(int I, int J, double density = 0.5);

} // namespace myra

#endif