test2.cpp File Reference

Detailed Description

Basic test of FlowBalloonModel.

This test configures a FlowBalloonModel with parameters and functions set according to the second experiment in Buxton et al. (1998). Results are output on stdout, tab delimited, with columns representing:

; time
$f_{in}(t)$
$f_{out}(t)$
$CMRO_2(t) = \frac{E(t)}{E_0}f_{in}(t)$ (Buxton et al. 2004).
; BOLD response

Results are as follows:

Results, c.f. Buxton et al. (1998), Figure 2

Note:: The only difference between this test and that in test1.cpp is the definition of F_OUT -- linear in the case of test1.cpp and nonlinear in the case of test1.cpp.
Look out for a double line in the graph for $f_{out}(v)$ which may indicate that the error bounds need to be tightened or that discontinuities in $f_{in}(t)$ are not being handled well. As the volume increases, peaks and then decreases with time in this test, this graph has two lines drawn. They should be exactly the same.

Definition in file test2.cpp.

Go to the source code of this file.

Functions

double F_IN (double t, const double y[], void *o)

$f_{in}(t)$ ; Trapezoidal function with a rise time of 4 s, duration of 4 s and fall time of 4 s.

double F_OUT (double t, const double y[], void *o)

$f_{out}(t) = \frac{7}{80}\left[\frac{20}{3}(v - 1)\right]^{3} + 1$ ;

int main (int argc, const char *argv[])

Run tests.

Variables

static const double TAU_0 = 2.0

$\tau_0$

static const double E_0 = 0.4

$E_0$

static const double V_0 = 0.01

$V_0$

static const double ALPHA = 0.5

$\alpha$

static const double END = 20.0

Duration of simulation (s).

Function Documentation

double F_IN	(	double	t,
		const double	y[],
		void *	o
	)

$f_{in}(t)$ ; Trapezoidal function with a rise time of 4 s, duration of 4 s and fall time of 4 s.

Definition at line 77 of file test2.cpp.

double F_OUT	(	double	t,
		const double	y[],
		void *	o
	)

$f_{out}(t) = \frac{7}{80}\left[\frac{20}{3}(v - 1)\right]^{3} + 1$ ;

i.e. nonlinear relative to $v$ . The exact form of the function used in Buxton et al. (1998) Figure 2 is unknown, but this seems a reasonable approximation.

Definition at line 99 of file test2.cpp.

int main	(	int	argc,
		const char *	argv[]
	)

Run tests.

Definition at line 106 of file test2.cpp.

Variable Documentation

const double ALPHA = 0.5 [static]

$\alpha$

Definition at line 66 of file test2.cpp.

const double E_0 = 0.4 [static]

$E_0$

Definition at line 56 of file test2.cpp.

const double END = 20.0 [static]

Duration of simulation (s).

Definition at line 71 of file test2.cpp.

const double TAU_0 = 2.0 [static]

$\tau_0$

Definition at line 51 of file test2.cpp.

const double V_0 = 0.01 [static]

$V_0$

Definition at line 61 of file test2.cpp.

Generated on Mon Aug 13 19:51:39 2007 for fmrii Hemodynamic Models Test Suite by

1.5.2


Functions
double	F_IN (double t, const double y[], void *o)
	$f_{in}(t)$ ; Trapezoidal function with a rise time of 4 s, duration of 4 s and fall time of 4 s.
double	F_OUT (double t, const double y[], void *o)
	$f_{out}(t) = \frac{7}{80}\left[\frac{20}{3}(v - 1)\right]^{3} + 1$ ;
int	main (int argc, const char *argv[])
	Run tests.
Variables
static const double	TAU_0 = 2.0
	$\tau_0$
static const double	E_0 = 0.4

static const double	V_0 = 0.01

static const double	ALPHA = 0.5
	$\alpha$
static const double	END = 20.0
	Duration of simulation (s).