#include <codi.hpp>
#include <iostream>
#include <omp.h>
template<typename Real>
void func(
const Real* x,
size_t l,
Real* y) {
y[0] = 0.0;
y[1] = 1.0;
for(size_t i = 0; i < l; ++i) {
y[0] += x[i];
y[1] *= x[i];
}
}
int main(int nargs, char** args) {
x[0] = 1.0;
x[1] = 2.0;
x[2] = 3.0;
x[3] = 4.0;
x[4] = 5.0;
tape.setActive();
for(size_t i = 0; i < 5; ++i) {
tape.registerInput(x[i]);
}
func(x, 5, y);
tape.registerOutput(y[0]);
tape.registerOutput(y[1]);
tape.setPassive();
#pragma omp parallel num_threads(2)
{
int tid = omp_get_thread_num();
for(size_t i = 0; i < 5; ++i) {
}
}
std::cout << "Custom adjoint vector helper:" << std::endl;
std::cout << "f(1 .. 5) = (" << y[0] << ", " << y[1] << ")" << std::endl;
std::cout << "df/dx (1 .. 5) = \n" << jacobian << std::endl;
tape.reset();
}
DataExtraction< Type >::Identifier getIdentifier(Type const &v)
Extract the identifiers from a type of aggregated active types.
Definition: realTraits.hpp:216
RealReverseGen< double > RealReverse
Definition: codi.hpp:120
static Tape & getTape()
Get a reference to the tape which manages this expression.
Definition: activeType.hpp:99
T_Tape Tape
See ActiveType.
Definition: activeType.hpp:55
Allows for an arbitrary adjoint evaluation of a recorded tape.
Definition: customAdjointVectorHelper.hpp:152
Gradient & gradient(Identifier const &identifier)
Get a reference to the gradient. Checked access.
Definition: customAdjointVectorHelper.hpp:254
Gradient const & getGradient(Identifier const &identifier) const
Get a constant reference to the gradient.
Definition: customAdjointVectorHelper.hpp:239
void evaluate(Position const &start, Position const &end)
Perform a full reverse evaluation of the tape.
Definition: customAdjointVectorHelper.hpp:207
Default implementation of the Jacobian interface.
Definition: jacobian.hpp:60