advantage 0.1.0

Automatic Differentiation Toolkit
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36

#ifndef _ADV_DRIVERS_HPP
#define _ADV_DRIVERS_HPP

#include "Tape.hpp"

#include <vector>

namespace adv
{

std::vector<double> zero_order(const Tape& tape, const std::vector<double>& x);
std::pair<std::vector<double>, std::vector<double>> first_order(const Tape& tape, const std::vector<double>& x, const std::vector<double>& dx);
std::pair<std::vector<double>, std::vector<double>> first_order_reverse(const Tape& tape, const std::vector<double>& x, const std::vector<double>& ybar);

std::vector<double> jacobian(const Tape& tape, const std::vector<double>& x);
std::vector<double> jacobian_reverse(const Tape& tape, const std::vector<double>& x);

struct AbsNormalForm {
	std::size_t n;
	std::size_t m;
	std::size_t s;

	std::vector<double> a;
	std::vector<double> Z;
	std::vector<double> L;

	std::vector<double> b;
	std::vector<double> J;
	std::vector<double> Y;
};

AbsNormalForm abs_normal(const Tape& tape, const std::vector<double>& x);

} // namespace adv

#endif // _ADV_DRIVERS_HPP