s_spm_norm.c File Reference

#include "common.h"
#include <lapacke.h>
#include <cblas.h>
#include "frobeniusupdate.h"

Go to the source code of this file.

#define	LAPACKE_slassq_work(_n_, _X_, _incx_, _scale_, _sumsq_)    __spm_slassq( (_n_), (_X_), (_incx_), (_scale_), (_sumsq_) )
	TODO.
static int	__spm_slassq (spm_int_t n, float *x, spm_int_t incx, float *scale, float *sumsq)
	Updates the values scale and sumsq such that.
static void	s_spm_frobenius_elt_sym_diag (spm_int_t dofs, const float *valptr, float *data)
	Compute the Frobenius norm of a diagonal element within a symmetric/symmetric matrix with column/row major storage.
static void	s_spm_frobenius_elt_sym_offd (spm_int_t nbelts, const float *valptr, float *data)
	Compute the Frobenius norm of an off-diagonal element matrix in the symmetric/symmetric case.
static void	s_spm_frobenius_elt_sym (spm_int_t row, spm_int_t dofi, spm_int_t col, spm_int_t dofj, const float *valptr, float *data)
	Compute the Frobenius norm of any element matrix in the symmetric/symmetric case.
static void	s_spmFrobeniusNorm_csc (const spmatrix_t *spm, float *data)
	Compute the Frobenius norm of a symmetrix/symmetric CSC matrix.
static void	s_spmFrobeniusNorm_csr (const spmatrix_t *spm, float *data)
	Compute the Frobenius norm of a symmetrix/symmetric CSR matrix.
static void	s_spmFrobeniusNorm_ijv (const spmatrix_t *spm, float *data)
	Compute the Frobenius norm of a symmetrix/symmetric IJV matrix.
float	s_spmFrobeniusNorm (const spmatrix_t *spm)
	Compute the Frobenius norm of the given spm structure.
float	s_spmMaxNorm (const spmatrix_t *spm)
	Compute the Max norm of the given spm structure.
static void	s_spm_oneinf_elt_sym_diag (spm_int_t row, spm_int_t dofi, const float *valptr, float *sumtab)
	Compute the sum array for the one/inf norms of a diagonal element within a symmetric/symmetric matrix with column/row major storage.
static void	s_spm_oneinf_elt_gen_A (spm_int_t dofi, spm_int_t dofj, const float *valptr, float *sumtab)
	Compute the sum array for the one/inf norms of a general element.
static void	s_spm_oneinf_elt_gen_B (spm_int_t dofi, spm_int_t dofj, const float *valptr, float *sumtab)
	Compute the sum array for the one/inf norms of a general element.
static void	s_spm_oneinf_elt_gen_AB (spm_int_t row, spm_int_t dofi, spm_int_t col, spm_int_t dofj, const float *valptr, float *sumtab)
	Compute the sum array for both the one and inf norms for the off-diagonal elements of symmetric/symmetric element matrices.
static void	s_spm_oneinf_elt_gen (spm_layout_t layout, spm_int_t row, spm_int_t dofi, spm_int_t col, spm_int_t dofj, const float *valptr, spm_normtype_t ntype, float *sumtab)
	Compute the sum array for the one and inf norms of a general element.
static void	s_spm_oneinf_elt_sym_offd (spm_layout_t layout, spm_int_t row, spm_int_t dofi, spm_int_t col, spm_int_t dofj, const float *valptr, float *sumtab)
	Compute the sum array for both the one and inf norms for the off-diagonal elements of symmetric/symmetric element matrices in either column or row major layout.
static void	s_spm_oneinf_elt (spm_mtxtype_t mtxtype, spm_layout_t layout, spm_int_t row, spm_int_t dofi, spm_int_t col, spm_int_t dofj, const float *valptr, spm_normtype_t ntype, float *sumtab)
	Compute the sum array for the one/inf norm for an element matrix.
static void	s_spmOneInfNorm_csc (spm_normtype_t ntype, const spmatrix_t *spm, float *sumtab)
	Compute the one/inf norm of an spm CSC structure.
static void	s_spmOneInfNorm_csr (spm_normtype_t ntype, const spmatrix_t *spm, float *sumtab)
	Compute the one/inf norm of an spm CSR structure.
