formula_parser.cpp

// Copyright (C) 2005, Ondra Kamenik

// $Id: formula_parser.cpp 2268 2008-11-22 10:38:03Z michel $

#include "utils/cc/pascal_triangle.h"
#include "utils/cc/exception.h"

#include "parser_exception.h"
#include "location.h"
#include "formula_parser.h"
#include "formula_tab.hh"

#include <cmath>

using namespace ogp;

extern location_type fmla_lloc;

FormulaParser::FormulaParser(const FormulaParser& fp, Atoms& a)
	: otree(fp.otree), atoms(a), formulas(fp.formulas), ders()
{
	// create derivatives
	for (unsigned int i = 0; i < fp.ders.size(); i++)
		ders.push_back(new FormulaDerivatives(*(fp.ders[i])));
}

FormulaParser::~FormulaParser()
{
	destroy_derivatives();
}

void FormulaParser::differentiate(int max_order)
{
	destroy_derivatives();
	vector<int> vars;
	vars = atoms.variables();
	for (unsigned int i = 0; i < formulas.size(); i++)
		ders.push_back(new FormulaDerivatives(otree, vars, formulas[i], max_order));
}

const FormulaDerivatives& FormulaParser::derivatives(int i) const
{
	if (i < (int)ders.size())
		return *(ders[i]);
	else
		throw ogu::Exception(__FILE__,__LINE__,
							 "Wrong formula index in FormulaParser::derivatives");
	return *(ders[0]); // just because of compiler
}



void FormulaParser::add_formula(int t)
{
	formulas.push_back(t);
}

int FormulaParser::add_binary(code_t code, int t1, int t2)
{
	return otree.add_binary(code, t1, t2);
}

int FormulaParser::add_unary(code_t code, int t)
{
	return otree.add_unary(code, t);
}

int FormulaParser::add_nulary(const char* str)
{
	int t = -1;