Codebase list cppad / upstream/2015.00.00.7 doc / exp_eps_for2.cpp.xml
upstream/2015.00.00.7

Tree @upstream/2015.00.00.7 (Download .tar.gz)

exp_eps_for2.cpp.xml @upstream/2015.00.00.7 — raw · history · blame 

  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
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
<?xml version='1.0'?>
<html xmlns='http://www.w3.org/1999/xhtml'
      xmlns:math='http://www.w3.org/1998/Math/MathML'
>
<head>
<title>exp_eps: Verify Second Order Forward Sweep</title>
<meta http-equiv='Content-Type' content='text/html' charset='utf-8'/>
<meta name="description" id="description" content="exp_eps: Verify Second Order Forward Sweep"/>
<meta name="keywords" id="keywords" content=" exp_eps: verify second order forward sweep first exp_2 "/>
<style type='text/css'>
body { color : black }
body { background-color : white }
A:link { color : blue }
A:visited { color : purple }
A:active { color : purple }
</style>
<script type='text/javascript' language='JavaScript' src='_exp_eps_for2.cpp_xml.js'>
</script>
</head>
<body>
<table><tr>
<td>
<a href="http://www.coin-or.org/CppAD/" target="_top"><img border="0" src="_image.gif"/></a>
</td>
<td><a href="exp_eps_for2.xml" target="_top">Prev</a>
</td><td><a href="exp_eps_rev2.xml" target="_top">Next</a>
</td><td>
<select onchange='choose_across0(this)'>
<option>Index-&gt;</option>
<option>contents</option>
<option>reference</option>
<option>index</option>
<option>search</option>
<option>external</option>
</select>
</td>
<td>
<select onchange='choose_up0(this)'>
<option>Up-&gt;</option>
<option>CppAD</option>
<option>Introduction</option>
<option>exp_eps</option>
<option>exp_eps_for2</option>
<option>exp_eps_for2.cpp</option>
</select>
</td>
<td>
<select onchange='choose_down3(this)'>
<option>Introduction-&gt;</option>
<option>get_started.cpp</option>
<option>exp_2</option>
<option>exp_eps</option>
<option>exp_apx_main.cpp</option>
</select>
</td>
<td>
<select onchange='choose_down2(this)'>
<option>exp_eps-&gt;</option>
<option>exp_eps.hpp</option>
<option>exp_eps.cpp</option>
<option>exp_eps_for0</option>
<option>exp_eps_for1</option>
<option>exp_eps_rev1</option>
<option>exp_eps_for2</option>
<option>exp_eps_rev2</option>
<option>exp_eps_cppad</option>
</select>
</td>
<td>
<select onchange='choose_down1(this)'>
<option>exp_eps_for2-&gt;</option>
<option>exp_eps_for2.cpp</option>
</select>
</td>
<td>exp_eps_for2.cpp</td>
<td>Headings</td>
</tr></table><br/>



<center><b><big><big>exp_eps: Verify Second Order Forward Sweep</big></big></b></center>
<code><font color='blue'><pre style='display:inline'> # include &lt;cmath&gt;                     // for fabs function
extern bool exp_eps_for0(double *v0); // computes zero order forward sweep
extern bool exp_eps_for1(double *v1); // computes first order forward sweep
bool exp_eps_for2(void)
{	bool ok = true;
	double v0[8], v1[8], v2[8];

	// set the value of v0[j], v1[j] for j = 1 , ... , 7
	ok &amp;= exp_eps_for0(v0);
	ok &amp;= exp_eps_for1(v1);

	v2[1] = 0.;                                      // v1 = x
	ok    &amp;= std::fabs( v2[1] - 0. ) &lt;= 1e-10;

	v2[2] = 1. * v2[1];                              // v2 = 1 * v1
	ok    &amp;= std::fabs( v2[2] - 0. ) &lt;= 1e-10;

	v2[3] = v2[2] / 1.;                              // v3 = v2 / 1
	ok    &amp;= std::fabs( v2[3] - 0. ) &lt;= 1e-10;

	v2[4] = v2[3];                                   // v4 = 1 + v3
	ok    &amp;= std::fabs( v2[4] - 0. ) &lt;= 1e-10;

	v2[5] = v2[3] * v0[1] + 2. * v1[3] * v1[1]       // v5 = v3 * v1
	      + v0[3] * v2[1];           
	ok    &amp;= std::fabs( v2[5] - 2. ) &lt;= 1e-10;

	v2[6] = v2[5] / 2.;                              // v6 = v5 / 2
	ok    &amp;= std::fabs( v2[6] - 1. ) &lt;= 1e-10;

	v2[7] = v2[4] + v2[6];                           // v7 = v4 + v6
	ok    &amp;= std::fabs( v2[7] - 1. ) &lt;= 1e-10;

	return ok;
}
</pre></font></code>


<hr/>Input File: introduction/exp_apx/exp_eps_for2.cpp

</body>
</html>