move old code into old-code-folder (#68)

* move old code into old-code-folder

* new folder structure

Author: dance858

include/affine.h include/atoms/affine.hinclude/affine.h renamed to include/atoms/affine.h

@@ -22,8 +22,6 @@
 #include "subexpr.h"
 #include "utils/CSR_Matrix.h"
 
-expr *new_linear(expr *u, const CSR_Matrix *A, const double *b);
-
 expr *new_add(expr *left, expr *right);
 expr *new_neg(expr *child);
 

include/bivariate_full_dom.h include/atoms/bivariate_full_dom.hinclude/bivariate_full_dom.h renamed to include/atoms/bivariate_full_dom.h

@@ -15,8 +15,8 @@
  * See the License for the specific language governing permissions and
  * limitations under the License.
  */
-#ifndef OTHER_H
-#define OTHER_H
+#ifndef NON_ELEMENTWISE_FULL_DOM_H
+#define NON_ELEMENTWISE_FULL_DOM_H
 
 #include "expr.h"
 #include "subexpr.h"
@@ -33,4 +33,4 @@ expr *new_prod_axis_zero(expr *child);
 /* product of entries along axis=1 (rowwise products) */
 expr *new_prod_axis_one(expr *child);
 
-#endif /* OTHER_H */
+#endif /* NON_ELEMENTWISE_FULL_DOM_H */

include/bivariate_restricted_dom.h include/atoms/bivariate_restricted_dom.hinclude/bivariate_restricted_dom.h renamed to include/atoms/bivariate_restricted_dom.h

@@ -0,0 +1,29 @@
+#ifndef OLD_CSR_H
+#define OLD_CSR_H
+
+#include "utils/CSR_Matrix.h"
+
+/* Build (I_p kron A) = blkdiag(A, A, ..., A) of size (p*A->m) x (p*A->n) */
+CSR_Matrix *block_diag_repeat_csr(const CSR_Matrix *A, int p);
+
+/* Build (A kron I_p) of size (A->m * p) x (A->n * p) with nnz = A->nnz * p. */
+CSR_Matrix *kron_identity_csr(const CSR_Matrix *A, int p);
+
+/* Computes values of the row matrix C = z^T A (column indices must have been
+   pre-computed) and transposed matrix AT must be provided) */
+void Ax_csr_fill_values(const CSR_Matrix *AT, const double *z, CSR_Matrix *C);
+
+/* Insert value into CSR matrix A with just one row at col_idx. Assumes that A
+has enough space and that A does not have an element at col_idx. It does update
+nnz. */
+void csr_insert_value(CSR_Matrix *A, int col_idx, double value);
+
+/* Compute C = diag(d) * A where d is an array and A, C are CSR matrices
+ * d must have length m
+ * C must be pre-allocated with same dimensions as A */
+void diag_csr_mult(const double *d, const CSR_Matrix *A, CSR_Matrix *C);
+
+/* y = Ax, where y is returned as dense (no column offset) */
+void Ax_csr_wo_offset(const CSR_Matrix *A, const double *x, double *y);
+
+#endif /* OLD_CSR_H */

include/elementwise_full_dom.h include/atoms/elementwise_full_dom.hinclude/elementwise_full_dom.h renamed to include/atoms/elementwise_full_dom.h

@@ -0,0 +1,44 @@
+#ifndef OLD_CSR_SUM_H
+#define OLD_CSR_SUM_H
+
+#include "utils/CSR_Matrix.h"
+
+/* Compute C = A + B where A, B, C are CSR matrices
+ * A and B must have same dimensions
+ * C must be pre-allocated with sufficient nnz capacity.
+ * C must be different from A and B */
+void sum_csr_matrices(const CSR_Matrix *A, const CSR_Matrix *B, CSR_Matrix *C);
+
+/* Compute C = diag(d1) * A + diag(d2) * B where A, B, C are CSR matrices */
+void sum_scaled_csr_matrices(const CSR_Matrix *A, const CSR_Matrix *B, CSR_Matrix *C,
+                             const double *d1, const double *d2);
+
+/* forward declaration */
+struct int_double_pair;
+
+/* Sum all rows of A into a single row matrix C */
+void sum_all_rows_csr(const CSR_Matrix *A, CSR_Matrix *C,
+                      struct int_double_pair *pairs);
+
+/* Sum blocks of rows of A into a matrix C */
+void sum_block_of_rows_csr(const CSR_Matrix *A, CSR_Matrix *C,
+                           struct int_double_pair *pairs, int row_block_size);
+
+/* Sum evenly spaced rows of A into a matrix C */
+void sum_evenly_spaced_rows_csr(const CSR_Matrix *A, CSR_Matrix *C,
+                                struct int_double_pair *pairs, int row_spacing);
+
+/* Sum evenly spaced rows of A starting at offset into a row matrix C */
+void sum_spaced_rows_into_row_csr(const CSR_Matrix *A, CSR_Matrix *C,
+                                  struct int_double_pair *pairs, int offset,
+                                  int spacing);
+
+/* Fill values of summed rows using precomputed idx_map and sparsity of C */
+void sum_all_rows_csr_fill_values(const CSR_Matrix *A, CSR_Matrix *C,
+                                  const int *idx_map);
+
+/* Fill values of summed block rows using precomputed idx_map */
+void sum_block_of_rows_csr_fill_values(const CSR_Matrix *A, CSR_Matrix *C,
+                                       const int *idx_map);
+
+#endif /* OLD_CSR_SUM_H */

include/elementwise_restricted_dom.h …clude/atoms/elementwise_restricted_dom.hinclude/elementwise_restricted_dom.h renamed to include/atoms/elementwise_restricted_dom.h include/elementwise_restricted_dom.h renamed to include/atoms/elementwise_restricted_dom.h

@@ -0,0 +1,9 @@
+#ifndef OLD_AFFINE_H
+#define OLD_AFFINE_H
+
+#include "expr.h"
+#include "utils/CSR_Matrix.h"
+
+expr *new_linear(expr *u, const CSR_Matrix *A, const double *b);
+
+#endif /* OLD_AFFINE_H */