src/math/ChebyshevPolynomials.cpp

/*
 * Copyright (C) 2000-2004  Object Oriented Parallel Simulation Engine (OOPSE) project
 * 
 * Contact: oopse@oopse.org
 * 
 * This program is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public License
 * as published by the Free Software Foundation; either version 2.1
 * of the License, or (at your option) any later version.
 * All we ask is that proper credit is given for our work, which includes
 * - but is not limited to - adding the above copyright notice to the beginning
 * of your source code files, and to any copyright notice that you may distribute
 * with programs based on this work.
 * 
 * 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 Lesser General Public License for more details.
 * 
 * You should have received a copy of the GNU Lesser 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 "math/ChebyshevPolynomials.hpp"

namespace oopse {
ChebyshevPolynomials::ChebyshevPolynomials(int maxPower) : maxPower_(maxPower){

    assert(maxPower >= 0);
    GeneratePolynomials(maxPower_);
}

ChebyshevPolynomials::GeneratePolynomials(int maxPower) {
    DoublePolynomial t0;
    t0.setCoefficient(0, 1.0);
    polyList_.push_back(t0);
    
    DoublePolynomial t1;
    t1.setCoefficient(1, 1.0);

    DoublePolynomial twoX;
    twoX.setCoefficient(1, 2.0);
    
    for (int i = 2; i <= maxPower; ++i) {
        DoublePolynomial tn;
        
        tn = polyList_[i] * twoX - polyList_[i-2];
        polyList_.push_back(tn);
    }
}

} //end namespace oopse
Revision:	1747
Committed:	Wed Nov 17 18:58:49 2004 UTC (19 years, 7 months ago) by tim
File size:	1821 byte(s)
Log Message:	more types stuff and Chebyshev Polynomial
#	User	Rev	Content
1	tim	1747	/*
2			* Copyright (C) 2000-2004 Object Oriented Parallel Simulation Engine (OOPSE) project
3			*
4			* Contact: oopse@oopse.org
5			*
6			* This program is free software; you can redistribute it and/or
7			* modify it under the terms of the GNU Lesser General Public License
8			* as published by the Free Software Foundation; either version 2.1
9			* of the License, or (at your option) any later version.
10			* All we ask is that proper credit is given for our work, which includes
11			* - but is not limited to - adding the above copyright notice to the beginning
12			* of your source code files, and to any copyright notice that you may distribute
13			* with programs based on this work.
14			*
15			* This program is distributed in the hope that it will be useful,
16			* but WITHOUT ANY WARRANTY; without even the implied warranty of
17			* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18			* GNU Lesser General Public License for more details.
19			*
20			* You should have received a copy of the GNU Lesser General Public License
21			* along with this program; if not, write to the Free Software
22			* Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
23			*
24			*/
25
26			#include "math/ChebyshevPolynomials.hpp"
27
28			namespace oopse {
29			ChebyshevPolynomials::ChebyshevPolynomials(int maxPower) : maxPower_(maxPower){
30
31			assert(maxPower >= 0);
32			GeneratePolynomials(maxPower_);
33			}
34
35			ChebyshevPolynomials::GeneratePolynomials(int maxPower) {
36			DoublePolynomial t0;
37			t0.setCoefficient(0, 1.0);
38			polyList_.push_back(t0);
39
40			DoublePolynomial t1;
41			t1.setCoefficient(1, 1.0);
42
43			DoublePolynomial twoX;
44			twoX.setCoefficient(1, 2.0);
45
46			for (int i = 2; i <= maxPower; ++i) {
47			DoublePolynomial tn;
48
49			tn = polyList_[i] * twoX - polyList_[i-2];
50			polyList_.push_back(tn);
51			}
52			}
53
54			} //end namespace oopse