Skip to content

Commit

Permalink
Add support for Kokkos::View element access in the rvs mode (#1068)
Browse files Browse the repository at this point in the history
  • Loading branch information
gojakuch authored Oct 9, 2024
1 parent 4f2a378 commit a14a3f6
Show file tree
Hide file tree
Showing 2 changed files with 220 additions and 5 deletions.
216 changes: 216 additions & 0 deletions include/clad/Differentiator/KokkosBuiltins.h
Original file line number Diff line number Diff line change
Expand Up @@ -30,6 +30,37 @@ constructor_pushforward(
Kokkos::View<DataType, ViewParams...>(
"_diff_" + name, idx0, idx1, idx2, idx3, idx4, idx5, idx6, idx7)};
}
template <class DataType, class... ViewParams>
clad::ValueAndAdjoint<::Kokkos::View<DataType, ViewParams...>,
::Kokkos::View<DataType, ViewParams...>>
constructor_reverse_forw(
clad::ConstructorReverseForwTag<::Kokkos::View<DataType, ViewParams...>>,
const ::std::string& name, const size_t& idx0, const size_t& idx1,
const size_t& idx2, const size_t& idx3, const size_t& idx4,
const size_t& idx5, const size_t& idx6, const size_t& idx7,
const ::std::string& /*d_name*/, const size_t& /*d_idx0*/,
const size_t& /*d_idx1*/, const size_t& /*d_idx2*/,
const size_t& /*d_idx3*/, const size_t& /*d_idx4*/,
const size_t& /*d_idx5*/, const size_t& /*d_idx6*/,
const size_t& /*d_idx7*/) {
return {::Kokkos::View<DataType, ViewParams...>(name, idx0, idx1, idx2, idx3,
idx4, idx5, idx6, idx7),
::Kokkos::View<DataType, ViewParams...>(
"_diff_" + name, idx0, idx1, idx2, idx3, idx4, idx5, idx6, idx7)};
}
template <class DataType, class... ViewParams>
void constructor_pullback(::Kokkos::View<DataType, ViewParams...>* v,
const ::std::string& name, const size_t& idx0,
const size_t& idx1, const size_t& idx2,
const size_t& idx3, const size_t& idx4,
const size_t& idx5, const size_t& idx6,
const size_t& idx7,
::Kokkos::View<DataType, ViewParams...>* d_v,
const ::std::string* /*d_name*/,
const size_t& /*d_idx0*/, const size_t* /*d_idx1*/,
const size_t* /*d_idx2*/, const size_t* /*d_idx3*/,
const size_t* /*d_idx4*/, const size_t* /*d_idx5*/,
const size_t* /*d_idx6*/, const size_t* /*d_idx7*/) {}

