clean up

Author: dance858

src/bivariate/left_matmul.c

@@ -200,7 +200,10 @@ expr *new_left_matmul(expr *param_node, expr *child, const CSR_Matrix *A)
     node->left = child;
     expr_retain(child);
 
-    /* Store small A (NOT block-diagonal) — block functions handle the rest */
+    /* Store small A (NOT block-diagonal) — block functions handle the rest
+       Allocate workspace. iwork is used for converting J_child csr to csc
+       (requring size node->n_vars). csc_to_csr_workspace is used for
+       converting J_CSC to CSR (requring node->size) */
     node->iwork = (int *) malloc(node->n_vars * sizeof(int));
     lin_node->AT_iwork = (int *) malloc(A->n * sizeof(int));
     lin_node->csc_to_csr_workspace = (int *) malloc(node->size * sizeof(int));

src/bivariate/vector_mult.c

@@ -27,9 +27,6 @@
 static void forward(expr *node, const double *u)
 {
     expr *child = node->left;
-    // vector_mult_expr *vn = (vector_mult_expr *) node;
-    // const double *a = vn->param_source->value;
-
     const double *a = ((vector_mult_expr *) node)->param_source->value;
 
     /* child's forward pass */
@@ -58,9 +55,6 @@ static void jacobian_init(expr *node)
 static void eval_jacobian(expr *node)
 {
     expr *x = node->left;
-    // vector_mult_expr *vn = (vector_mult_expr *) node;
-    // const double *a = vn->param_source->value;
-
     const double *a = ((vector_mult_expr *) node)->param_source->value;
 
     /* evaluate x */
@@ -94,9 +88,6 @@ static void wsum_hess_init(expr *node)
 static void eval_wsum_hess(expr *node, const double *w)
 {
     expr *x = node->left;
-    // vector_mult_expr *vn = (vector_mult_expr *) node;
-    // const double *a = vn->param_source->value;
-
     const double *a = ((vector_mult_expr *) node)->param_source->value;
 
     /* scale weights w by a */