//
//  $filename
//  $projectname
//
//  Created by Joshua moerman on $TODAY.
//  Copyright 2010 Joshua Moerman. All rights reserved.
//

#include "PolynomialA3D.hpp"
#include <algorithm>

PolynomialA3D::PolynomialA3D():
		KernelBasicParameter(3),
		KernelBasicVector(3){
    init();
}

void PolynomialA3D::init() {

}

void PolynomialA3D::operator()() {
	std::swap(vectorNew, vectorOld);

    vectorNew[0] = parameters[0] + vectorOld[1] - vectorOld[1]*vectorOld[2];
    vectorNew[1] = parameters[1] + vectorOld[2] - vectorOld[2]*vectorOld[0];
    vectorNew[2] = parameters[2] + vectorOld[0] - vectorOld[0]*vectorOld[1];
}
/*
N dimensional:
    for ( unsigned int i = 0; i < dim; i++ ) {
        new_point[i] = param[i] + point[(i+1) % dim] - point[(i+1) % dim]*point[(i+2) % dim];
    }
*/