LMU/musrfit
LMU/musrfit
5
0
Fork 0
Code Issues 6 Pull Requests Actions Packages Projects Releases 3 Activity

preventing more nan's by replacing func(a) -> func(|a|) for func: sqrt(), ln(), log()

Browse Source
This commit is contained in:
nemu 2010-03-09 10:13:24 +00:00
parent 78d0469993
commit 901f12b5c5
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/classes/PFunction.cpp
View File

@ -458,13 +458,13 @@ Double_t PFunction::EvalNode(PFuncTreeNode &node)
} else if (node.fFunctionTag == FUN_ATANH) { } else if (node.fFunctionTag == FUN_ATANH) {
return atanh(EvalNode(node.children[0])); return atanh(EvalNode(node.children[0]));
} else if (node.fFunctionTag == FUN_LOG) { } else if (node.fFunctionTag == FUN_LOG) {
return log(EvalNode(node.children[0]))/log(10.0); return log(fabs(EvalNode(node.children[0])))/log(10.0);
} else if (node.fFunctionTag == FUN_LN) { } else if (node.fFunctionTag == FUN_LN) {
return log(EvalNode(node.children[0])); return log(fabs(EvalNode(node.children[0])));
} else if (node.fFunctionTag == FUN_EXP) { } else if (node.fFunctionTag == FUN_EXP) {
return exp(EvalNode(node.children[0])); return exp(EvalNode(node.children[0]));
} else if (node.fFunctionTag == FUN_SQRT) { } else if (node.fFunctionTag == FUN_SQRT) {
return sqrt(EvalNode(node.children[0])); return sqrt(fabs(EvalNode(node.children[0])));
} else { } else {
cerr << endl << "**PANIC ERROR**: PFunction::EvalNode: node.fID == PFunctionGrammar::functionID: you never should have reached this point!"; cerr << endl << "**PANIC ERROR**: PFunction::EvalNode: node.fID == PFunctionGrammar::functionID: you never should have reached this point!";
cerr << endl; cerr << endl;