fixed bug in test and syntax for thrust::transform

Author: cnpetra

src/LinAlg/VectorCudaKernels.cu

@@ -1062,7 +1062,7 @@ double log_barr_wei_obj_kernel(int n, double* d1, const double* id, const double
   // compute v[i] = ( log(v[i]) if v[i]>0, otherwise 0 )
   thrust::transform(thrust::device, v_temp, v_temp+n, v_temp, log_select_op); 
   // compute v[i] = w[i]*v[i] 
-  thrust::transform(thrust::device, v_temp, v_temp+n, v_temp, mult_op); 
+  thrust::transform(thrust::device, v_temp, v_temp+n, wei_v, v_temp, mult_op); 
   // sum up
   const double sum = thrust::reduce(thrust::device, v_temp, v_temp+n, 0.0, plus_op);
 

tests/LinAlg/vectorTests.hpp

@@ -1240,8 +1240,8 @@ class VectorTests : public TestBase
     setLocalElement(&weight, N-1, 100);
 
     real_type expected = (N - 1) * std::log(x_val) / N;
-    real_type result = x.logBarrier_local(pattern);
-
+    real_type result = x.logBarrierWeighted_local(pattern, weight);
+    
     int fail = !isEqual(result, expected);
 
     // Make sure pattern eliminates the correct elements
@@ -1253,7 +1253,7 @@ class VectorTests : public TestBase
     setLocalElement(&weight, N - 1, w_val);
 
     expected = std::log(x_val)*w_val;
-    result = x.logBarrier_local(pattern);
+    result = x.logBarrierWeighted_local(pattern, weight);
     fail += !isEqual(result, expected);
 
     printMessage(fail, __func__, rank);