arts-docs-2.4/doxygen/interpolation__poly_8h_source.html

/* Copyright (C) 2012 Stefan Buehler <sbuehler(at)ltu.se>


   This program is free software; you can redistribute it and/or modify it

   under the terms of the GNU General Public License as published by the

   Free Software Foundation; either version 2, or (at your option) any

   later version.


   This program is distributed in the hope that it will be useful,

   but WITHOUT ANY WARRANTY; without even the implied warranty of

   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

   GNU General Public License for more details.


   You should have received a copy of the GNU General Public License

   along with this program; if not, write to the Free Software

   Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307,

   USA. */


#ifndef interpolation_poly_h

#define interpolation_poly_h


#include "interpolation.h"

#include "matpackI.h"


struct GridPosPoly {

  ArrayOfIndex idx;

  Vector w;

};


typedef Array<GridPosPoly> ArrayOfGridPosPoly;


ostream& operator<<(ostream& os, const GridPosPoly& gp);


void gridpos_poly(ArrayOfGridPosPoly& gp,

                  ConstVectorView old_grid,

                  ConstVectorView new_grid,

                  const Index order,

                  const Numeric& extpolfac = 0.5);


void gridpos_poly(GridPosPoly& gp,

                  ConstVectorView old_grid,

                  const Numeric& new_grid,

                  const Index order,

                  const Numeric& extpolfac = 0.5);


void gridpos_poly_longitudinal(const String& error_msg,

                               ArrayOfGridPosPoly& gp,

                               ConstVectorView old_grid,

                               ConstVectorView new_grid,

                               const Index order,

                               const Numeric& extpolfac = 0.5);


void gridpos_poly_cyclic_longitudinal(ArrayOfGridPosPoly& gp,

                                      ConstVectorView old_grid,

                                      ConstVectorView new_grid,

                                      const Index order,

                                      const Numeric& extpolfac = 0.5);


void gridpos_poly_longitudinal(GridPosPoly& gp,

                               ConstVectorView old_grid,

                               const Numeric& new_grid,

                               const Index order,

                               const Numeric& extpolfac = 0.5);


void gridpos_poly_cyclic_longitudinal(GridPosPoly& gp,

                                      ConstVectorView old_grid,

                                      const Numeric& new_grid,

                                      const Index order,

                                      const Numeric& extpolfac = 0.5);


//                      Red Interpolation


void interpweights(VectorView itw, const GridPosPoly& tc);


void interpweights(VectorView itw,

                   const GridPosPoly& tr,

                   const GridPosPoly& tc);


void interpweights(VectorView itw,

                   const GridPosPoly& tp,

                   const GridPosPoly& tr,

                   const GridPosPoly& tc);


void interpweights(VectorView itw,

                   const GridPosPoly& tb,

                   const GridPosPoly& tp,

                   const GridPosPoly& tr,

                   const GridPosPoly& tc);


void interpweights(VectorView itw,

                   const GridPosPoly& ts,

                   const GridPosPoly& tb,

                   const GridPosPoly& tp,

                   const GridPosPoly& tr,

                   const GridPosPoly& tc);


void interpweights(VectorView itw,

                   const GridPosPoly& tv,

                   const GridPosPoly& ts,

                   const GridPosPoly& tb,

                   const GridPosPoly& tp,

                   const GridPosPoly& tr,

                   const GridPosPoly& tc);


Numeric interp(ConstVectorView itw, ConstVectorView a, const GridPosPoly& tc);


Numeric interp(ConstVectorView itw,

               ConstMatrixView a,

               const GridPosPoly& tr,

               const GridPosPoly& tc);


Numeric interp(ConstVectorView itw,

               ConstTensor3View a,

               const GridPosPoly& tp,

               const GridPosPoly& tr,

               const GridPosPoly& tc);


Numeric interp(ConstVectorView itw,

               ConstTensor4View a,

               const GridPosPoly& tb,

               const GridPosPoly& tp,

               const GridPosPoly& tr,

               const GridPosPoly& tc);


Numeric interp(ConstVectorView itw,

               ConstTensor5View a,

               const GridPosPoly& ts,

               const GridPosPoly& tb,

               const GridPosPoly& tp,

               const GridPosPoly& tr,

               const GridPosPoly& tc);


Numeric interp(ConstVectorView itw,

               ConstTensor6View a,

               const GridPosPoly& tv,

               const GridPosPoly& ts,

               const GridPosPoly& tb,

               const GridPosPoly& tp,

               const GridPosPoly& tr,

               const GridPosPoly& tc);


//                      Blue interpolation


void interpweights(MatrixView itw, const ArrayOfGridPosPoly& cgp);


void interpweights(MatrixView itw,

                   const ArrayOfGridPosPoly& rgp,

                   const ArrayOfGridPosPoly& cgp);


void interpweights(MatrixView itw,

                   const ArrayOfGridPosPoly& pgp,

                   const ArrayOfGridPosPoly& rgp,

                   const ArrayOfGridPosPoly& cgp);


void interpweights(MatrixView itw,

                   const ArrayOfGridPosPoly& bgp,

                   const ArrayOfGridPosPoly& pgp,

                   const ArrayOfGridPosPoly& rgp,

                   const ArrayOfGridPosPoly& cgp);


