Evaluate a single component.
| 81 | |
| 82 | // Evaluate a single component. |
| 83 | double value(const libMesh::Point& p, const unsigned int component_n) const |
| 84 | { |
| 85 | TBOX_ASSERT(component_n < parsers.size()); |
| 86 | for (unsigned int var_n = 0; var_n < vars.size(); ++var_n) vars[var_n] = p(var_n); |
| 87 | |
| 88 | try |
| 89 | { |
| 90 | return parsers[component_n].Eval(); |
| 91 | } |
| 92 | catch (mu::ParserError& e) |
| 93 | { |
| 94 | std::cerr << "Message: <" << e.GetMsg() << ">\n"; |
| 95 | std::cerr << "Formula: <" << e.GetExpr() << ">\n"; |
| 96 | std::cerr << "Token: <" << e.GetToken() << ">\n"; |
| 97 | std::cerr << "Position: <" << e.GetPos() << ">\n"; |
| 98 | std::cerr << "Errc: <" << e.GetCode() << ">" << std::endl; |
| 99 | throw e; |
| 100 | } |
| 101 | return 1.0; |
| 102 | } |
| 103 | |
| 104 | // Evaluate all components. Only makes sense if the number of variables is |
| 105 | // equal to the number of functions. |