quadpts.h

Go to the documentation of this file.

//

//--------------------------------------------------------------------------
//--------------------------------------------------------------------------
//  
//  libMoM Library - a library for stochastic simulations in engineering.
// 
//  Copyright (C) 2010 Rochan R. Upadhyay
// 
//  This library is free software; you can redistribute it and/or
//  modify it under the terms of the Version 2.1 GNU Lesser General
//  Public License as published by the Free Software Foundation.
// 
//  This library 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
//  Lesser General Public License for more details.
// 
//  You should have received a copy of the GNU Lesser General Public
//  License along with this library; if not, write to the Free Software
//  Foundation, Inc. 51 Franklin Street, Fifth Floor, 
//  Boston, MA  02110-1301  USA
// 
//--------------------------------------------------------------------------

#ifndef _QUADPTS_H
#define _QUADPTS_H

#include <stdio.h>
#include <stdlib.h>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_odeiv.h>
#include "gsl_vector_integer.h"
#include "gsl_vector_double.h"
#include "gsl_matrix_double.h"
#include "gsl_matrix_integer.h"
#include "gsl_vector_matrix_ops.h"
#include "gsl_ode_solver.h"
#include "moments.h"

class moments;

//

class quadpts : public vector
{
 
public:

//

    quadpts();

// 
//

    quadpts ( size_t a );

//
//
//

    moments getmoments ( size_t a, size_t b ) const;

//
//
//

    moments get_fixed_number_moments_1D ( size_t num_moments, size_t num_points ) const;

//
//
//

    moments get_user_defined_moments_GLD_1D ( size_t dimension, size_t num_points ) const;
    
     moments get_user_defined_moments_1D ( size_t dimension, size_t num_points, size_t number ) const;

//
//
    
    matrix get_non_singular_Jacobian ( size_t a, size_t b ) const;

//
//

    vector get_source ( size_t dimension, size_t num_points ) const;

//
//
    
    void solve ( size_t dimension, size_t num_points, double dt,
                                    long int num_steps, int print_interval );

//
//
//

    void solve_with_gsl ( size_t dimension, size_t num_points, double dt,
                                    long int num_steps, int print_interval );                                

//
//
    
    void post_process_quadpts (size_t dimension, size_t num_points, double time) const;
    
//
//

    matrix get_cdf_empirical (size_t dimension, size_t num_points, size_t num_intervals) const;

//
//

    void return_quadwts_quadpts_vector (size_t dimension, size_t num_points, double y[]) const;

private:

//
    vector quads;
    
};
  

class my_ode_function : public ode_function
{
   public:   

   vector operator()(double time, const vector& x);   

};

#endif // QUADPTS_H