add alpha parameter for spmScal function

spm.c

@@ -1085,32 +1085,35 @@ spmCheckAxb( int nrhs,
  *
  * @ingroup pastix_spm
  *
- * @brief Scal a matrix with ||A||_2
+ * @brief Scal the spm: A = alpha * A
  *
  *******************************************************************************
  *
+ * @param[in] alpha
+ *           The scaling parameter.
+ *
  * @param[in,out] spm
  *          The sparse matrix to scal.
  *
  *******************************************************************************/
 void
-spmScal(pastix_spm_t* spm)
+spmScal(const pastix_complex64_t alpha, pastix_spm_t* spm)
 {
     switch(spm->flttype)
     {
     case PastixPattern:
         break;
     case PastixFloat:
-        s_spmScal(spm);
+        s_spmScal(alpha, spm);
         break;
     case PastixComplex32:
-        c_spmScal(spm);
+        c_spmScal(alpha, spm);
         break;
     case PastixComplex64:
-        z_spmScal(spm);
+        z_spmScal(alpha, spm);
         break;
     case PastixDouble:
     default:
-        d_spmScal(spm);
+        d_spmScal(alpha, spm);
     }
 }

spm.h

@@ -150,7 +150,7 @@ void *        spm2Dense( const pastix_spm_t *spm );
 pastix_int_t  spmFindBase( const pastix_spm_t *spm );
 double        spmNorm( int ntype, const pastix_spm_t *spm );
 int           spmMatVec(const pastix_trans_t trans, const void *alpha, const pastix_spm_t *spm, const void *x, const void *beta, void *y );
-void          spmScal( pastix_spm_t* spm );
+void          spmScal( const pastix_complex64_t alpha, pastix_spm_t* spm );
 
 int           spmSort( pastix_spm_t *spm );
 pastix_int_t  spmMergeDuplicate( pastix_spm_t *spm );

z_spm.h

@@ -66,7 +66,7 @@ void z_spmPrint( FILE *f, const pastix_spm_t *spm );
 
 pastix_spm_t *z_spmExpand(const pastix_spm_t *spm);
 void          z_spmDofExtend(pastix_spm_t *spm);
-void          z_spmScal( pastix_spm_t *spm );
+void          z_spmScal( const pastix_complex64_t alpha, pastix_spm_t *spm );
 
 
 #endif /* _z_spm_H_ */

z_spm_scal.c

@@ -24,26 +24,27 @@
  *
  * @ingroup pastix_spm_internal
  *
- * z_spmScal - Scal the matrix with ||A||_2
+ * z_spmScal - Scal the spm: A = alpha * A
  *
  *******************************************************************************
  *
+ * @param[in] alpha
+ *           The scaling parameter.
+ *
  * @param[in,out] spm
  *           The spm which needs to be scaled.
  *
  *******************************************************************************/
 void
-z_spmScal( pastix_spm_t *spm )
+z_spmScal( const pastix_complex64_t alpha, pastix_spm_t *spm )
 {
-    double              norm;
     pastix_int_t        nnz, i;
     pastix_complex64_t *values;
 
     nnz    = spm->nnz;
     values = spm->values;
-    norm   = z_spmNorm( PastixFrobeniusNorm, spm );
 
     for (i=0; i<nnz; i++){
-        values[i] /= norm;
+        values[i] *= alpha;
     }
 }