void interpweights(MatrixView itw,

                   const ArrayOfGridPosPoly& sgp,

                   const ArrayOfGridPosPoly& bgp,

                   const ArrayOfGridPosPoly& pgp,

                   const ArrayOfGridPosPoly& rgp,

                   const ArrayOfGridPosPoly& cgp);


void interpweights(MatrixView itw,

                   const ArrayOfGridPosPoly& vgp,

                   const ArrayOfGridPosPoly& sgp,

                   const ArrayOfGridPosPoly& bgp,

                   const ArrayOfGridPosPoly& pgp,

                   const ArrayOfGridPosPoly& rgp,

                   const ArrayOfGridPosPoly& cgp);


void interp(VectorView ia,

            ConstMatrixView itw,

            ConstVectorView a,

            const ArrayOfGridPosPoly& cgp);


void interp(VectorView ia,

            ConstMatrixView itw,

            ConstMatrixView a,

            const ArrayOfGridPosPoly& rgp,

            const ArrayOfGridPosPoly& cgp);


void interp(VectorView ia,

            ConstMatrixView itw,

            ConstTensor3View a,

            const ArrayOfGridPosPoly& pgp,

            const ArrayOfGridPosPoly& rgp,

            const ArrayOfGridPosPoly& cgp);


void interp(VectorView ia,

            ConstMatrixView itw,

            ConstTensor4View a,

            const ArrayOfGridPosPoly& bgp,

            const ArrayOfGridPosPoly& pgp,

            const ArrayOfGridPosPoly& rgp,

            const ArrayOfGridPosPoly& cgp);


void interp(VectorView ia,

            ConstMatrixView itw,

            ConstTensor5View a,

            const ArrayOfGridPosPoly& sgp,

            const ArrayOfGridPosPoly& bgp,

            const ArrayOfGridPosPoly& pgp,

            const ArrayOfGridPosPoly& rgp,

            const ArrayOfGridPosPoly& cgp);


void interp(VectorView ia,

            ConstMatrixView itw,

            ConstTensor6View a,

            const ArrayOfGridPosPoly& vgp,

            const ArrayOfGridPosPoly& sgp,

            const ArrayOfGridPosPoly& bgp,

            const ArrayOfGridPosPoly& pgp,

            const ArrayOfGridPosPoly& rgp,

            const ArrayOfGridPosPoly& cgp);


//                      Green interpolation


void interpweights(Tensor3View itw,

                   const ArrayOfGridPosPoly& rgp,

                   const ArrayOfGridPosPoly& cgp);


void interpweights(Tensor4View itw,

                   const ArrayOfGridPosPoly& pgp,

                   const ArrayOfGridPosPoly& rgp,

                   const ArrayOfGridPosPoly& cgp);


void interpweights(Tensor5View itw,

                   const ArrayOfGridPosPoly& bgp,

                   const ArrayOfGridPosPoly& pgp,

                   const ArrayOfGridPosPoly& rgp,

                   const ArrayOfGridPosPoly& cgp);


void interpweights(Tensor6View itw,

                   const ArrayOfGridPosPoly& sgp,

                   const ArrayOfGridPosPoly& bgp,

                   const ArrayOfGridPosPoly& pgp,

                   const ArrayOfGridPosPoly& rgp,

                   const ArrayOfGridPosPoly& cgp);


void interpweights(Tensor7View itw,

                   const ArrayOfGridPosPoly& vgp,

                   const ArrayOfGridPosPoly& sgp,

                   const ArrayOfGridPosPoly& bgp,

                   const ArrayOfGridPosPoly& pgp,

                   const ArrayOfGridPosPoly& rgp,

                   const ArrayOfGridPosPoly& cgp);


void interp(MatrixView ia,

            ConstTensor3View itw,

            ConstMatrixView a,

            const ArrayOfGridPosPoly& rgp,

            const ArrayOfGridPosPoly& cgp);


void interp(Tensor3View ia,

            ConstTensor4View itw,

            ConstTensor3View a,

            const ArrayOfGridPosPoly& pgp,

            const ArrayOfGridPosPoly& rgp,

            const ArrayOfGridPosPoly& cgp);


void interp(Tensor4View ia,

            ConstTensor5View itw,

            ConstTensor4View a,

            const ArrayOfGridPosPoly& bgp,

            const ArrayOfGridPosPoly& pgp,

            const ArrayOfGridPosPoly& rgp,

            const ArrayOfGridPosPoly& cgp);


void interp(Tensor5View ia,

            ConstTensor6View itw,

            ConstTensor5View a,

            const ArrayOfGridPosPoly& sgp,

            const ArrayOfGridPosPoly& bgp,

            const ArrayOfGridPosPoly& pgp,

            const ArrayOfGridPosPoly& rgp,

            const ArrayOfGridPosPoly& cgp);


void interp(Tensor6View ia,

            ConstTensor7View itw,

            ConstTensor6View a,

            const ArrayOfGridPosPoly& vgp,

            const ArrayOfGridPosPoly& sgp,

            const ArrayOfGridPosPoly& bgp,

            const ArrayOfGridPosPoly& pgp,

            const ArrayOfGridPosPoly& rgp,

            const ArrayOfGridPosPoly& cgp);


#endif  // interpolation_poly_h