static void	s_spmOneInfNorm_ijv (spm_normtype_t ntype, const spmatrix_t *spm, float *sumtab)
	Compute the one/inf norm of an spm IJV structure.
static float	s_spmOneInfNorm (spm_normtype_t ntype, const spmatrix_t *spm)
	Compute the one/inf norm of the given spm structure given by the maximum row sum.
float	s_spmNorm (spm_normtype_t ntype, const spmatrix_t *spm)
	Compute the norm of an spm matrix.
float	s_spmNormMat (spm_normtype_t ntype, const spmatrix_t *spm, spm_int_t n, const float *A, spm_int_t lda)
	Compute the norm of a dense matrix that follows the distribution of an spm matrix.

Detailed Description

SParse Matrix package norm routine.

Copyright

2016-2024 Bordeaux INP, CNRS (LaBRI UMR 5800), Inria, Univ. Bordeaux. All rights reserved.

Version

1.2.4

Author

Mathieu Faverge

Pierre Ramet

Tony Delarue

Matias Hastaran

Alycia Lisito

Date

2024-06-25

Generated arithmetic file from /builds/2mk6rsew/0/fpruvost/spm/src/z_spm_norm.c, normal z -> s, Fri Nov 29 11:34:30 2024

Definition in file s_spm_norm.c.

Macro Definition Documentation

LAPACKE_slassq_work

#define LAPACKE_slassq_work	(		_n_,
			_X_,
			_incx_,
			_scale_,
			_sumsq_
	)		__spm_slassq( (_n_), (_X_), (_incx_), (_scale_), (_sumsq_) )

TODO.

Definition at line 88 of file s_spm_norm.c.

Function Documentation

__spm_slassq()

static int __spm_slassq ( spm_int_t  n, float *  x, spm_int_t  incx, float *  scale, float *  sumsq  )

static

Updates the values scale and sumsq such that.

( scale**2 )*sumsq = x( 1 )**2 +...+ x( n )**2 + ( scale**2 )*sumsq,

This routine is inspired from LAPACK slassq function.

Parameters

[in]	n	The number of elements in the vector
[in]	x	The vector of size abs(n * incx)
[in]	incx	The increment between two elements in the vector x.
[in,out]	scale	On entry, the former scale On exit, the update scale to take into account the value
[in,out]	sumsq	On entry, the former sumsq On exit, the update sumsq to take into account the value

Returns

SPM_SUCESS to match the prototype of LAPACKE_slassq_work

Definition at line 63 of file s_spm_norm.c.

s_spm_frobenius_elt_sym_diag()

static void s_spm_frobenius_elt_sym_diag ( spm_int_t  dofs, const float *  valptr, float *  data  )

static

Compute the Frobenius norm of a diagonal element within a symmetric/symmetric matrix with column/row major storage.

Note that column major is using the low triangular part only of the diagonal element matrices, and row major, by symmetry, is using only the upper triangular part.

The comments in the code are made for column major storage.

Parameters

[in]	dofs	TODO
[in]	valptr	TODO
[in,out]	data	TODO

Definition at line 152 of file s_spm_norm.c.

Referenced by s_spm_frobenius_elt_sym().

s_spm_frobenius_elt_sym_offd()

static void s_spm_frobenius_elt_sym_offd ( spm_int_t  nbelts, const float *  valptr, float *  data  )

static

Compute the Frobenius norm of an off-diagonal element matrix in the symmetric/symmetric case.

Parameters

[in]	nbelts	TODO
[in]	valptr	TODO
[in,out]	data

Definition at line 206 of file s_spm_norm.c.

Referenced by s_spm_frobenius_elt_sym().

s_spm_frobenius_elt_sym()

static void s_spm_frobenius_elt_sym ( spm_int_t  row, spm_int_t  dofi, spm_int_t  col, spm_int_t  dofj, const float *  valptr, float *  data  )

static

Compute the Frobenius norm of any element matrix in the symmetric/symmetric case.

Parameters

[in]	row	TODO
[in]	dofi	TODO
[in]	col	TODO
[in]	dofj	TODO
[in]	valptr	TODO
[in,out]	data	TODO

Definition at line 252 of file s_spm_norm.c.

References s_spm_frobenius_elt_sym_diag(), and s_spm_frobenius_elt_sym_offd().

Referenced by s_spmFrobeniusNorm_csc(), s_spmFrobeniusNorm_csr(), and s_spmFrobeniusNorm_ijv().

