remove default parameter from dense matrix constructor

Author: Ahdhn

apps/NeoHookean/neo_hookean.cu

@@ -45,7 +45,7 @@ void neo_hookean(RXMeshStatic& rx, T dx)
     const T        stiffness_coef = 4e4;
     const T        tol            = 0.01;
     const T        inv_time_step  = T(1) / time_step;
-    DenseMatrix<T> dbc_stiff(1, 1);
+    DenseMatrix<T> dbc_stiff(1, 1, LOCATION_ALL);
     dbc_stiff(0) = 1000;
     dbc_stiff.move(HOST, DEVICE);
     const T dhat  = 0.01;
@@ -169,10 +169,10 @@ void neo_hookean(RXMeshStatic& rx, T dx)
         problem, contact_area, time_step, ground_n, ground_o, dhat, kappa);
 
     DenseMatrix<T, Eigen::RowMajor> dir(
-        rx, problem.grad.rows(), problem.grad.cols());
+        rx, problem.grad.rows(), problem.grad.cols(), LOCATION_ALL);
 
     DenseMatrix<T, Eigen::RowMajor> grad(
-        rx, problem.grad.rows(), problem.grad.cols());
+        rx, problem.grad.rows(), problem.grad.cols(), LOCATION_ALL);
 
     T line_search_init_step = 0;
 

include/rxmesh/diff/candidate_pairs.h

@@ -44,10 +44,13 @@ struct CandidatePairs
           m_variable_dim(hess.K_),
           m_pairs_id(DenseMatrix<IndexT, Eigen::ColMajor>(
               max_capacity * m_variable_dim * m_variable_dim * 2,
-              2)),
-          m_pairs_handle(DenseMatrix<PairT, Eigen::ColMajor>(max_capacity, 1)),
-          m_current_num_pairs(DenseMatrix<int>(1, 1)),
-          m_current_num_index(DenseMatrix<int>(1, 1)),
+              2,
+              LOCATION_ALL)),
+          m_pairs_handle(DenseMatrix<PairT, Eigen::ColMajor>(max_capacity,
+                                                             1,
+                                                             LOCATION_ALL)),
+          m_current_num_pairs(DenseMatrix<int>(1, 1, LOCATION_ALL)),
+          m_current_num_index(DenseMatrix<int>(1, 1, LOCATION_ALL)),
           m_context(ctx)
     {
         reset();

include/rxmesh/diff/diff_scalar_problem.h

@@ -54,7 +54,10 @@ struct DiffScalarProblem
                       bool          assmble_hessian,
                       int           expected_vv_candidate_pairs = 0)
         : rx(rx),
-          grad(DenseMatT(rx, rx.get_num_elements<ObjHandleT>(), VariableDim)),
+          grad(DenseMatT(rx,
+                         rx.get_num_elements<ObjHandleT>(),
+                         VariableDim,
+                         LOCATION_ALL)),
           objective(rx.add_vertex_attribute<T>("objective", VariableDim))
     {
         grad.reset(0, LOCATION_ALL);

include/rxmesh/diff/newton_solver.h

@@ -33,7 +33,7 @@ struct NetwtonSolver
      */
     NetwtonSolver(DiffProblemT& p, SolverT* s)
         : problem(p),
-          dir(DenseMatT(p.rx, p.grad.rows(), p.grad.cols())),
+          dir(DenseMatT(p.rx, p.grad.rows(), p.grad.cols(), LOCATION_ALL)),
           temp_objective(
               p.rx.add_attribute_like("temp_objective", *p.objective)),
           solver(s),

include/rxmesh/matrix/dense_matrix.h

@@ -61,10 +61,10 @@ struct DenseMatrix
      * of vertices, edges, or faces. With constructor, the user can access
      * the matrix using mesh handles
      */
-    DenseMatrix(const RXMesh& rx,
-                IndexT        num_rows,
-                IndexT        num_cols,
-                locationT     location = LOCATION_ALL)
+    explicit DenseMatrix(const RXMesh& rx,
+                         IndexT        num_rows,
+                         IndexT        num_cols,
+                         locationT     location)
         : m_context(rx.get_context()),
           m_num_rows(num_rows),
           m_num_cols(num_cols),
@@ -84,9 +84,7 @@ struct DenseMatrix
      * So the size (num_rows) can be anything. However, using this constructor
      * means that you can not access the matrix using handles
      */
-    DenseMatrix(IndexT    num_rows,
-                IndexT    num_cols,
-                locationT location = LOCATION_ALL)
+    explicit DenseMatrix(IndexT num_rows, IndexT num_cols, locationT location)
         : m_num_rows(num_rows),
           m_num_cols(num_cols),
           m_dendescr(NULL),