void MUMPS_CALL
mumps_c(MUMPS_STRUC_C * mumps_par)
{
/*
* The following local variables will
* be passed to the F77 interface.
*/
MUMPS_INT *icntl;
MUMPS_REAL *cntl;
MUMPS_INT *irn; MUMPS_INT *jcn; MUMPS_COMPLEX *a;
MUMPS_INT *irn_loc; MUMPS_INT *jcn_loc; MUMPS_COMPLEX *a_loc;
MUMPS_INT *eltptr, *eltvar; MUMPS_COMPLEX *a_elt;
MUMPS_INT *perm_in; MUMPS_INT perm_in_avail;
MUMPS_INT *listvar_schur; MUMPS_INT listvar_schur_avail;
MUMPS_COMPLEX *schur; MUMPS_INT schur_avail;
MUMPS_COMPLEX *rhs; MUMPS_COMPLEX *redrhs;
MUMPS_COMPLEX *wk_user; MUMPS_INT wk_user_avail;
MUMPS_REAL *colsca; MUMPS_REAL *rowsca;
MUMPS_COMPLEX *rhs_sparse, *sol_loc;
MUMPS_INT *irhs_sparse, *irhs_ptr, *isol_loc;
MUMPS_INT irn_avail, jcn_avail, a_avail, rhs_avail, redrhs_avail;
/* These are actually used
* as booleans, but we stick
* to simple types for the
* C-F77 interface */
MUMPS_INT irn_loc_avail, jcn_loc_avail, a_loc_avail;
MUMPS_INT eltptr_avail, eltvar_avail, a_elt_avail;
MUMPS_INT colsca_avail, rowsca_avail;
MUMPS_INT irhs_ptr_avail, rhs_sparse_avail, sol_loc_avail;
MUMPS_INT irhs_sparse_avail, isol_loc_avail;
MUMPS_INT *info; MUMPS_INT *infog;
MUMPS_REAL *rinfo; MUMPS_REAL *rinfog;
MUMPS_INT ooc_tmpdir[255]; MUMPS_INT ooc_prefix[63];
MUMPS_INT write_problem[255];
/* Other local variables */
MUMPS_INT idummy; MUMPS_INT *idummyp;
MUMPS_REAL rdummy; MUMPS_REAL *rdummyp;
MUMPS_COMPLEX cdummy; MUMPS_COMPLEX *cdummyp;
/* String lengths to be passed to Fortran by address */
int ooc_tmpdirlen;
int ooc_prefixlen;
int write_problemlen;
int i;
const static MUMPS_INT no = 0;
const static MUMPS_INT yes = 1;
idummyp = &idummy;
cdummyp = &cdummy;
rdummyp = &rdummy;
#ifdef return_scaling
/* Don't forget to initialize those two before
* each call to mumps as we may copy values from
* old instances otherwise ! */
MUMPS_COLSCA_STATIC=0;
MUMPS_ROWSCA_STATIC=0;
#endif
/* Initialize pointers to zero for job == -1 */
if ( mumps_par->job == -1 )
{ /* job = -1: we just reset all pointers to 0 */
mumps_par->irn=0; mumps_par->jcn=0; mumps_par->a=0; mumps_par->rhs=0; mumps_par->wk_user=0;
mumps_par->redrhs=0;
mumps_par->eltptr=0; mumps_par->eltvar=0; mumps_par->a_elt=0; mumps_par->perm_in=0; mumps_par->sym_perm=0; mumps_par->uns_perm=0; mumps_par->irn_loc=0;mumps_par->jcn_loc=0;mumps_par->a_loc=0; mumps_par->listvar_schur=0;mumps_par->schur=0;mumps_par->mapping=0;mumps_par->pivnul_list=0;mumps_par->colsca=0;mumps_par->rowsca=0; mumps_par->rhs_sparse=0; mumps_par->irhs_sparse=0; mumps_par->sol_loc=0; mumps_par->irhs_ptr=0; mumps_par->isol_loc=0;
strcpy(mumps_par->ooc_tmpdir,"NAME_NOT_INITIALIZED");
strcpy(mumps_par->ooc_prefix,"NAME_NOT_INITIALIZED");
strcpy(mumps_par->write_problem,"NAME_NOT_INITIALIZED");
strncpy(mumps_par->version_number,MUMPS_VERSION,MUMPS_VERSION_MAX_LEN);
mumps_par->version_number[MUMPS_VERSION_MAX_LEN+1] = '\0';
/* Next line initializes scalars to arbitrary values.
* Some of those will anyway be overwritten during the
* call to Fortran routine [SDCZ]MUMPS_INIT_PHASE */
mumps_par->n=0; mumps_par->nz=0; mumps_par->nz_loc=0; mumps_par->nelt=0;mumps_par->instance_number=0;mumps_par->deficiency=0;mumps_par->lwk_user=0;mumps_par->size_schur=0;mumps_par->lrhs=0; mumps_par->lredrhs=0; mumps_par->nrhs=0; mumps_par->nz_rhs=0; mumps_par->lsol_loc=0;
mumps_par->schur_mloc=0; mumps_par->schur_nloc=0; mumps_par->schur_lld=0; mumps_par->mblock=0; mumps_par->nblock=0; mumps_par->nprow=0; mumps_par->npcol=0;
}
ooc_tmpdirlen=(int)strlen(mumps_par->ooc_tmpdir);
ooc_prefixlen=(int)strlen(mumps_par->ooc_prefix);
write_problemlen=(int)strlen(mumps_par->write_problem);
/* Avoid the use of strnlen which may not be
* available on all systems. Allow strings without
* \0 at the end, if the file is not found, the
* Fortran layer is responsible for raising an
* error. */
if(ooc_tmpdirlen > 255){
ooc_tmpdirlen=255;
}
if(ooc_prefixlen > 63){
ooc_prefixlen=63;
}
if(write_problemlen > 255){
write_problemlen=255;
}
/*
* Extract info from the C structure to call the F77 interface. The
* following macro avoids repeating the same code with risks of errors.
*/
#define EXTRACT_POINTERS(component,dummypointer) \
if ( mumps_par-> component == 0) \
{ component = dummypointer; \
component ## _avail = no; } \
else \
{ component = mumps_par-> component; \
component ## _avail = yes; }
/*
* For example, EXTRACT_POINTERS(irn,idummyp) produces the following line of code:
if (mumps_par->irn== 0) {irn= idummyp;irn_avail = no; } else { irn = mumps_par->irn;irn_avail = yes; } ;
* which says that irn is set to mumps_par->irn except if
* mumps_par->irn is 0, which means that it is not available.
*/
EXTRACT_POINTERS(irn,idummyp);
EXTRACT_POINTERS(jcn,idummyp);
EXTRACT_POINTERS(rhs,cdummyp);
EXTRACT_POINTERS(wk_user,cdummyp);
EXTRACT_POINTERS(redrhs,cdummyp);
EXTRACT_POINTERS(irn_loc,idummyp);
EXTRACT_POINTERS(jcn_loc,idummyp);
EXTRACT_POINTERS(a_loc,cdummyp);
EXTRACT_POINTERS(a,cdummyp);
EXTRACT_POINTERS(eltptr,idummyp);
EXTRACT_POINTERS(eltvar,idummyp);
EXTRACT_POINTERS(a_elt,cdummyp);
EXTRACT_POINTERS(perm_in,idummyp);
EXTRACT_POINTERS(listvar_schur,idummyp);
EXTRACT_POINTERS(schur,cdummyp);
EXTRACT_POINTERS(colsca,rdummyp);
EXTRACT_POINTERS(rowsca,rdummyp);
EXTRACT_POINTERS(rhs_sparse,cdummyp);
EXTRACT_POINTERS(sol_loc,cdummyp);
EXTRACT_POINTERS(irhs_sparse,idummyp);
EXTRACT_POINTERS(isol_loc,idummyp);
EXTRACT_POINTERS(irhs_ptr,idummyp);
/* printf("irn_avail,jcn_avail, rhs_avail, a_avail, eltptr_avail, eltvar_avail,a_elt_avail,perm_in_avail= %d %d %d %d %d %d %d \n", irn_avail,jcn_avail, rhs_avail, a_avail, eltptr_avail, eltvar_avail, a_elt_avail, perm_in_avail);*/
/*
* Extract integers (input) or pointers that are
* always allocated (such as ICNTL, INFO, ...)
*/
/* size_schur = mumps_par->size_schur; */
/* instance_number = mumps_par->instance_number; */
icntl = mumps_par->icntl;
cntl = mumps_par->cntl;
info = mumps_par->info;
infog = mumps_par->infog;
rinfo = mumps_par->rinfo;
rinfog = mumps_par->rinfog;
for(i=0;i<ooc_tmpdirlen;i++){
ooc_tmpdir[i]=(int)mumps_par->ooc_tmpdir[i];
}
for(i=0;i<ooc_prefixlen;i++){
ooc_prefix[i]=(int)mumps_par->ooc_prefix[i];
}
for(i=0;i<write_problemlen;i++){
write_problem[i]=(int)mumps_par->write_problem[i];
}
/* Call F77 interface */
MUMPS_F77(&(mumps_par->job), &(mumps_par->sym), &(mumps_par->par), &(mumps_par->comm_fortran),
&(mumps_par->n), icntl, cntl,
&(mumps_par->nz), irn, &irn_avail, jcn, &jcn_avail, a, &a_avail,
&(mumps_par->nz_loc), irn_loc, &irn_loc_avail, jcn_loc, &jcn_loc_avail,
a_loc, &a_loc_avail,
&(mumps_par->nelt), eltptr, &eltptr_avail, eltvar, &eltvar_avail, a_elt, &a_elt_avail,
perm_in, &perm_in_avail,
rhs, &rhs_avail, redrhs, &redrhs_avail, info, rinfo, infog, rinfog,
&(mumps_par->deficiency), &(mumps_par->lwk_user), &(mumps_par->size_schur), listvar_schur, &listvar_schur_avail, schur,
&schur_avail, wk_user, &wk_user_avail, colsca, &colsca_avail, rowsca, &rowsca_avail,
&(mumps_par->instance_number), &(mumps_par->nrhs), &(mumps_par->lrhs),
&(mumps_par->lredrhs),
rhs_sparse, &rhs_sparse_avail, sol_loc, &sol_loc_avail, irhs_sparse,
&irhs_sparse_avail, irhs_ptr, &irhs_ptr_avail, isol_loc,
&isol_loc_avail, &(mumps_par->nz_rhs), &(mumps_par->lsol_loc)
, &(mumps_par->schur_mloc)
, &(mumps_par->schur_nloc)
, &(mumps_par->schur_lld)
, &(mumps_par->mblock)
, &(mumps_par->nblock)
, &(mumps_par->nprow)
, &(mumps_par->npcol)
, ooc_tmpdir
, ooc_prefix
, write_problem
, &ooc_tmpdirlen
, &ooc_prefixlen
, &write_problemlen
);
/*
* mapping and pivnul_list are usually 0 except if
* MUMPS_AFFECT_MAPPING/MUMPS_AFFECT_PIVNUL_LIST was called.
*/
mumps_par->mapping=mumps_get_mapping();
mumps_par->pivnul_list=mumps_get_pivnul_list();
/* to get permutations computed during analysis */
mumps_par->sym_perm=mumps_get_sym_perm();
mumps_par->uns_perm=mumps_get_uns_perm();
#ifdef return_scaling
/*
* colsca/rowsca can either be user data or have been
* modified within mumps by calls to MUMPS_AFFECT_COLSCA/ROWSCA.
*/
if (colsca_avail == no) mumps_par->colsca = MUMPS_COLSCA_STATIC;
if (rowsca_avail == no) mumps_par->rowsca = MUMPS_ROWSCA_STATIC;
#endif
}
void MUMPS_CALL
mumps_c(MUMPS_STRUC_C * mumps_par)
{
/*
* The following local variables will
* be passed to the F77 interface.
*/
MUMPS_INT *icntl;
MUMPS_REAL *cntl;
MUMPS_INT *keep;
MUMPS_REAL *dkeep;
MUMPS_INT8 *keep8;
MUMPS_INT *irn; MUMPS_INT *jcn; MUMPS_COMPLEX *a;
MUMPS_INT *irn_loc; MUMPS_INT *jcn_loc; MUMPS_COMPLEX *a_loc;
MUMPS_INT *eltptr, *eltvar; MUMPS_COMPLEX *a_elt;
MUMPS_INT *perm_in; MUMPS_INT perm_in_avail;
MUMPS_INT *listvar_schur; MUMPS_INT listvar_schur_avail;
MUMPS_COMPLEX *schur; MUMPS_INT schur_avail;
MUMPS_COMPLEX *rhs; MUMPS_COMPLEX *redrhs;
MUMPS_COMPLEX *wk_user; MUMPS_INT wk_user_avail;
MUMPS_REAL *colsca; MUMPS_REAL *rowsca;
MUMPS_COMPLEX *rhs_sparse, *sol_loc;
MUMPS_INT *irhs_sparse, *irhs_ptr, *isol_loc;
MUMPS_INT irn_avail, jcn_avail, a_avail, rhs_avail, redrhs_avail;
/* These are actually used
* as booleans, but we stick
* to simple types for the
* C-F77 interface */
MUMPS_INT irn_loc_avail, jcn_loc_avail, a_loc_avail;
MUMPS_INT eltptr_avail, eltvar_avail, a_elt_avail;
MUMPS_INT colsca_avail, rowsca_avail;
MUMPS_INT irhs_ptr_avail, rhs_sparse_avail, sol_loc_avail;
MUMPS_INT irhs_sparse_avail, isol_loc_avail;
MUMPS_INT *info; MUMPS_INT *infog;
MUMPS_REAL *rinfo; MUMPS_REAL *rinfog;
MUMPS_INT ooc_tmpdir[255]; MUMPS_INT ooc_prefix[63];
MUMPS_INT write_problem[255];
/* Other local variables */
MUMPS_INT idummy; MUMPS_INT *idummyp;
MUMPS_REAL rdummy; MUMPS_REAL *rdummyp;
MUMPS_COMPLEX cdummy; MUMPS_COMPLEX *cdummyp;
/* String lengths to be passed to Fortran by address */
MUMPS_INT ooc_tmpdirlen;
MUMPS_INT ooc_prefixlen;
MUMPS_INT write_problemlen;
int i;
static const MUMPS_INT no = 0;
static const MUMPS_INT yes = 1;
idummyp = &idummy;
cdummyp = &cdummy;
rdummyp = &rdummy;
/* [SDCZ]MUMPS_F77 always calls either
* MUMPS_NULLIFY_C_COLSCA or MUMPS_ASSIGN_C_COLSCA
* (and ROWSCA). The next two lines are thus not
* strictly necessary. */
MUMPS_COLSCA_STATIC=0;
MUMPS_ROWSCA_STATIC=0;
/* Initialize pointers to zero for job == -1 */
if ( mumps_par->job == -1 )
{ /* job = -1: we just reset all pointers to 0 */
mumps_par->irn=0; mumps_par->jcn=0; mumps_par->a=0; mumps_par->rhs=0; mumps_par->wk_user=0;
mumps_par->redrhs=0;
mumps_par->eltptr=0; mumps_par->eltvar=0; mumps_par->a_elt=0; mumps_par->perm_in=0; mumps_par->sym_perm=0; mumps_par->uns_perm=0; mumps_par->irn_loc=0;mumps_par->jcn_loc=0;mumps_par->a_loc=0; mumps_par->listvar_schur=0;mumps_par->schur=0;mumps_par->mapping=0;mumps_par->pivnul_list=0;mumps_par->colsca=0;mumps_par->colsca_from_mumps=0;mumps_par->rowsca=0;mumps_par->colsca_from_mumps=0; mumps_par->rhs_sparse=0; mumps_par->irhs_sparse=0; mumps_par->sol_loc=0; mumps_par->irhs_ptr=0; mumps_par->isol_loc=0;
strcpy(mumps_par->ooc_tmpdir,"NAME_NOT_INITIALIZED");
strcpy(mumps_par->ooc_prefix,"NAME_NOT_INITIALIZED");
strcpy(mumps_par->write_problem,"NAME_NOT_INITIALIZED");
strncpy(mumps_par->version_number,MUMPS_VERSION,MUMPS_VERSION_MAX_LEN);
mumps_par->version_number[MUMPS_VERSION_MAX_LEN+1] = '\0';
/* Next line initializes scalars to arbitrary values.
* Some of those will anyway be overwritten during the
* call to Fortran routine [SDCZ]MUMPS_INIT_PHASE */
mumps_par->n=0; mumps_par->nz=0; mumps_par->nz_loc=0; mumps_par->nelt=0;mumps_par->instance_number=0;mumps_par->deficiency=0;mumps_par->lwk_user=0;mumps_par->size_schur=0;mumps_par->lrhs=0; mumps_par->lredrhs=0; mumps_par->nrhs=0; mumps_par->nz_rhs=0; mumps_par->lsol_loc=0;
mumps_par->schur_mloc=0; mumps_par->schur_nloc=0; mumps_par->schur_lld=0; mumps_par->mblock=0; mumps_par->nblock=0; mumps_par->nprow=0; mumps_par->npcol=0;
}
ooc_tmpdirlen=(int)strlen(mumps_par->ooc_tmpdir);
ooc_prefixlen=(int)strlen(mumps_par->ooc_prefix);
write_problemlen=(int)strlen(mumps_par->write_problem);
/* Avoid the use of strnlen which may not be
* available on all systems. Allow strings without
* \0 at the end, if the file is not found, the
* Fortran layer is responsible for raising an
* error. */
if(ooc_tmpdirlen > 255){
ooc_tmpdirlen=255;
}
if(ooc_prefixlen > 63){
ooc_prefixlen=63;
}
if(write_problemlen > 255){
write_problemlen=255;
}
/*
* Extract info from the C structure to call the F77 interface. The
* following macro avoids repeating the same code with risks of errors.
*/
#define EXTRACT_POINTERS(component,dummypointer) \
if ( mumps_par-> component == 0) \
{ component = dummypointer; \
component ## _avail = no; } \
else \
{ component = mumps_par-> component; \
component ## _avail = yes; }
/*
* For example, EXTRACT_POINTERS(irn,idummyp) produces the following line of code:
if (mumps_par->irn== 0) {irn= idummyp;irn_avail = no; } else { irn = mumps_par->irn;irn_avail = yes; } ;
* which says that irn is set to mumps_par->irn except if
* mumps_par->irn is 0, which means that it is not available.
*/
EXTRACT_POINTERS(irn,idummyp);
EXTRACT_POINTERS(jcn,idummyp);
EXTRACT_POINTERS(rhs,cdummyp);
EXTRACT_POINTERS(wk_user,cdummyp);
EXTRACT_POINTERS(redrhs,cdummyp);
EXTRACT_POINTERS(irn_loc,idummyp);
EXTRACT_POINTERS(jcn_loc,idummyp);
EXTRACT_POINTERS(a_loc,cdummyp);
EXTRACT_POINTERS(a,cdummyp);
EXTRACT_POINTERS(eltptr,idummyp);
EXTRACT_POINTERS(eltvar,idummyp);
EXTRACT_POINTERS(a_elt,cdummyp);
EXTRACT_POINTERS(perm_in,idummyp);
EXTRACT_POINTERS(listvar_schur,idummyp);
EXTRACT_POINTERS(schur,cdummyp);
/* EXTRACT_POINTERS not adapted to rowsca and colsca */
if ( mumps_par->rowsca != 0 && mumps_par->rowsca_from_mumps == 0 )
{
/* has been set by user and was not allocated in mumps */
rowsca = mumps_par-> rowsca;
rowsca_avail = yes;
}
else
{
/* FIXME: changing rowsca in C after an earlier call
where rowsca was computed by mumps is not possible. */
rowsca = rdummyp;
rowsca_avail = no;
}
if ( mumps_par->colsca != 0 && mumps_par->colsca_from_mumps == 0 )
/* has been changed by user and was not allocated in mumps */
{
colsca = mumps_par-> colsca;
colsca_avail = yes;
}
else
{
/* FIXME: changing colsca in C after an earlier call
where colsca was computed by mumps is not possible */
colsca = rdummyp;
colsca_avail = no;
}
EXTRACT_POINTERS(rhs_sparse,cdummyp);
EXTRACT_POINTERS(sol_loc,cdummyp);
EXTRACT_POINTERS(irhs_sparse,idummyp);
EXTRACT_POINTERS(isol_loc,idummyp);
EXTRACT_POINTERS(irhs_ptr,idummyp);
/* printf("irn_avail,jcn_avail, rhs_avail, a_avail, eltptr_avail, eltvar_avail,a_elt_avail,perm_in_avail= %d %d %d %d %d %d %d \n", irn_avail,jcn_avail, rhs_avail, a_avail, eltptr_avail, eltvar_avail, a_elt_avail, perm_in_avail); */
/*
* Extract integers (input) or pointers that are
* always allocated (such as ICNTL, INFO, ...)
*/
/* size_schur = mumps_par->size_schur; */
/* instance_number = mumps_par->instance_number; */
icntl = mumps_par->icntl;
cntl = mumps_par->cntl;
keep = mumps_par->keep;
dkeep = mumps_par->dkeep;
keep8 = mumps_par->keep8;
info = mumps_par->info;
infog = mumps_par->infog;
rinfo = mumps_par->rinfo;
rinfog = mumps_par->rinfog;
for(i=0;i<ooc_tmpdirlen;i++){
ooc_tmpdir[i]=(int)mumps_par->ooc_tmpdir[i];
}
for(i=0;i<ooc_prefixlen;i++){
ooc_prefix[i]=(int)mumps_par->ooc_prefix[i];
}
for(i=0;i<write_problemlen;i++){
write_problem[i]=(int)mumps_par->write_problem[i];
}
/* Call F77 interface */
MUMPS_F77(&(mumps_par->job), &(mumps_par->sym), &(mumps_par->par), &(mumps_par->comm_fortran),
&(mumps_par->n), icntl, cntl, keep, dkeep, keep8,
&(mumps_par->nz), irn, &irn_avail, jcn, &jcn_avail, a, &a_avail,
&(mumps_par->nz_loc), irn_loc, &irn_loc_avail, jcn_loc, &jcn_loc_avail,
a_loc, &a_loc_avail,
&(mumps_par->nelt), eltptr, &eltptr_avail, eltvar, &eltvar_avail, a_elt, &a_elt_avail,
perm_in, &perm_in_avail,
rhs, &rhs_avail, redrhs, &redrhs_avail, info, rinfo, infog, rinfog,
&(mumps_par->deficiency), &(mumps_par->lwk_user), &(mumps_par->size_schur), listvar_schur, &listvar_schur_avail, schur,
&schur_avail, wk_user, &wk_user_avail, colsca, &colsca_avail, rowsca, &rowsca_avail,
&(mumps_par->instance_number), &(mumps_par->nrhs), &(mumps_par->lrhs),
&(mumps_par->lredrhs),
rhs_sparse, &rhs_sparse_avail, sol_loc, &sol_loc_avail, irhs_sparse,
&irhs_sparse_avail, irhs_ptr, &irhs_ptr_avail, isol_loc,
&isol_loc_avail, &(mumps_par->nz_rhs), &(mumps_par->lsol_loc)
, &(mumps_par->schur_mloc)
, &(mumps_par->schur_nloc)
, &(mumps_par->schur_lld)
, &(mumps_par->mblock)
, &(mumps_par->nblock)
, &(mumps_par->nprow)
, &(mumps_par->npcol)
, ooc_tmpdir
, ooc_prefix
, write_problem
, &ooc_tmpdirlen
, &ooc_prefixlen
, &write_problemlen
);
/*
* mapping and pivnul_list are usually 0 except if
* MUMPS_ASSIGN_MAPPING/MUMPS_ASSIGN_PIVNUL_LIST was called.
*/
mumps_par->mapping=mumps_get_mapping();
mumps_par->pivnul_list=mumps_get_pivnul_list();
/* to get permutations computed during analysis */
mumps_par->sym_perm=mumps_get_sym_perm();
mumps_par->uns_perm=mumps_get_uns_perm();
/*
* colsca/rowsca can either be user data or have been modified
* within mumps by calls to MUMPS_ASSIGN_COLSCA and/or
* MUMPS_ASSIGN_ROWSCA. In all cases their address is contained
* in MUMPS_COLSCA_STATIC and/or MUMPS_ROWSCA_STATIC.
*
* In case of a null pointer, we also reset mumps_par->rowsca/colsca
* to 0 (case of JOB=-2, the Fortran pointer will be NULL but the
* C pointer should also be null.
*/
if (rowsca_avail == no) {
mumps_par->rowsca = MUMPS_ROWSCA_STATIC;
if (MUMPS_ROWSCA_STATIC) {
/* remember that row Scaling was computed by MUMPS */
mumps_par->rowsca_from_mumps=1;
}
}
if (colsca_avail == no) {
mumps_par->colsca = MUMPS_COLSCA_STATIC;
if (MUMPS_COLSCA_STATIC) {
/* remember that column Scaling was computed by MUMPS */
mumps_par->colsca_from_mumps=1;
}
}
}
