GreenPolyProducts/cantera
5
0
Fork 0
mirror of https://github.com/Cantera/cantera.git synced 2025-02-25 18:55:29 -06:00
Code Issues Packages Projects Releases Wiki Activity

*** empty log message ***

Browse Source
This commit is contained in:
Dave Goodwin
2007-04-25 00:14:36 +00:00
parent 92eb3d250e
commit 6221d28fe0
3 changed files with 250 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

240
Cantera/src/Func1.cpp Normal file
View File

@@ -0,0 +1,240 @@
#include "Func1.h"
#include "stringUtils.h"
using namespace std;
namespace Cantera {
Func1& Sin1::derivative() const {
Func1* c = new Cos1(m_c);
return *(new TimesConstant1(*c, m_c));
}
Func1& Cos1::derivative() const {
Func1* s = new Sin1(m_c);
return *(new TimesConstant1(*s, -m_c));
}
Func1& Exp1::derivative() const {
Func1* f = new Exp1(m_c);
if (m_c != 1.0)
return *(new TimesConstant1(*f, m_c));
else
return *f;
}
Func1& Pow1::derivative() const {
Func1* f = new Pow1(m_c - 1.0);
if (m_c != 1.0)
return *(new TimesConstant1(*f, m_c));
else
return *(new Const1(1.0));
}
string Func1::write(std::string arg) const {
return "<unknown " + int2str(ID()) + ">("+arg+")";
}
string Sin1::write(string arg) const {
string c = "";
if (m_c != 1.0) c = fp2str(m_c);
return "\\sin{"+c+arg+"}";
}
string Cos1::write(string arg) const {
string c = "";
if (m_c != 1.0) c = fp2str(m_c);
return "\\cos{"+c+arg+"}";
}
string Exp1::write(string arg) const {
string c = "";
if (m_c != 1.0) c = fp2str(m_c);
return "\\exp{"+c+arg+"}";
}
string Ratio1::write(string arg) const {
return "\\frac{" + m_f1->write(arg) + "}{"
+ m_f2->write(arg) + "}";
}
string Product1::write(string arg) const {
string s = m_f1->write(arg);
if (m_f1->order() < order()) s = "(" + s + ")";
string s2 = m_f2->write(arg);
if (m_f2->order() < order()) s2 = "(" + s2 + ")";
return s + " " + s2;
}
string Sum1::write(string arg) const {
string s1 = m_f1->write(arg);
string s2 = m_f2->write(arg);
if (s2[0] == '-') return s1 + " - " + s2.substr(1,s2.size());
else return s1 + " + " + s2;
}
string Diff1::write(string arg) const {
string s1 = m_f1->write(arg);
string s2 = m_f2->write(arg);
if (s2[0] == '-') return s1 + " + " + s2.substr(1,s2.size());
else return s1 + " - " + s2;
}
string Composite1::write(string arg) const {
string g = m_f2->write(arg);
return m_f1->write(g);
}
string TimesConstant1::write(string arg) const {
string s = m_f1->write(arg);
if (m_f1->order() < order()) s = "(" + s + ")";
if (m_c == 1.0) return s;
if (m_c == -1.0) return "-"+s;
return fp2str(m_c) + s;
}
string PlusConstant1::write(string arg) const {
if (m_c == 0.0) return m_f1->write(arg);
return m_f1->write(arg) + " + " + fp2str(m_c);
}
doublereal Func1::isProportional(TimesConstant1& other) {
if (isIdentical(other.func1())) return other.c();
return 0.0;
}
doublereal Func1::isProportional(Func1& other) {
if (isIdentical(other)) return 1.0;
else return 0.0;
}
static bool isZero(Func1& f) {
if (f.ID() == ConstFuncType && f.c() == 0.0)
return true;
else
return false;
}
static bool isOne(Func1& f) {
if (f.ID() == ConstFuncType && f.c() == 1.0)
return true;
else
return false;
}
static bool isTimesConst(Func1& f) {
if (f.ID() == TimesConstantFuncType)
return true;
else
return false;
}
Func1& newSumFunction(Func1& f1, Func1& f2) {
if (isZero(f1)) {
delete &f1;
return f2;
}
if (isZero(f2)) {
delete &f2;
return f1;
}
doublereal c = f1.isProportional(f2);
if (c != 0) {
if (c == -1.0) return *(new Const1(0.0));
else return newTimesConstFunction(f1, c + 1.0);
}
return *(new Sum1(f1, f2));
}
Func1& newDiffFunction(Func1& f1, Func1& f2) {
if (isZero(f2)) {
delete &f2; return f1;
}
if (f1.isIdentical(f2)) {
delete &f1; delete &f2;
return *(new Const1(0.0));
}
doublereal c = f1.isProportional(f2);
if (c != 0) {
if (c == 1.0) return *(new Const1(0.0));
else return newTimesConstFunction(f1, 1.0 - c);
}
return *(new Diff1(f1, f2));
}
Func1& newProdFunction(Func1& f1, Func1& f2) {
if (isOne(f1)) {
delete &f1; return f2;
}
if (isOne(f2)) {
delete &f2; return f1;
}
if (isZero(f1) || isZero(f2)) {
delete &f1; delete &f2;
return *(new Const1(0.0));
}
bool tc1 = isTimesConst(f1);
bool tc2 = isTimesConst(f2);
//cout << "product: " << f1.write("t") << " " << f2.write("t") << endl;
//cout << tc1 << " " << tc2 << " " << f1.ID() << " " << f2.ID() << endl;
if (tc1 || tc2) {
doublereal c1 = 1.0, c2 = 1.0;
Func1 *ff1 = 0, *ff2 = 0;
if (tc1) {
c1 = f1.c();
ff1 = &f1.func1_dup();
}
else ff1 = &f1;
if (tc2) {
c2 = f2.c();
ff2 = &f2.func1_dup();
}
else ff2 = &f2;
Func1& p = newProdFunction(*ff1, *ff2);
cout << "p = " << p.write("t") << endl;
if (c1*c2 != 1.0) {
return newTimesConstFunction(p, c1*c2);
}
else
return p;
}
else
return *(new Product1(f1, f2));
}
Func1& newRatioFunction(Func1& f1, Func1& f2) {
if (isOne(f2)) return f1;
if (isZero(f1)) return *(new Const1(0.0));
if (f1.isIdentical(f2)) {
delete &f1; delete &f2;
return *(new Const1(1.0));
}
return *(new Ratio1(f1, f2));
}
Func1& newCompositeFunction(Func1& f1, Func1& f2) {
if (isZero(f1)) {
delete &f1; delete &f2;
return *(new Const1(0.0));
}
return *(new Composite1(f1, f2));
}
Func1& newTimesConstFunction(Func1& f, doublereal c) {
if (c == 0.0) {
delete &f;
return *(new Const1(0.0));
}
if (c == 1.0) {
return f;
}
return *(new TimesConstant1(f, c));
}
Func1& newPlusConstFunction(Func1& f, doublereal c) {
if (c == 0.0) {
return f;
}
return *(new PlusConstant1(f, c));
}
}

5
ext/f2c_math/printstring.c Normal file
View File

@@ -0,0 +1,5 @@
#include <stdio.h>
void printstring_(char* s) {
printf("%s",s);
}

5
ext/math/printstring.c Normal file
View File

@@ -0,0 +1,5 @@
#include <stdio.h>
void printstring_(char* s) {
printf("%s",s);
}