From 0ecbe60bff10060e0d3c0e251903db0f139f6b11 Mon Sep 17 00:00:00 2001
From: Mathieu Faverge <mathieu.faverge@inria.fr>
Date: Fri, 11 Nov 2016 13:37:05 +0100
Subject: [PATCH] Start getting some norms to work in cmpressed csc
---
z_spm_norm.c | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/z_spm_norm.c b/z_spm_norm.c
index be0ed77f..7ce9d3e3 100644
--- a/z_spm_norm.c
+++ b/z_spm_norm.c
@@ -50,7 +50,7 @@ z_spmFrobeniusNorm( const pastix_spm_t *spm )
double sumsq = 0.;
if (spm->mtxtype == PastixGeneral) {
- for(i=0; i <spm->nnz; i++, valptr++) {
+ for(i=0; i <spm->nnzexp; i++, valptr++) {
frobenius_update( 1, &scale, &sumsq, valptr );
#if defined(PRECISION_z) || defined(PRECISION_c)
@@ -139,7 +139,7 @@ z_spmMaxNorm( const pastix_spm_t *spm )
pastix_complex64_t *valptr = (pastix_complex64_t*)spm->values;
double tmp, norm = 0.;
- for(i=0; i <spm->nnz; i++, valptr++) {
+ for(i=0; i <spm->nnzexp; i++, valptr++) {
tmp = cabs( *valptr );
norm = norm > tmp ? norm : tmp;
}
--
GitLab