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
37
38
39
40
41
42
43
44
45
|
#pragma once
#include "dynmlp.h"
#include "ode_solver.h"
#include "utils.h"
typedef struct {
DynMLP net;
int nparams;
double trace_eps; /* finite-difference epsilon for trace computation */
} CNF;
typedef struct {
double *z1;
double delta_logp;
int nfe;
} CNFSampleResult;
typedef struct {
double *z0;
double delta_logp;
int nfe;
} CNFLogProbResult;
typedef struct {
double *dL_dz0;
double *dL_dtheta;
int nfe;
} CNFBackwardResult;
void cnf_init(CNF *cnf, int D, int H, double *theta, RNG *r);
CNFSampleResult cnf_sample(const CNF *cnf, const double *theta,
const double *z0, double t0, double t1,
double atol, double rtol);
CNFLogProbResult cnf_log_prob(const CNF *cnf, const double *theta,
const double *z1, double t0, double t1,
double atol, double rtol);
CNFBackwardResult cnf_backward(const CNF *cnf, const double *theta,
const double *z1,
double dL_dlogp, const double *dL_dz1_in,
double t0, double t1,
double atol, double rtol);
|