s_spmFrobeniusNorm_csc()

static void s_spmFrobeniusNorm_csc ( const spmatrix_t *  spm, float *  data  )

static

Compute the Frobenius norm of a symmetrix/symmetric CSC matrix.

Parameters

[in]	spm	The spm from which the norm need to be computed.
[in,out]	data	TODO

Definition at line 283 of file s_spm_norm.c.

References spmatrix_s::baseval, spmatrix_s::colptr, spmatrix_s::dof, spmatrix_s::dofs, spmatrix_s::flttype, spmatrix_s::fmttype, spmatrix_s::loc2glob, spmatrix_s::n, spmatrix_s::replicated, spmatrix_s::rowptr, s_spm_frobenius_elt_sym(), SpmCSC, SpmFloat, and spmatrix_s::values.

Referenced by s_spmFrobeniusNorm().

s_spmFrobeniusNorm_csr()

static void s_spmFrobeniusNorm_csr ( const spmatrix_t *  spm, float *  data  )

static

Compute the Frobenius norm of a symmetrix/symmetric CSR matrix.

Parameters

[in]	spm	The spm from which the norm need to be computed.
[in,out]	data	TODO

Definition at line 350 of file s_spm_norm.c.

References spmatrix_s::baseval, spmatrix_s::colptr, spmatrix_s::dof, spmatrix_s::dofs, spmatrix_s::flttype, spmatrix_s::fmttype, spmatrix_s::loc2glob, spmatrix_s::n, spmatrix_s::replicated, spmatrix_s::rowptr, s_spm_frobenius_elt_sym(), SpmCSR, SpmFloat, and spmatrix_s::values.

Referenced by s_spmFrobeniusNorm().

s_spmFrobeniusNorm_ijv()

static void s_spmFrobeniusNorm_ijv ( const spmatrix_t *  spm, float *  data  )

static

Compute the Frobenius norm of a symmetrix/symmetric IJV matrix.

Parameters

[in]	spm	The spm from which the norm need to be computed.
[in,out]	data	TODO

Definition at line 418 of file s_spm_norm.c.

References spmatrix_s::baseval, spmatrix_s::colptr, spmatrix_s::dof, spmatrix_s::dofs, spmatrix_s::flttype, spmatrix_s::fmttype, spmatrix_s::nnz, spmatrix_s::rowptr, s_spm_frobenius_elt_sym(), SpmFloat, SpmIJV, and spmatrix_s::values.

Referenced by s_spmFrobeniusNorm().

s_spmFrobeniusNorm()

float s_spmFrobeniusNorm

(

const spmatrix_t *

spm

)

Compute the Frobenius norm of the given spm structure.

||A|| = sqrtf( sum( a_ij ^ 2 ) )

Parameters

[in]

spm

The spm from which the norm need to be computed.

Returns

The computed frobenius norm

Definition at line 480 of file s_spm_norm.c.

References spmatrix_s::clustnbr, spmatrix_s::comm, spmatrix_s::fmttype, spmatrix_s::mtxtype, spmatrix_s::nnzexp, spmatrix_s::replicated, s_spmFrobeniusNorm_csc(), s_spmFrobeniusNorm_csr(), s_spmFrobeniusNorm_ijv(), SpmCSC, SpmCSR, SpmGeneral, SpmIJV, and spmatrix_s::values.

Referenced by s_spmNorm().

s_spmMaxNorm()

float s_spmMaxNorm

(

const spmatrix_t *

spm

)

Compute the Max norm of the given spm structure.

||A|| = max( abs(a_ij) )

Parameters

[in]

spm

The spm from which the norm need to be computed.

Returns

The computed max norm

Definition at line 543 of file s_spm_norm.c.

References spmatrix_s::clustnbr, spmatrix_s::comm, spmatrix_s::nnzexp, spmatrix_s::replicated, and spmatrix_s::values.

Referenced by s_spmNorm().

s_spm_oneinf_elt_sym_diag()

static void s_spm_oneinf_elt_sym_diag ( spm_int_t  row, spm_int_t  dofi, const float *  valptr, float *  sumtab  )

static

Compute the sum array for the one/inf norms of a diagonal element within a symmetric/symmetric matrix with column/row major storage.

