master: valvector: combine split and join in split_join.hpp and finis…

…h out AD version of reverse.

CMakeLists.txt

@@ -16,7 +16,7 @@ IF( POLICY CMP0054 )
 ENDIF( POLICY CMP0054 )
 #
 # cppad_version is used by version.sh to get the version number.
-SET(cppad_version "20240301")
+SET(cppad_version "20240302")
 SET(cppad_url          "https://coin-or.github.io/CppAD" )
 SET(cppad_description  "Differentiation of C++ Algorithms" )
 IF( NOT DEFINED CMAKE_BUILD_TYPE)

appendix/whats_new/2024.xrst

@@ -28,6 +28,11 @@ Release Notes for 2024
 mm-dd
 *****
 
+03-02
+=====
+Add the :ref:`valvector_ad_split-name` and :ref:`valvector_ad_join-name`
+examples.
+
 03-01
 =====
 Fix bug in the optimizer :ref:`optimize@options@val_graph` option.

example/valvector/ad_join.cpp

@@ -2,7 +2,7 @@
 // SPDX-FileCopyrightText: Bradley M. Bell <[email protected]>
 // SPDX-FileContributor: 2024 Bradley M. Bell
 // ---------------------------------------------------------------------------
-# include <cppad/example/valvector/ad_join.hpp>
+# include <cppad/example/valvector/split_join.hpp>
 # include <cppad/cppad.hpp>
 /*
 {xrst_begin valvector_ad_join.cpp}
@@ -153,6 +153,19 @@ bool ad_join(void)
    for(size_t j = 0; j < n; ++j)
       ok &= dz[0][j] == 2.0 * x[j][0];
    //
+   // h
+   CPPAD_TESTVECTOR( ad_valvector ) aw(m), adw(n);
+   aw[0] = valvector( 1.0 );
+   CppAD::Independent(ax);
+   af.Forward(0, ax);
+   adw = af.Reverse(1, aw);
+   CppAD::ADFun<valvector> h(ax, adw);
+   //
+   // ok
+   dw = h.Forward(0, x);
+   for(size_t j = 0; j < n; ++j)
+      ok &= dw[j][0] == 2.0 * x[j][0];
+   //
    // ok
    f.optimize();
    z = f.Forward(0, x);

example/valvector/ad_split.cpp

@@ -2,7 +2,7 @@
 // SPDX-FileCopyrightText: Bradley M. Bell <[email protected]>
 // SPDX-FileContributor: 2024 Bradley M. Bell
 // ---------------------------------------------------------------------------
-# include <cppad/example/valvector/ad_split.hpp>
+# include <cppad/example/valvector/split_join.hpp>
 # include <cppad/cppad.hpp>
 /*
 {xrst_begin valvector_ad_split.cpp}
@@ -135,6 +135,20 @@ bool ad_split(void)
    for(size_t i = 0; i < m; ++i)
       ok &= dy[i][0] == 2.0 * x[0][i];
    //
+   // h
+   CPPAD_TESTVECTOR( ad_valvector ) aw(m), adw(n);
+   for(size_t i = 0; i < m; ++i)
+      aw[i] = valvector( 1.0 );
+   CppAD::Independent(ax);
+   af.Forward(0, ax);
+   adw = af.Reverse(1, aw);
+   CppAD::ADFun<valvector> h(ax, adw);
+   //
+   // ok
+   dw = h.Forward(0, x);
+   for(size_t i = 0; i < m; ++i)
+      ok &= dw[0][i] == 2.0 * x[0][i];
+   //
    // ok
    f.optimize();
    y = f.Forward(0, x);