From 1e1fdfc2e3892667bccaf317a01defd8832041c7 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Mon, 16 Jan 2017 01:31:10 -0500
Subject: fixed voltage and torus conditions, current and free boundaries and
 broken right now

---
 src/factor_update.c | 30 ++++++++++++++++++++++++++++++
 1 file changed, 30 insertions(+)

(limited to 'src/factor_update.c')

diff --git a/src/factor_update.c b/src/factor_update.c
index f27971b..a51705a 100644
--- a/src/factor_update.c
+++ b/src/factor_update.c
@@ -36,3 +36,33 @@ void factor_update(cholmod_factor *factor, uint_t v1, uint_t v2, cholmod_common
 	CHOL_F(free_sparse)(&perm_update_mat, c);
 	CHOL_F(free_sparse)(&update_mat, c);
 }
+
+void factor_update2(cholmod_factor *factor, uint_t v, cholmod_common *c) {
+	uint_t n = factor->n;
+
+	cholmod_sparse *update_mat =
+			CHOL_F(allocate_sparse)(n, n, 1, true, true, 0, CHOLMOD_REAL, c);
+
+	int_t *pp = (int_t *)update_mat->p;
+	int_t *ii = (int_t *)update_mat->i;
+	double *xx = (double *)update_mat->x;
+
+	for (uint_t i = 0; i <= v; i++) {
+		pp[i] = 0;
+	}
+
+	for (uint_t i = v + 1; i <= n; i++) {
+		pp[i] = 1;
+	}
+
+	ii[0] = v;
+	xx[0] = 1;
+
+	cholmod_sparse *perm_update_mat = CHOL_F(submatrix)(
+			update_mat, factor->Perm, factor->n, NULL, -1, true, true, c);
+
+	CHOL_F(updown)(false, perm_update_mat, factor, c);
+
+	CHOL_F(free_sparse)(&perm_update_mat, c);
+	CHOL_F(free_sparse)(&update_mat, c);
+}
-- 
cgit v1.2.3-70-g09d2