space.cpp

/*
    Copyright 2017 Philippe Grandclement
 
    This file is part of Kadath.
 
    Kadath 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 3 of the License, or
    (at your option) any later version.
 
    Kadath 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 Kadath.  If not, see <http://www.gnu.org/licenses/>.
*/
 
#include "headcpp.hpp"
#include "space.hpp"
#include <assert.h>
namespace Kadath {
double eta_lim_chi (double chi, double rext, double a, double eta_c) ;
double chi_lim_eta (double chi, double rext, double a, double chi_c) ;
 
Space::Space () : nbr_domains(0), ndim(0), type_base(0), domains(0x0) {
}
 
 
Space::~Space () {
   if (domains != 0x0) {
       for (int i=0 ; i<nbr_domains ; i++)
         delete domains[i] ;
       delete [] domains ;
    }
}
 
ostream& operator<< (ostream& o, const Space& so) {
    o << "Space of " << so.nbr_domains << " domains" << endl ;
    for (int i=0 ; i<so.nbr_domains ; i++) {
       o << "Domain " << i+1 << endl ;
      o << *so.domains[i] << endl ;
       o << "********************" << endl ;
    }
    switch (so.type_base) {
    case CHEB_TYPE :
        o << "Chebyshev polynomials are used." << endl ;
        break ;
    case LEG_TYPE :
        o << "Legendre polynomials are used." << endl ;
        break ;
    default :
        cerr << "Unknown type of base" << endl ;
        abort() ;
    }
    return o ;
}
 
void Space::save (FILE*) const {
    cerr << "save not implemented for" << endl ;
    cerr << *this << endl ;
    abort() ;
}
 
Array<int> Space::get_indices_matching_non_std(int, int) const {
    cerr << "get_indices_matching_non_std not implemented for" << endl ;
    cerr << *this << endl ;
    abort() ;
}
}