/*
 * YieldPanel.cpp
 *
 * Takes few pairs of (maturity, forward rates) from the term structure
 * of zero coupon bonds, and performs a fit using the Spline class
 */

#include "Curve.cpp"
#include "Spline.cpp"
#include <vector>
#include <iostream>
#include <fstream>
#include <string>
#include <math.h>
/**
 *
 * @author  horacio aliaga
 */
class YieldPanel {
    vector<double> t;
    vector<double> y;
    double* time;
    double* yield;
    int npoints;    
    /** Creates new form DataTextField */
public: YieldPanel() {
        t.push_back(0.0);t.push_back(0.2);t.push_back(0.4);
		t.push_back(0.8);t.push_back(1.0);
        y.push_back(0.10);y.push_back(0.14);y.push_back(0.24);
		y.push_back(0.15);y.push_back(0.11);
        Spline s;
        Curve data;
        data=s.spline(t,y,400,1.0);
    }
public: void set_npoints(int _npoints){
        npoints = _npoints;
    }
public: int get_npoints(){
        return npoints;
    }    
public:	void Go() {                                           
        Spline s;
        Curve data;
        data=s.spline(t,y,400,1.0);
        return;
    }                                          

public: void Add(double t1, double y1) {                                            
        t.push_back(t1);
        y.push_back(y1);

    }                                           

public: void ClearAll() {
        Curve data;
        if(t.size()>0){
            t.clear();
            y.clear();
        }
    } 
public: void set_time(double maturity){
        double dt = maturity/(1.0*npoints);
        double* time=new double[npoints];
        double t=0;
        for (int i=0;i<npoints;i++){
            time[i] = t;
            t = t + dt;
        }
        return;
    }
public: void set_yield(double maturity){
        double* yield=new double[npoints];
        Spline s;
        yield = s.spline_array(get_t(),get_y(),get_npoints(),maturity);
        return;
    }    
public: double* get_time(){
        return time;
    }
public: double* get_yield(){
        return yield;
    }
public: vector<double> get_t(){
        return t;
    }
public: vector<double> get_y(){
        return y;
    }
};