fdurastante/psblas3
1
0
Fork 0
mirror of https://github.com/sfilippone/psblas3.git synced 2026-10-20 12:52:14 +00:00
Code Issues Packages Projects Releases Wiki Activity

First version of complex SuperLU preconditioners.

Browse Source
This commit is contained in:
Salvatore Filippone
2006-04-21 14:05:47 +00:00
parent 66955001f2
commit 41d7626b66
3 changed files with 15 additions and 16 deletions
Download Patch File Download Diff File Expand all files Collapse all files
src/prec/psb_slu_impl.c
+6 -6
View File
@@ -131,7 +131,7 @@ typedef struct {
void
psb_dslu_factor_(int *n, int *nnz,
double *values, int *rowind, int *colptr,
double *values, int *rowptr, int *colind,
#ifdef Have_SLU_
fptr *f_factors, /* a handle containing the address
pointing to the factored matrices */
@@ -177,10 +177,10 @@ psb_dslu_factor_(int *n, int *nnz,
StatInit(&stat);
/* Adjust to 0-based indexing */
for (i = 0; i < *nnz; ++i) --colptr[i];
for (i = 0; i <= *n; ++i) --rowind[i];
for (i = 0; i < *nnz; ++i) --colind[i];
for (i = 0; i <= *n; ++i) --rowptr[i];
dCreate_CompRow_Matrix(&A, *n, *n, *nnz, values, colptr, rowind,
dCreate_CompRow_Matrix(&A, *n, *n, *nnz, values, colind, rowptr,
SLU_NR, SLU_D, SLU_GE);
L = (SuperMatrix *) SUPERLU_MALLOC( sizeof(SuperMatrix) );
U = (SuperMatrix *) SUPERLU_MALLOC( sizeof(SuperMatrix) );
@@ -230,8 +230,8 @@ psb_dslu_factor_(int *n, int *nnz,
}
/* Restore to 1-based indexing */
for (i = 0; i < *nnz; ++i) ++colptr[i];
for (i = 0; i <= *n; ++i) ++rowind[i];
for (i = 0; i < *nnz; ++i) ++colind[i];
for (i = 0; i <= *n; ++i) ++rowptr[i];
/* Save the LU factors in the factors handle */
LUfactors = (factors_t*) SUPERLU_MALLOC(sizeof(factors_t));
src/prec/psb_zslu_bld.f90
+1 -1
View File
@@ -177,7 +177,7 @@ subroutine psb_zslu_bld(a,desc_a,p,info)
endif
call psb_zslu_factor(atmp%m,nzt,&
& atmp%aspk,atmp%ia2,atmp%ia1,p%iprcparm(slu_ptr_),info)
& atmp%aspk,atmp%ia1,atmp%ia2,p%iprcparm(slu_ptr_),info)
if(info /= 0) then
info=4010
ch_err='psb_slu_fact'
src/prec/psb_zslu_impl.c
+8 -9
View File
@@ -85,7 +85,6 @@ SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
#ifdef Have_SLU_
#include "zsp_defs.h"
#undef Have_SLU_
#define HANDLE_SIZE 8
/* kind of integer to hold a pointer. Use int.
This might need to be changed on 64-bit systems. */
@@ -132,11 +131,11 @@ typedef struct {
void
psb_zslu_factor_(int *n, int *nnz,
#ifdef Have_SLU_
doublecomplex *values, int *rowind, int *colptr,
doublecomplex *values, int *colind, int *rowptr,
fptr *f_factors, /* a handle containing the address
pointing to the factored matrices */
#else
void *values, int *rowind, int *colptr,
void *values, int *colind, int *rowptr,
void *f_factors,
#endif
int *info)
@@ -178,11 +177,11 @@ psb_zslu_factor_(int *n, int *nnz,
StatInit(&stat);
/* Adjust to 0-based indexing */
for (i = 0; i < *nnz; ++i) --colptr[i];
for (i = 0; i <= *n; ++i) --rowind[i];
for (i = 0; i < *nnz; ++i) --colind[i];
for (i = 0; i <= *n; ++i) --rowptr[i];
zCreate_CompRow_Matrix(&A, *n, *n, *nnz, values, colptr, rowind,
SLU_NC, SLU_Z, SLU_GE);
zCreate_CompRow_Matrix(&A, *n, *n, *nnz, values, colind, rowptr,
SLU_NR, SLU_Z, SLU_GE);
L = (SuperMatrix *) SUPERLU_MALLOC( sizeof(SuperMatrix) );
U = (SuperMatrix *) SUPERLU_MALLOC( sizeof(SuperMatrix) );
if ( !(perm_r = intMalloc(*n)) ) ABORT("Malloc fails for perm_r[].");
@@ -231,8 +230,8 @@ psb_zslu_factor_(int *n, int *nnz,
}
/* Restore to 1-based indexing */
for (i = 0; i < *nnz; ++i) ++colptr[i];
for (i = 0; i <= *n; ++i) ++rowind[i];
for (i = 0; i < *nnz; ++i) ++colind[i];
for (i = 0; i <= *n; ++i) ++rowptr[i];
/* Save the LU factors in the factors handle */
LUfactors = (factors_t*) SUPERLU_MALLOC(sizeof(factors_t));