1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
From: Ghislain Antony Vaillant <ghisvail@gmail.com>
Date: Fri, 23 Sep 2016 12:26:08 +0100
Subject: superlu 5.x
---
Src/superlu3module.c | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/Src/superlu3module.c b/Src/superlu3module.c
index 810191b..cb96ae9 100644
--- a/Src/superlu3module.c
+++ b/Src/superlu3module.c
@@ -193,6 +193,7 @@ static PyObject *newSuperLUObject(int n, CSRMatObject *matrix, PyObject *issym,
SuperLUObject *self;
SuperMatrix A; /* A in NC format used by the factorization routine. */
SuperMatrix AC; /* Matrix postmultiplied by Pc */
+ GlobalLU_t Glu;
mem_usage_t mem_usage;
int *etree;
int info = 0;
@@ -267,9 +268,9 @@ static PyObject *newSuperLUObject(int n, CSRMatObject *matrix, PyObject *issym,
//printf("Factorizing with drop_tol = %g\n", drop_tol);
/* Perform factorization (perm_c and perm_r are swapped because our matrix
* is stored in compressed-row format and not in compressed-column format! */
- dgstrf(&self->options, &AC, drop_tol, relax, panel_size,
+ dgstrf(&self->options, &AC, relax, panel_size,
etree, NULL, 0, self->perm_c, self->perm_r,
- &self->L, &self->U, &self->stat, &info);
+ &self->L, &self->U, &Glu, &self->stat, &info);
//dPrint_SuperNode_Matrix("L", &self->L);
//dPrint_CompCol_Matrix("U", &self->U);
|
