00001 # ifndef CPPAD_REVERSE_INCLUDED
00002 # define CPPAD_REVERSE_INCLUDED
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039 # include <algorithm>
00040
00041
00042 namespace CppAD {
00043
00044 template <typename Base>
00045 template <typename VectorBase>
00046 VectorBase ADFun<Base>::Reverse(size_t p, const VectorBase &w) const
00047 {
00048 size_t i, j, k;
00049
00050
00051 size_t n = ind_taddr.size();
00052
00053
00054 size_t m = dep_taddr.size();
00055
00056 Base *Partial = CPPAD_NULL;
00057 Partial = CPPAD_TRACK_NEW_VEC(totalNumVar * p, Partial);
00058
00059
00060
00061
00062
00063
00064
00065 CheckSimpleVector<Base, VectorBase>();
00066
00067 CPPAD_ASSERT_KNOWN(
00068 w.size() == m,
00069 "Argument w to Reverse does not have length equal to\n"
00070 "the dimension of the range for the corresponding ADFun."
00071 );
00072 CPPAD_ASSERT_KNOWN(
00073 p > 0,
00074 "The first argument to Reverse must be greater than zero."
00075 );
00076 CPPAD_ASSERT_KNOWN(
00077 taylor_per_var >= p,
00078 "Less that p Taylor coefficients are currently stored"
00079 " in this ADFun object."
00080 );
00081
00082
00083 for(i = 0; i < totalNumVar; i++)
00084 for(j = 0; j < p; j++)
00085 Partial[i * p + j] = Base(0);
00086
00087
00088
00089 for(i = 0; i < m; i++)
00090 { CPPAD_ASSERT_UNKNOWN( dep_taddr[i] < totalNumVar );
00091 Partial[dep_taddr[i] * p + p - 1] += w[i];
00092 }
00093
00094
00095 ReverseSweep(
00096 p - 1, totalNumVar, &Rec, TaylorColDim, Taylor, p, Partial
00097 );
00098
00099
00100 VectorBase value(n * p);
00101 for(j = 0; j < n; j++)
00102 { CPPAD_ASSERT_UNKNOWN( ind_taddr[j] < totalNumVar );
00103
00104
00105 CPPAD_ASSERT_UNKNOWN( Rec.GetOp( ind_taddr[j] ) == InvOp );
00106
00107
00108
00109
00110 for(k = 0; k < p; k++)
00111 value[j * p + k ] =
00112 Partial[ind_taddr[j] * p + p - 1 - k];
00113 }
00114
00115
00116 CPPAD_TRACK_DEL_VEC(Partial);
00117
00118 return value;
00119 }
00120
00121
00122 }
00123
00124
00125 # endif