codac2_Parallelepiped_eval.h

Go to the documentation of this file.

 
#pragma once
 
#include "codac2_AnalyticFunction.h"
#include "codac2_Parallelepiped.h"
#include "codac2_peibos.h"
 
namespace codac2
{
  template<typename T>
    requires std::is_base_of_v<AnalyticTypeBase,T>
  template<typename... Args>
  inline Parallelepiped AnalyticFunction<T>::parallelepiped_eval(const Args&... x) const
  {
    this->check_valid_inputs(x...);
    assert_release(this->input_size() < this->output_size() &&
                   "Parallelepiped evaluation requires more outputs than inputs.");
    assert_release(this->input_size() > 0 &&
                "Parallelepiped evaluation requires at least one input.");
 
    auto Y = this->eval(x.mid()...);
    auto z = Y.mid();
 
    Matrix A = this->diff(x.mid()...).mid();
 
    // Maximum error computation
    double rho = error(Y, z, this->diff(x...), A, cart_prod(x...));
 
    // Inflation of the parallelepiped
    auto A_inf = inflate_flat_parallelepiped(A, cart_prod(x...).rad(), rho);
 
    return Parallelepiped(z, A_inf);
  }
}