Function getErrorNorms

test/testingutilities.cpp:785–854

}

DenseMatrix getErrorNorms(const Function *exact, const Function *approx, const std::vector<std::vector<double>> &points)
{
    assert(exact->getNumVariables() == approx->getNumVariables());

    DenseVector normOneValVec(points.size());
    DenseVector normTwoValVec(points.size());
    DenseVector normInfValVec(points.size());

    DenseVector normOneJacVec(points.size());
    DenseVector normTwoJacVec(points.size());
    DenseVector normInfJacVec(points.size());

    DenseVector normOneHesVec(points.size());
    DenseVector normTwoHesVec(points.size());
    DenseVector normInfHesVec(points.size());

    int i = 0;
    for (auto &point : points) {
        DenseVector x = vecToDense(point);

        {
            DenseMatrix exactValue(1,1);
            exactValue(0,0) = exact->eval(x);
            DenseMatrix approxValue(1,1);
            approxValue(0,0) = approx->eval(x);
            DenseMatrix error = exactValue - approxValue;

            normOneValVec(i) = getOneNorm(error);
            normTwoValVec(i) = getTwoNorm(error);
            normInfValVec(i) = getInfNorm(error);
        }

        {
            auto exactJacobian = exact->evalJacobian(x);
            auto approxJacobian = approx->evalJacobian(x);
            auto errorJacobian = exactJacobian - approxJacobian;

            normOneJacVec(i) = getOneNorm(errorJacobian);
            normTwoJacVec(i) = getTwoNorm(errorJacobian);
            normInfJacVec(i) = getInfNorm(errorJacobian);
        }

        {
            auto exactHessian = exact->evalHessian(x);
            auto approxHessian = approx->evalHessian(x);
            auto errorHessian = exactHessian - approxHessian;

            normOneHesVec(i) = getOneNorm(errorHessian);
            normTwoHesVec(i) = getTwoNorm(errorHessian);
            normInfHesVec(i) = getInfNorm(errorHessian);
        }

        i++;
    }

    DenseMatrix errorNorms(3,3);
    errorNorms(0,0) = getOneNorm(normOneValVec);
    errorNorms(0,1) = getTwoNorm(normTwoValVec);