log_expr.cc

由matlab开发的hybrid系统的描述语言

/*
	HYSDEL
	Copyright (C) 1999-2002  Fabio D. Torrisi

	This file is part of HYSDEL.
    
	HYSDEL is free software; you can redistribute it and/or
	modify it under the terms of the GNU General Public
	License as published by the Free Software Foundation; either
	version 2 of the License, or (at your option) any later version.

	HYSDEL 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
	General Public License for more details.

	You should have received a copy of the GNU General Public
	License along with this library; if not, write to the Free Software
	Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA

	CONTACT INFORMATION
	===================
	Fabio D. Torrisi
	ETH Zentrum
	Physikstrasse. 3 ETL,
	CH-8032 Zurich
	Switzerland
	mailto:[e-mail] (preferred)
*/

#include "Log_expr.h"
#include "Affine_func.h"
#include "Real_number.h"
#include <stdio.h>
#include "Globals.h"
#include "Problem.h"
#include "Problem_handler.h"

Log_expr::Log_expr(Expr * c) : Unary_expr(c) {
}

void Log_expr::semantic_checks() {
	child->semantic_checks();
	assert(child->is_real()); //fine
}


Affine_func * Log_expr::compute_affine() const {
	Affine_func * aff;

	aff = new Affine_func(clone(), globals); // put this as constant
	return aff;
}

bool Log_expr::is_affine() const {
	return child->is_const();
}

bool Log_expr::is_number() const {
	return child->is_number();
}

Number * Log_expr::compute_number() const {
	Real_number * num;

	num = (Real_number *) child->compute_number();
	num->log();
	return num;
}

string Log_expr::arg_range_check_matlab() const {
	string res;
	string msg;
	char buf[100];

	res += Unary_expr::arg_range_check_matlab();
	if (child->is_number()) {
		if (((Real_number *) child->compute_number())->get_double() <= 0.0) {
			sprintf(buf, "line %d: ", get_source_line());
			msg = buf;
			msg += string("log of non positive (") +
				child->get_source() + string(")");
			globals->problem_handler->process(
				new Problem(ERROR, msg));
		}
	} else if (child->is_const()) {
		res += string("if ") + child->to_matlab() + string(" <= 0.0") +
			string(" then error(\'");
		sprintf(buf, "line %d: ", get_source_line());
		res += buf;
		res += string("log of non positive (") +
			child->get_source() + string(")");
		res += string("\'); end\n");
	}
	return res;
}