diff --git a/NewKernel_d/include/CGAL/NewKernel_d/Types/Weighted_point.h b/NewKernel_d/include/CGAL/NewKernel_d/Types/Weighted_point.h
index 8f0a80280b1..815766787b6 100644
--- a/NewKernel_d/include/CGAL/NewKernel_d/Types/Weighted_point.h
+++ b/NewKernel_d/include/CGAL/NewKernel_d/Types/Weighted_point.h
@@ -186,11 +186,10 @@ template <class R_> struct Power_center : Store_kernel<R_> {
     for(i=0; ++f!=e; ++i) {
       WPoint const& wp=*f;
       Point const& p=pdw(wp);
-      FT const& np = sdo(p) - pw(wp);
       for(int j=0;j<d;++j) {
 	m(i,j)=2*(c(p,j)-c(p0,j));
-	b[i] = np - n0;
       }
+      b[i] = sdo(p) - pw(wp) - n0;
     }
     CGAL_assertion (i == d);
     Vec res = typename CVec::Dimension()(d);;
diff --git a/NewKernel_d/include/CGAL/NewKernel_d/function_objects_cartesian.h b/NewKernel_d/include/CGAL/NewKernel_d/function_objects_cartesian.h
index cf82d299a45..66e8f75c275 100644
--- a/NewKernel_d/include/CGAL/NewKernel_d/function_objects_cartesian.h
+++ b/NewKernel_d/include/CGAL/NewKernel_d/function_objects_cartesian.h
@@ -640,8 +640,8 @@ template <class R_> struct Construct_circumcenter : Store_kernel<R_> {
 	Point const& p=*f;
 	for(int j=0;j<d;++j) {
 	  m(i,j)=2*(c(p,j)-c(p0,j));
-	  b[i] = sdo(p) - n0;
 	}
+	b[i] = sdo(p) - n0;
       }
       CGAL_assertion (i == d);
       Vec res = typename CVec::Dimension()(d);;