test_cia.cc

Go to the documentation of this file.

/* Copyright (C) 2012 Stefan Buehler <sbuehler@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. */

#include "cia.h"
#include "make_vector.h"

void test01()
{
    cout << "Testing CIA Interpolation.\n";
    GriddedField2 cia_data;
    
    Matrix A(5,3,0.);
    A(2,1) = 1;
    //    cout << "A:" << A << endl;
    
    cia_data.data = A;
    cia_data.set_grid(0, MakeVector(1,2,3,4,5));
    cia_data.set_grid(1, MakeVector(100,200,300));
    
    cout << "cia_data:" << cia_data << endl;
    
    
    // Output frequencies and temperature:
    Vector f_out(1,9,0.5);
    cout << "f_out:" << f_out << endl;
    Numeric T_out = 150;
    cout << "T_out:" << T_out << endl;

    Vector result(9);
    cia_interpolation(result,
                      f_out,
                      T_out,
                      cia_data,
                      0.5,
                      0,
                      Verbosity(3,3,0));
    cout << "result:" << result << endl;

}

int main()
{
  test01();
  return 0;
}