Note that column major is using the low triangular part only of the diagonal element matrices, and row major, by symmetry, is using only the upper triangular part.

The comments in the code are made for column major storage.

Parameters

[in]	row	TODO
[in]	dofi	TODO
[in]	valptr	TODO
[in,out]	sumtab	TODO

Definition at line 591 of file s_spm_norm.c.

Referenced by s_spm_oneinf_elt().

s_spm_oneinf_elt_gen_A()

static void s_spm_oneinf_elt_gen_A ( spm_int_t  dofi, spm_int_t  dofj, const float *  valptr, float *  sumtab  )

static

Compute the sum array for the one/inf norms of a general element.

We can observe two cases A and B;

---

| | | | | | |________| | | | | | | |A| OR |________| |B| |__|__|__| |_| |________| |_|

---

|___B____| |___A____|

          | One Norm | Inf norm |

----------—+-------—+-------—+ Column Major | B | A | ----------—+-------—+-------—+ Row Major | A | B | ----------—+-------—+-------—+

Warning

: The sumtab must be shifted at the right place on input

Parameters

[in]	dofi	TODO
[in]	dofj	TODO
[in]	valptr	TODO
[in,out]	sumtab	TODO

Definition at line 659 of file s_spm_norm.c.

Referenced by s_spm_oneinf_elt_gen().

s_spm_oneinf_elt_gen_B()

static void s_spm_oneinf_elt_gen_B ( spm_int_t  dofi, spm_int_t  dofj, const float *  valptr, float *  sumtab  )

static

Compute the sum array for the one/inf norms of a general element.

See s_spm_oneinf_elt_gen_A()

Parameters

[in]	dofi	TODO
[in]	dofj	TODO
[in]	valptr	TODO
[in,out]	sumtab	TODO

Definition at line 698 of file s_spm_norm.c.

Referenced by s_spm_oneinf_elt_gen().

s_spm_oneinf_elt_gen_AB()

static void s_spm_oneinf_elt_gen_AB ( spm_int_t  row, spm_int_t  dofi, spm_int_t  col, spm_int_t  dofj, const float *  valptr, float *  sumtab  )

static

Compute the sum array for both the one and inf norms for the off-diagonal elements of symmetric/symmetric element matrices.

See s_spm_oneinf_elt_gen_A()

Parameters

[in]	row	TODO
[in]	dofi	TODO
[in]	col	TODO
[in]	dofj	TODO
[in]	valptr	TODO
[in,out]	sumtab	TODO

Definition at line 744 of file s_spm_norm.c.

Referenced by s_spm_oneinf_elt_sym_offd().

s_spm_oneinf_elt_gen()

static void s_spm_oneinf_elt_gen ( spm_layout_t  layout, spm_int_t  row, spm_int_t  dofi, spm_int_t  col, spm_int_t  dofj, const float *  valptr, spm_normtype_t  ntype, float *  sumtab  )

static

Compute the sum array for the one and inf norms of a general element.

Parameters

[in]	layout	TODO
[in]	row	TODO
[in]	dofi	TODO
[in]	col	TODO
[in]	dofj	TODO
[in]	valptr	TODO
[in]	ntype	TODO
[in,out]	sumtab	TODO

Definition at line 799 of file s_spm_norm.c.

References s_spm_oneinf_elt_gen_A(), s_spm_oneinf_elt_gen_B(), SpmColMajor, SpmInfNorm, and SpmOneNorm.

Referenced by s_spm_oneinf_elt().

s_spm_oneinf_elt_sym_offd()

static void s_spm_oneinf_elt_sym_offd ( spm_layout_t  layout, spm_int_t  row, spm_int_t  dofi, spm_int_t  col, spm_int_t  dofj, const float *  valptr, float *  sumtab  )

static

Compute the sum array for both the one and inf norms for the off-diagonal elements of symmetric/symmetric element matrices in either column or row major layout.

Parameters

[in]	layout	TODO
[in]	row	TODO
[in]	dofi	TODO
[in]	col	TODO
[in]	dofj	TODO
[in]	valptr	TODO
[in,out]	sumtab	TODO

Definition at line 860 of file s_spm_norm.c.

References s_spm_oneinf_elt_gen_AB(), and SpmColMajor.

Referenced by s_spm_oneinf_elt().

s_spm_oneinf_elt()

