Tested the c++ codes run on GPUs and fixed unit tests associated with them

What is the bug

The c++ codes were not running in parallel using Kokkos and were not tested properly

How do you fix it

How urgent is the bugfix

• I need it as soon as possible
• I can wait for a couple of days
• None of my current codes is directly affected

Mandatory steps before review

• Gitlab CI passes (Hint: use make format for linting)
• Bugfix is covered by additional unit tests
• Mark the merge request as ready by removing Draft:

Mandatory steps before merge

• Reviewed by a maintainer
• Incorporate review suggestions
• Remember to edit the commit message and select the proper changelog category (feature/bugfix/other)
• Prior to merging, please remove any boilerplate from the MR description, retaining only the What is the bug and How do you fix it section to maintain

You are not supposed to merge this request by yourself, the maintainers of libiconmath take care of this action!

test/c/test_horizontal_div.cpp

@@ -117,7 +117,6 @@ TYPED_TEST(HorizontalDivTest, TestDiv3DSpecific) {
       cell_edge_blk_h[cell_edge_at(i, 0, j)] = 0;
     }
 
-    // Geometric factors
     geofac_div_h[geofac_div_at(i, 0, 0)] = static_cast<TypeParam>(0.5);
     geofac_div_h[geofac_div_at(i, 1, 0)] = static_cast<TypeParam>(0.3);
     geofac_div_h[geofac_div_at(i, 2, 0)] = static_cast<TypeParam>(0.2);
@@ -153,8 +152,6 @@ TYPED_TEST(HorizontalDivTest, TestDiv3DSpecific) {
   EXPECT_NEAR(div_vec_c_h(div_vec_c_at(2, 1, 0)), static_cast<TypeParam>(4.4), 1e-6);
 }
 
-// write another unit test for the same function where the input arrays are filled with random values
-
 TYPED_TEST(HorizontalDivTest, TestDiv3DRandom) {
   constexpr int nproma = this->nproma;
   constexpr int nlev = this->nlev;
@@ -479,8 +476,6 @@ TYPED_TEST(HorizontalDivTest, TestDiv3D2FRandom) {
   }
 }
 
-// Third unit test: div4d
-
 TYPED_TEST(HorizontalDivTest, TestDiv4DSpecific) {
   constexpr int nproma = this->nproma;
   constexpr int nlev = this->nlev;