source: trunk/example/example.cpp @ 3746

Last change on this file since 3746 was 3746, checked in by bradbell, 4 years ago

merge to branch: trunk
from repository: https://github.com/coin-or/CppAD
start hash code: 57d3215cae5e9da7b4c92e89f038d70870ef7327
end hash code: 9aebc1ca2c0949dec7c2d156517db26e60f28159

commit 9aebc1ca2c0949dec7c2d156517db26e60f28159
Author: Brad Bell <bradbell@…>
Date: Sun Nov 8 20:15:38 2015 -0800

Remove invisible white space.

commit a92ac50e9f4c8d0007ea5a245b3e23145dfcebfe
Author: Brad Bell <bradbell@…>
Date: Sun Nov 8 20:15:31 2015 -0800

Use vectorBool with partial sparsity patterns per pass to reduce memory requirements.


solve_callback.hpp: remove invisible white space.
rev_sparse_jac.hpp: fix bug (argument transposed).
bool_sparsity.cpp: remove invisible white space.

commit c09744b13ba2c70d6ffa857206d45560154d800a
Author: Brad Bell <bradbell@…>
Date: Sun Nov 8 03:22:57 2015 -0800

check_for_nan.hpp: fix minor type in documentation.

  • Property svn:keywords set to Id
File size: 13.6 KB
Line 
1/* $Id: example.cpp 3746 2015-11-09 04:50:27Z bradbell $ */
2/* --------------------------------------------------------------------------
3CppAD: C++ Algorithmic Differentiation: Copyright (C) 2003-15 Bradley M. Bell
4
5CppAD is distributed under multiple licenses. This distribution is under
6the terms of the
7                    Eclipse Public License Version 1.0.
8
9A copy of this license is included in the COPYING file of this distribution.
10Please visit http://www.coin-or.org/CppAD/ for information on other licenses.
11-------------------------------------------------------------------------- */
12
13/*
14$begin example.cpp$$
15$spell
16        Cpp
17$$
18
19$section CppAD Examples and Tests$$
20
21$head Running Tests$$
22To build this program and run its correctness tests see $cref cmake_check$$.
23
24$code
25$verbatim%example/example.cpp%0%// BEGIN C++%// END C++%1%$$
26$$
27
28$end
29*/
30// BEGIN C++
31
32// system include files used for I/O
33# include <iostream>
34
35// C style asserts
36# include <cassert>
37
38// standard string
39# include <string>
40
41// memory utility
42# include <cppad/thread_alloc.hpp>
43
44// prototype external compiled tests (this line expected by bin/new_test.sh)
45extern bool abort_recording(void);
46extern bool ad_assign(void);
47extern bool ad_ctor(void);
48extern bool abs(void);
49extern bool acos(void);
50extern bool acosh(void);
51extern bool Add(void);
52extern bool AddEq(void);
53extern bool ad_fun(void);
54extern bool ad_in_c(void);
55extern bool ad_input(void);
56extern bool ad_output(void);
57extern bool asin(void);
58extern bool asinh(void);
59extern bool atan(void);
60extern bool atanh(void);
61extern bool atan2(void);
62extern bool azmul(void);
63extern bool base_require(void);
64extern bool BenderQuad(void);
65extern bool BoolFun(void);
66extern bool bool_sparsity(void);
67extern bool capacity_order(void);
68extern bool change_param(void);
69extern bool check_for_nan(void);
70extern bool CheckNumericType(void);
71extern bool CheckSimpleVector(void);
72extern bool colpack_hes(void);
73extern bool colpack_jac(void);
74extern bool Compare(void);
75extern bool compare_change(void);
76extern bool complex_poly(void);
77extern bool CondExp(void);
78extern bool conj_grad(void);
79extern bool Cos(void);
80extern bool Cosh(void);
81extern bool CppAD_vector(void);
82extern bool dependency(void);
83extern bool Div(void);
84extern bool DivEq(void);
85extern bool eigen_array(void);
86extern bool eigen_det(void);
87extern bool EqualOpSeq(void);
88extern bool Erf(void);
89extern bool ErrorHandler(void);
90extern bool exp(void);
91extern bool expm1(void);
92extern bool ForOne(void);
93extern bool ForTwo(void);
94extern bool ForSparseJac(void);
95extern bool Forward(void);
96extern bool forward_dir(void);
97extern bool forward_order(void);
98extern bool fun_assign(void);
99extern bool FunCheck(void);
100extern bool HesLagrangian(void);
101extern bool HesLuDet(void);
102extern bool HesMinorDet(void);
103extern bool Hessian(void);
104extern bool HesTimesDir(void);
105extern bool Independent(void);
106extern bool index_sort(void);
107extern bool Integer(void);
108extern bool Interface2C(void);
109extern bool interp_onetape(void);
110extern bool interp_retape(void);
111extern bool JacLuDet(void);
112extern bool JacMinorDet(void);
113extern bool Jacobian(void);
114extern bool log(void);
115extern bool log10(void);
116extern bool log1p(void);
117extern bool LuFactor(void);
118extern bool LuInvert(void);
119extern bool LuRatio(void);
120extern bool LuSolve(void);
121extern bool LuVecADOk(void);
122extern bool Mul(void);
123extern bool MulEq(void);
124extern bool mul_level(void);
125extern bool mul_level_adolc(void);
126extern bool mul_level_ode(void);
127extern bool mul_level_adolc_ode(void);
128extern bool nan(void);
129extern bool Near_Equal(void);
130extern bool NearEqualExt(void);
131extern bool number_skip(void);
132extern bool NumericType(void);
133extern bool num_limits(void);
134extern bool OdeErrControl(void);
135extern bool OdeErrMaxabs(void);
136extern bool OdeGear(void);
137extern bool OdeGearControl(void);
138extern bool OdeStiff(void);
139extern bool ode_taylor(void);
140extern bool opt_val_hes(void);
141extern bool optimize(void);
142extern bool ParVar(void);
143extern bool Poly(void);
144extern bool pow(void);
145extern bool pow_int(void);
146extern bool print_for(void);
147extern bool reverse_any(void);
148extern bool reverse_one(void);
149extern bool reverse_three(void);
150extern bool reverse_two(void);
151extern bool RevOne(void);
152extern bool RevSparseHes(void);
153extern bool RevSparseJac(void);
154extern bool RevTwo(void);
155extern bool RombergMul(void);
156extern bool RombergOne(void);
157extern bool Rosen34(void);
158extern bool runge_45_1(void);
159extern bool runge_45_2(void);
160extern bool seq_property(void);
161extern bool sign(void);
162extern bool SimpleVector(void);
163extern bool Sin(void);
164extern bool Sinh(void);
165extern bool sparse_hessian(void);
166extern bool sparse_jacobian(void);
167extern bool sparse_sub_hes(void);
168extern bool sparsity_sub(void);
169extern bool Sqrt(void);
170extern bool StackMachine(void);
171extern bool Sub(void);
172extern bool SubEq(void);
173extern bool sub_sparse_hes(void);
174extern bool Tan(void);
175extern bool Tanh(void);
176extern bool TapeIndex(void);
177extern bool thread_alloc(void);
178extern bool UnaryMinus(void);
179extern bool UnaryPlus(void);
180extern bool Value(void);
181extern bool Var2Par(void);
182extern bool vec_ad(void);
183extern bool vectorBool(void);
184
185namespace {
186        // function that runs one test
187        static size_t Run_ok_count    = 0;
188        static size_t Run_error_count = 0;
189        bool Run(bool TestOk(void), std::string name)
190        {       bool ok      = true;
191                std::streamsize width =  20;
192                std::cout.width( width );
193                std::cout.setf( std::ios_base::left );
194                std::cout << name;
195                //
196                ok &= name.size() < size_t(width);
197                ok &= TestOk();
198                if( ok )
199                {       std::cout << "OK" << std::endl;
200                        Run_ok_count++;
201                }
202                else
203                {       std::cout << "Error" << std::endl;
204                        Run_error_count++;
205                }
206                return ok;
207        }
208}
209
210// main program that runs all the tests
211int main(void)
212{       bool ok = true;
213
214        // This line is used by test_one.sh
215
216        // run external compiled tests (this line expected by bin/new_test.sh)
217        ok &= Run( abort_recording,   "abort_recording"  );
218        ok &= Run( ad_assign,         "ad_assign"        );
219        ok &= Run( ad_ctor,           "ad_ctor"          );
220        ok &= Run( abs,               "abs"              );
221        ok &= Run( acos,              "acos"             );
222        ok &= Run( acosh,             "acosh"            );
223        ok &= Run( Add,               "Add"              );
224        ok &= Run( AddEq,             "AddEq"            );
225        ok &= Run( ad_fun,            "ad_fun"           );
226        ok &= Run( ad_in_c,           "ad_in_c"          );
227        ok &= Run( ad_input,          "ad_input"         );
228        ok &= Run( ad_output,         "ad_output"        );
229        ok &= Run( asin,              "asin"             );
230        ok &= Run( asinh,             "asinh"            );
231        ok &= Run( atan,              "atan"             );
232        ok &= Run( atanh,             "atanh"            );
233        ok &= Run( atan2,             "atan2"            );
234        ok &= Run( azmul,             "azmul"            );
235        ok &= Run( BenderQuad,        "BenderQuad"       );
236        ok &= Run( BoolFun,           "BoolFun"          );
237        ok &= Run( bool_sparsity,     "bool_sparsity"    );
238        ok &= Run( capacity_order,    "capacity_order"   );
239        ok &= Run( change_param,      "change_param"     );
240        ok &= Run( check_for_nan,     "check_for_nan"    );
241        ok &= Run( CheckNumericType,  "CheckNumericType" );
242        ok &= Run( CheckSimpleVector, "CheckSimpleVector");
243        ok &= Run( Compare,           "Compare"          );
244        ok &= Run( compare_change,    "compare_change"   );
245        ok &= Run( complex_poly,      "complex_poly"     );
246        ok &= Run( CondExp,           "CondExp"          );
247        ok &= Run( conj_grad,         "conj_grad"        );
248        ok &= Run( Cos,               "Cos"              );
249        ok &= Run( Cosh,              "Cosh"             );
250        ok &= Run( CppAD_vector,      "CppAD_vector"     );
251        ok &= Run( dependency,        "dependency"       );
252        ok &= Run( Div,               "Div"              );
253        ok &= Run( DivEq,             "DivEq"            );
254        ok &= Run( EqualOpSeq,        "EqualOpSeq"       );
255        ok &= Run( Erf,               "Erf"              );
256        ok &= Run( ErrorHandler,      "ErrorHandler"     );
257        ok &= Run( exp,               "exp"              );
258        ok &= Run( expm1,             "expm1"            );
259        ok &= Run( ForOne,            "ForOne"           );
260        ok &= Run( ForTwo,            "ForTwo"           );
261        ok &= Run( Forward,           "Forward"          );
262        ok &= Run( forward_dir,       "forward_dir"      );
263        ok &= Run( forward_order,     "forward_order"    );
264        ok &= Run( ForSparseJac,      "ForSparseJac"     );
265        ok &= Run( fun_assign,        "fun_assign"       );
266        ok &= Run( FunCheck,          "FunCheck"         );
267        ok &= Run( HesLagrangian,     "HesLagrangian"    );
268        ok &= Run( HesLuDet,          "HesLuDet"         );
269        ok &= Run( HesMinorDet,       "HesMinorDet"      );
270        ok &= Run( Hessian,           "Hessian"          );
271        ok &= Run( HesTimesDir,       "HesTimesDir"      );
272        ok &= Run( Independent,       "Independent"      );
273        ok &= Run( index_sort,        "index_sort"       );
274        ok &= Run( Integer,           "Integer"          );
275        ok &= Run( Interface2C,       "Interface2C"      );
276        ok &= Run( interp_onetape,    "interp_onetape"   );
277        ok &= Run( interp_retape,     "interp_retape"    );
278        ok &= Run( JacLuDet,          "JacLuDet"         );
279        ok &= Run( JacMinorDet,       "JacMinorDet"      );
280        ok &= Run( Jacobian,          "Jacobian"         );
281        ok &= Run( log,               "log"              );
282        ok &= Run( log10,             "log10"            );
283        ok &= Run( log1p,             "log1p"            );
284        ok &= Run( LuFactor,          "LuFactor"         );
285        ok &= Run( LuInvert,          "LuInvert"         );
286        ok &= Run( LuRatio,           "LuRatio"          );
287        ok &= Run( LuSolve,           "LuSolve"          );
288        ok &= Run( LuVecADOk,         "LuVecADOk"        );
289        ok &= Run( Mul,               "Mul"              );
290        ok &= Run( MulEq,             "MulEq"            );
291        ok &= Run( mul_level,         "mul_level"        );
292        ok &= Run( mul_level_ode,     "mul_level_ode"    );
293        ok &= Run( nan,               "nan"              );
294        ok &= Run( Near_Equal,        "Near_Equal"       );
295        ok &= Run( NearEqualExt,      "NearEqualExt"     );
296        ok &= Run( number_skip,       "number_skip"      );
297        ok &= Run( NumericType,       "NumericType"      );
298        ok &= Run( num_limits,        "num_limits"       );
299        ok &= Run( OdeErrControl,     "OdeErrControl"    );
300        ok &= Run( OdeErrMaxabs,      "OdeErrMaxabs"     );
301        ok &= Run( OdeGear,           "OdeGear"          );
302        ok &= Run( OdeGearControl,    "OdeGearControl"   );
303        ok &= Run( OdeStiff,          "OdeStiff"         );
304        ok &= Run( ode_taylor,        "ode_taylor"       );
305        ok &= Run( opt_val_hes,       "opt_val_hes"      );
306        ok &= Run( optimize,          "optimize"         );
307        ok &= Run( ParVar,            "ParVar"           );
308        ok &= Run( Poly,              "Poly"             );
309        ok &= Run( pow,               "pow"              );
310        ok &= Run( pow_int,           "pow_int"          );
311        ok &= Run( reverse_any,       "reverse_any"      );
312        ok &= Run( reverse_one,       "reverse_one"      );
313        ok &= Run( reverse_three,     "reverse_three"    );
314        ok &= Run( reverse_two,       "reverse_two"      );
315        ok &= Run( RevOne,            "RevOne"           );
316        ok &= Run( RevSparseHes,      "RevSparseHes"     );
317        ok &= Run( RevSparseJac,      "RevSparseJac"     );
318        ok &= Run( RevTwo,            "RevTwo"           );
319        ok &= Run( RombergMul,        "RombergMul"       );
320        ok &= Run( RombergOne,        "RombergOne"       );
321        ok &= Run( Rosen34,           "Rosen34"          );
322        ok &= Run( runge_45_1,        "runge_45_1"       );
323        ok &= Run( runge_45_2,        "runge_45_2"       );
324        ok &= Run( seq_property,      "seq_property"     );
325        ok &= Run( sign,              "sign"             );
326        ok &= Run( SimpleVector,      "SimpleVector"     );
327        ok &= Run( Sin,               "Sin"              );
328        ok &= Run( Sinh,              "Sinh"             );
329        ok &= Run( sparse_hessian,    "sparse_hessian"   );
330        ok &= Run( sparse_jacobian,   "sparse_jacobian"  );
331        ok &= Run( sparse_sub_hes,    "sparse_sub_hes"   );
332        ok &= Run( sparsity_sub,      "sparsity_sub"     );
333        ok &= Run( Sqrt,              "Sqrt"             );
334        ok &= Run( StackMachine,      "StackMachine"     );
335        ok &= Run( Sub,               "Sub"              );
336        ok &= Run( SubEq,             "SubEq"            );
337        ok &= Run( sub_sparse_hes,    "sub_sparse_hes"   );
338        ok &= Run( Tan,               "Tan"              );
339        ok &= Run( Tanh,              "Tanh"             );
340        ok &= Run( TapeIndex,         "TapeIndex"        );
341        ok &= Run( thread_alloc,      "thread_alloc"     );
342        ok &= Run( UnaryMinus,        "UnaryMinus"       );
343        ok &= Run( UnaryPlus,         "UnaryPlus"        );
344        ok &= Run( Value,             "Value"            );
345        ok &= Run( Var2Par,           "Var2Par"          );
346        ok &= Run( vec_ad,            "vec_ad"           );
347        ok &= Run( vectorBool,        "vectorBool"       );
348# ifdef CPPAD_ADOLC_EXAMPLES
349        ok &= Run( mul_level_adolc,      "mul_level_adolc"     );
350        ok &= Run( mul_level_adolc_ode,  "mul_level_adolc_ode" );
351# endif
352# ifdef CPPAD_COLPACK_EXAMPLES
353        ok &= Run( colpack_hes,       "colpack_hes"      );
354        ok &= Run( colpack_jac,       "colpack_jac"      );
355# endif
356# ifdef CPPAD_EIGEN_EXAMPLES
357        ok &= Run( eigen_array,       "eigen_array"      );
358        ok &= Run( eigen_det,         "eigen_det"        );
359# endif
360
361        // check for errors
362        using std::cout;
363        using std::endl;
364        assert( ok || (Run_error_count > 0) );
365        if( CppAD::thread_alloc::free_all() )
366        {       Run_ok_count++;
367                cout << "OK:    " << "No memory leak detected" << endl;
368        }
369        else
370        {       ok = false;
371                Run_error_count++;
372                cout << "Error: " << "memory leak detected" << endl;
373        }
374        // Run base_require after memory leak check because base_alloc.hpp uses
375        // thread_alloc to allocate memory for static copies of nan.
376        ok &= Run( base_require,      "base_require"     );
377        // convert int(size_t) to avoid warning on _MSC_VER systems
378        if( ok )
379                cout << "All " << int(Run_ok_count) << " tests passed." << endl;
380        else    cout << int(Run_error_count) << " tests failed." << endl;
381
382        return static_cast<int>( ! ok );
383}
384// END C++
Note: See TracBrowser for help on using the repository browser.