static void s_spm_oneinf_elt ( spm_mtxtype_t  mtxtype, spm_layout_t  layout, spm_int_t  row, spm_int_t  dofi, spm_int_t  col, spm_int_t  dofj, const float *  valptr, spm_normtype_t  ntype, float *  sumtab  )

static

Compute the sum array for the one/inf norm for an element matrix.

Parameters

[in]	mtxtype	TODO
[in]	layout	TODO
[in]	row	TODO
[in]	dofi	TODO
[in]	col	TODO
[in]	dofj	TODO
[in]	valptr	TODO
[in]	ntype	TODO
[in,out]	sumtab	TODO

Definition at line 912 of file s_spm_norm.c.

References s_spm_oneinf_elt_gen(), s_spm_oneinf_elt_sym_diag(), s_spm_oneinf_elt_sym_offd(), and SpmGeneral.

Referenced by s_spmOneInfNorm_csc(), s_spmOneInfNorm_csr(), and s_spmOneInfNorm_ijv().

s_spmOneInfNorm_csc()

static void s_spmOneInfNorm_csc ( spm_normtype_t  ntype, const spmatrix_t *  spm, float *  sumtab  )

static

Compute the one/inf norm of an spm CSC structure.

Parameters

[in]	ntype	TODO
[in]	spm	TODO
[in,out]	sumtab	TODO

Definition at line 953 of file s_spm_norm.c.

References spmatrix_s::baseval, spmatrix_s::colptr, spmatrix_s::dof, spmatrix_s::dofs, spmatrix_s::layout, spmatrix_s::loc2glob, spmatrix_s::mtxtype, spmatrix_s::n, spmatrix_s::replicated, spmatrix_s::rowptr, s_spm_oneinf_elt(), and spmatrix_s::values.

Referenced by s_spmOneInfNorm().

s_spmOneInfNorm_csr()

static void s_spmOneInfNorm_csr ( spm_normtype_t  ntype, const spmatrix_t *  spm, float *  sumtab  )

static

Compute the one/inf norm of an spm CSR structure.

Parameters

[in]	ntype	TODO
[in]	spm	TODO
[in,out]	sumtab	TODO

Definition at line 1019 of file s_spm_norm.c.

References spmatrix_s::baseval, spmatrix_s::colptr, spmatrix_s::dof, spmatrix_s::dofs, spmatrix_s::layout, spmatrix_s::loc2glob, spmatrix_s::mtxtype, spmatrix_s::n, spmatrix_s::replicated, spmatrix_s::rowptr, s_spm_oneinf_elt(), and spmatrix_s::values.

Referenced by s_spmOneInfNorm().

s_spmOneInfNorm_ijv()

static void s_spmOneInfNorm_ijv ( spm_normtype_t  ntype, const spmatrix_t *  spm, float *  sumtab  )

static

Compute the one/inf norm of an spm IJV structure.

Parameters

[in]	ntype	TODO
[in]	spm	TODO
[in,out]	sumtab	TODO

Definition at line 1085 of file s_spm_norm.c.

References spmatrix_s::baseval, spmatrix_s::colptr, spmatrix_s::dof, spmatrix_s::dofs, spmatrix_s::layout, spmatrix_s::mtxtype, spmatrix_s::nnz, spmatrix_s::rowptr, s_spm_oneinf_elt(), and spmatrix_s::values.

Referenced by s_spmOneInfNorm().

s_spmOneInfNorm()

static float s_spmOneInfNorm ( spm_normtype_t  ntype, const spmatrix_t *  spm  )

static

Compute the one/inf norm of the given spm structure given by the maximum row sum.

• SpmOneNorm: ||A|| = max_j( sum_i(|a_ij|) )
• SpmInfNorm: ||A|| = max_i( sum_j(|a_ij|) )

Parameters

[in]	ntype	The type of norm to compute.
[in]	spm	The spm from which the norm need to be computed.

Returns

The computed one norm

Definition at line 1149 of file s_spm_norm.c.

References spmatrix_s::clustnbr, spmatrix_s::comm, spmatrix_s::fmttype, spmatrix_s::gNexp, spmatrix_s::replicated, s_spmOneInfNorm_csc(), s_spmOneInfNorm_csr(), s_spmOneInfNorm_ijv(), SpmCSC, SpmCSR, and SpmIJV.

Referenced by s_spmNorm().