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