recsAccuracy Evaluates approximation accuracy
SE = recsAccuracy(MODEL,INTERP,S) evaluates the accuracy of the approximation
defined in the interpolition structure INTERP for the model defined in the
object MODEL. The accuracy is assessed over the state variables contained
in the n-by-d-by-y array S, as output by resSimul. recsAccuracy returns a
n*y-by-d matrix SE containing the state variables on which the accuracy was
evaluated.
MODEL is an object created by recsmodel.
INTERP is a structure, which includes the following fields:
cx : coefficient matrix of the interpolation of the response variables
fspace : a definition structure for the interpolation family (created by
the function fundef)
SE = recsAccuracy(MODEL,INTERP,S,OPTIONS) evaluates the accuracy with the
parameters defined by the structure OPTIONS. The fields of the
structure are
extrapolate : 1 if extrapolation is allowed outside the
interpolation space or 0 to forbid it (default: 1)
[SE,LEE] = recsAccuracy(MODEL,INTERP,S,...) returns the matrix LEE containing
the value of Euler equation error (in log10) evaluated on the grid points SE.
[SE,LEE,lEF] = recsAccuracy(MODEL,INTERP,S,...) returns the nxy-by-m matrix
lEF containing the value of equilibrium equation error (in log10) evaluated on
the grid points SE. For models with finite bounds, this error is assessed
using a minmax formulation: Ef = abs(min(max(-f,LB-x),UB-x)).