/// View indexing
template <typename View, typename Idx>
Expand Down Expand Up @@ -107,6 +138,191 @@ operator_call_pushforward(const View* v, Idx0 i0, Idx1 i1, Idx2 i2, Idx3 i3,
return {(*v)(i0, i1, i2, i3, i4, i5, i6, i7),
(*d_v)(i0, i1, i2, i3, i4, i5, i6, i7)};
}
template <class DataType, class... ViewParams, typename Idx>
clad::ValueAndAdjoint<
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&,
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&>
operator_call_reverse_forw(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx i0,
const ::Kokkos::View<DataType, ViewParams...>* d_v,
Idx /*d_i0*/) {
return {(*v)(i0), (*d_v)(i0)};
}
template <class DataType, class... ViewParams, typename Diff, typename Idx,
typename dIdx>
void operator_call_pullback(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx i0, Diff d_y,
::Kokkos::View<DataType, ViewParams...>* d_v,
dIdx* /*d_i0*/) {
(*d_v)(i0) += d_y;
}
template <class DataType, class... ViewParams, typename Idx0, typename Idx1>
clad::ValueAndAdjoint<
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&,
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&>
operator_call_reverse_forw(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx0 i0, Idx1 i1,
const ::Kokkos::View<DataType, ViewParams...>* d_v,
Idx0 /*d_i0*/, Idx1 /*d_i1*/) {
return {(*v)(i0, i1), (*d_v)(i0, i1)};
}
template <class DataType, class... ViewParams, typename Diff, typename Idx0,
typename dIdx0, typename Idx1, typename dIdx1>
void operator_call_pullback(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx0 i0, Idx1 i1, Diff d_y,
::Kokkos::View<DataType, ViewParams...>* d_v,
dIdx0* /*d_i0*/, dIdx1* /*d_i1*/) {
(*d_v)(i0, i1) += d_y;
}
template <class DataType, class... ViewParams, typename Idx0, typename Idx1,
typename Idx2>
clad::ValueAndAdjoint<
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&,
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&>
operator_call_reverse_forw(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx0 i0, Idx1 i1, Idx2 i2,
const ::Kokkos::View<DataType, ViewParams...>* d_v,
Idx0 /*d_i0*/, Idx1 /*d_i1*/, Idx2 /*d_i2*/) {
return {(*v)(i0, i1, i2), (*d_v)(i0, i1, i2)};
}
template <class DataType, class... ViewParams, typename Diff, typename Idx0,
typename dIdx0, typename Idx1, typename dIdx1, typename Idx2,
typename dIdx2>
void operator_call_pullback(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx0 i0, Idx1 i1, Idx2 i2, Diff d_y,
::Kokkos::View<DataType, ViewParams...>* d_v,
dIdx0* /*d_i0*/, dIdx1* /*d_i1*/, dIdx2* /*d_i2*/) {
(*d_v)(i0, i1, i2) += d_y;
}
template <class DataType, class... ViewParams, typename Idx0, typename Idx1,
typename Idx2, typename Idx3>
clad::ValueAndAdjoint<
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&,
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&>
operator_call_reverse_forw(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx0 i0, Idx1 i1, Idx2 i2, Idx3 i3,
const ::Kokkos::View<DataType, ViewParams...>* d_v,
Idx0 /*d_i0*/, Idx1 /*d_i1*/, Idx2 /*d_i2*/,
Idx3 /*d_i3*/) {
return {(*v)(i0, i1, i2, i3), (*d_v)(i0, i1, i2, i3)};
}
template <class DataType, class... ViewParams, typename Diff, typename Idx0,
typename dIdx0, typename Idx1, typename dIdx1, typename Idx2,
typename dIdx2, typename Idx3, typename dIdx3>
void operator_call_pullback(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx0 i0, Idx1 i1, Idx2 i2, Idx3 i3, Diff d_y,
::Kokkos::View<DataType, ViewParams...>* d_v,
dIdx0* /*d_i0*/, dIdx1* /*d_i1*/, dIdx2* /*d_i2*/,
dIdx3* /*d_i3*/) {
(*d_v)(i0, i1, i2, i3) += d_y;
}
template <class DataType, class... ViewParams, typename Idx0, typename Idx1,
typename Idx2, typename Idx3, typename Idx4>
clad::ValueAndAdjoint<
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&,
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&>
operator_call_reverse_forw(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx0 i0, Idx1 i1, Idx2 i2, Idx3 i3, Idx4 i4,
const ::Kokkos::View<DataType, ViewParams...>* d_v,
Idx0 /*d_i0*/, Idx1 /*d_i1*/, Idx2 /*d_i2*/,
Idx3 /*d_i3*/, Idx4 /*d_i4*/) {
return {(*v)(i0, i1, i2, i3, i4), (*d_v)(i0, i1, i2, i3, i4)};
}
template <class DataType, class... ViewParams, typename Diff, typename Idx0,
typename dIdx0, typename Idx1, typename dIdx1, typename Idx2,
typename dIdx2, typename Idx3, typename dIdx3, typename Idx4,
typename dIdx4>
void operator_call_pullback(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx0 i0, Idx1 i1, Idx2 i2, Idx3 i3, Idx4 i4,
Diff d_y,
::Kokkos::View<DataType, ViewParams...>* d_v,
dIdx0* /*d_i0*/, dIdx1* /*d_i1*/, dIdx2* /*d_i2*/,
dIdx3* /*d_i3*/, dIdx4* /*d_i4*/) {
(*d_v)(i0, i1, i2, i3, i4) += d_y;
}
template <class DataType, class... ViewParams, typename Idx0, typename Idx1,
typename Idx2, typename Idx3, typename Idx4, typename Idx5>
clad::ValueAndAdjoint<
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&,
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&>
operator_call_reverse_forw(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx0 i0, Idx1 i1, Idx2 i2, Idx3 i3, Idx4 i4, Idx5 i5,
const ::Kokkos::View<DataType, ViewParams...>* d_v,
Idx0 /*d_i0*/, Idx1 /*d_i1*/, Idx2 /*d_i2*/,
Idx3 /*d_i3*/, Idx4 /*d_i4*/, Idx5 /*d_i5*/) {
return {(*v)(i0, i1, i2, i3, i4, i5), (*d_v)(i0, i1, i2, i3, i4, i5)};
}
template <class DataType, class... ViewParams, typename Diff, typename Idx0,
typename dIdx0, typename Idx1, typename dIdx1, typename Idx2,
typename dIdx2, typename Idx3, typename dIdx3, typename Idx4,
typename dIdx4, typename Idx5, typename dIdx5>
void operator_call_pullback(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx0 i0, Idx1 i1, Idx2 i2, Idx3 i3, Idx4 i4,
Idx5 i5, Diff d_y,
::Kokkos::View<DataType, ViewParams...>* d_v,
dIdx0* /*d_i0*/, dIdx1* /*d_i1*/, dIdx2* /*d_i2*/,
dIdx3* /*d_i3*/, dIdx4* /*d_i4*/, dIdx5* /*d_i5*/) {
(*d_v)(i0, i1, i2, i3, i4, i5) += d_y;
}
template <class DataType, class... ViewParams, typename Idx0, typename Idx1,
typename Idx2, typename Idx3, typename Idx4, typename Idx5,
typename Idx6>
clad::ValueAndAdjoint<
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&,
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&>
operator_call_reverse_forw(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx0 i0, Idx1 i1, Idx2 i2, Idx3 i3, Idx4 i4, Idx5 i5,
Idx6 i6,
const ::Kokkos::View<DataType, ViewParams...>* d_v,
Idx0 /*d_i0*/, Idx1 /*d_i1*/, Idx2 /*d_i2*/,
Idx3 /*d_i3*/, Idx4 /*d_i4*/, Idx5 /*d_i5*/,
Idx6 /*d_i6*/) {
return {(*v)(i0, i1, i2, i3, i4, i5, i6), (*d_v)(i0, i1, i2, i3, i4, i5, i6)};
}
template <class DataType, class... ViewParams, typename Diff, typename Idx0,
typename dIdx0, typename Idx1, typename dIdx1, typename Idx2,
typename dIdx2, typename Idx3, typename dIdx3, typename Idx4,
typename dIdx4, typename Idx5, typename dIdx5, typename Idx6,
typename dIdx6>
void operator_call_pullback(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx0 i0, Idx1 i1, Idx2 i2, Idx3 i3, Idx4 i4,
Idx5 i5, Idx6 i6, Diff d_y,
::Kokkos::View<DataType, ViewParams...>* d_v,
dIdx0* /*d_i0*/, dIdx1* /*d_i1*/, dIdx2* /*d_i2*/,
dIdx3* /*d_i3*/, dIdx4* /*d_i3*/, dIdx5* /*d_i3*/,
dIdx6* /*d_i3*/) {
(*d_v)(i0, i1, i2, i3, i4, i5, i6) += d_y;
}
template <class DataType, class... ViewParams, typename Idx0, typename Idx1,
typename Idx2, typename Idx3, typename Idx4, typename Idx5,
typename Idx6, typename Idx7>
clad::ValueAndAdjoint<
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&,
typename ::Kokkos::View<DataType, ViewParams...>::reference_type&>
operator_call_reverse_forw(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx0 i0, Idx1 i1, Idx2 i2, Idx3 i3, Idx4 i4, Idx5 i5,
Idx6 i6, Idx7 i7,
const ::Kokkos::View<DataType, ViewParams...>* d_v,
Idx0 /*d_i0*/, Idx1 /*d_i1*/, Idx2 /*d_i2*/,
Idx3 /*d_i3*/, Idx4 /*d_i4*/, Idx5 /*d_i5*/,
Idx6 /*d_i6*/, Idx7 /*d_i7*/) {
return {(*v)(i0, i1, i2, i3, i4, i5, i6, i7),
(*d_v)(i0, i1, i2, i3, i4, i5, i6, i7)};
}
template <class DataType, class... ViewParams, typename Diff, typename Idx0,
typename dIdx0, typename Idx1, typename dIdx1, typename Idx2,
typename dIdx2, typename Idx3, typename dIdx3, typename Idx4,
typename dIdx4, typename Idx5, typename dIdx5, typename Idx6,
typename dIdx6, typename Idx7, typename dIdx7>
void operator_call_pullback(const ::Kokkos::View<DataType, ViewParams...>* v,
Idx0 i0, Idx1 i1, Idx2 i2, Idx3 i3, Idx4 i4,
Idx5 i5, Idx6 i6, Idx7 i7, Diff d_y,
::Kokkos::View<DataType, ViewParams...>* d_v,
dIdx0* /*d_i0*/, dIdx1* /*d_i1*/, dIdx2* /*d_i2*/,
dIdx3* /*d_i3*/, dIdx4* /*d_i3*/, dIdx5* /*d_i3*/,
dIdx6* /*d_i3*/, dIdx7* /*d_i3*/) {
(*d_v)(i0, i1, i2, i3, i4, i5, i6, i7) += d_y;
}
} // namespace class_functions

/// Kokkos functions (view utils)
Expand Down
9 changes: 4 additions & 5 deletions unittests/Kokkos/ViewAccess.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -60,11 +60,10 @@ TEST(ViewAccess, Test2) {
double dx_f_2_FD = finite_difference_tangent(f_2_tmp, 3., epsilon);
EXPECT_NEAR(f_2_x.execute(3, 4), dx_f_2_FD, tolerance * dx_f_2_FD);

// TODO: uncomment this once it has been implemented
// auto f_grad_exe = clad::gradient(f);
// double dx, dy;
// f_grad_exe.execute(3., 4., &dx, &dy);
// EXPECT_NEAR(f_x.execute(3, 4),dx,tolerance*dx);
auto f_grad_exe = clad::gradient(f);
double dx, dy;
f_grad_exe.execute(3., 4., &dx, &dy);
EXPECT_NEAR(f_x.execute(3, 4), dx, tolerance * dx);

// double dx_2, dy_2;
// auto f_2_grad_exe = clad::gradient(f_2);
Expand Down

0 comments on commit a14a3f6

Please sign in to comment.