#include <math/scmat/repl.h>

Inheritance diagram for sc::ReplSymmSCMatrix:

Public Member Functions
	ReplSymmSCMatrix (const RefSCDimension &, ReplSCMatrixKit *)

void	assign_val (double)
	Overridden to implement the assign functions.

void	assign_s (SymmSCMatrix *)

void	assign_p (const double *)

void	assign_pp (const double **)

void	convert_p (double *) const

void	convert_pp (double **) const

double	get_element (int, int) const
	Return or modify an element.

void	set_element (int, int, double)

void	accumulate_element (int, int, double)

void	scale (double)
	Multiply all elements by val.

SCMatrix *	get_subblock (int, int, int, int)
	Return a subblock of this. More...

SymmSCMatrix *	get_subblock (int, int)

void	assign_subblock (SCMatrix *, int, int, int, int)
	Assign m to a subblock of this.

void	assign_subblock (SymmSCMatrix *, int, int)

void	accumulate_subblock (SCMatrix *, int, int, int, int)
	Sum m into a subblock of this.

void	accumulate_subblock (SymmSCMatrix *, int, int)

SCVector *	get_row (int i)
	Return a row of this.

void	assign_row (SCVector *v, int i)
	Assign v to a row of this.

void	accumulate_row (SCVector *v, int i)
	Sum v to a row of this.

void	accumulate_product_rr (SCMatrix , SCMatrix )

void	accumulate (const SymmSCMatrix *)
	Sum m into this.

double	invert_this ()
	Invert this.

double	solve_this (SCVector *)

double	trace ()
	Return the trace.

double	determ_this ()
	Return the determinant of this. this is overwritten.

void	gen_invert_this (double condition_number_threshold=1e8)
	Return the generalized inverse of this using SVD decomposition. More...

double	scalar_product (SCVector *)
	Return the scalar obtained by multiplying this on the left and right by v.

void	diagonalize (DiagSCMatrix , SCMatrix )
	Diagonalize this, placing the eigenvalues in d and the eigenvectors in m.

void	eigensystem (SymmSCMatrix , DiagSCMatrix , SCMatrix *)
	Solve generalized eigensystem for this with metric s, placing the eigenvalues in d and the eigenvectors in m. More...

void	accumulate_symmetric_outer_product (SCVector *)

void	accumulate_symmetric_product (SCMatrix *)
	Sum into this a * a.t()

void	accumulate_symmetric_sum (SCMatrix *)
	Sum into a + a.t()

void	accumulate_transform (SCMatrix , SymmSCMatrix , SCMatrix::Transform=SCMatrix::NormalTransform)

void	accumulate_transform (SCMatrix , DiagSCMatrix , SCMatrix::Transform=SCMatrix::NormalTransform)

void	accumulate_transform (SymmSCMatrix , SymmSCMatrix )

void	element_op (const Ref< SCElementOp > &)
	only applied to the unique elements of this. More...

void	element_op (const Ref< SCElementOp2 > &, SymmSCMatrix *)

void	element_op (const Ref< SCElementOp3 > &, SymmSCMatrix , SymmSCMatrix )

void	vprint (const char *title=0, std::ostream &out=ExEnv::out0(), int=10) const

double *	get_data ()

double **	get_rows ()

Ref< SCMatrixSubblockIter >	local_blocks (SCMatrixSubblockIter::Access)
	Returns iterators for the local (rapidly accessible) blocks used in this matrix.

Ref< SCMatrixSubblockIter >	all_blocks (SCMatrixSubblockIter::Access)
	Returns iterators for the all blocks used in this matrix.

Ref< ReplSCMatrixKit >	skit ()

Public Member Functions inherited from sc::SymmSCMatrix
	SymmSCMatrix (const RefSCDimension &, SCMatrixKit *)

Ref< SCMatrixKit >	kit () const
	Return the SCMatrixKit object that created this object.

virtual void	save (StateOut &)
	Save and restore this in an implementation independent way.

virtual void	restore (StateIn &)

virtual double	maxabs () const
	Return the maximum absolute value element of this vector.

virtual void	randomize ()
	Assign each element to a random number between -1 and 1.

void	assign (double val)
	Set all elements to val.

void	assign (const double *m)
	Assign element i, j to m[i*(i+1)/2+j].

void	assign (const double **m)
	Assign element i, j to m[i][j].

void	assign (SymmSCMatrix *m)
	Make this have the same elements as m. More...

void	convert (double *a) const
	Like the assign members, but these write values to the arguments.

void	convert (double **a) const

void	convert (SymmSCMatrix *a)

virtual void	convert_s (SymmSCMatrix *)

virtual void	convert_accumulate (SymmSCMatrix *)

virtual void	scale_diagonal (double)
	Scale the diagonal elements by val.

virtual void	shift_diagonal (double)
	Shift the diagonal elements by val.

virtual void	unit ()
	Make this equal to the unit matrix.

int	n () const
	Return the dimension.

virtual SymmSCMatrix *	copy ()
	Return a matrix with the same dimension and same elements.

virtual SymmSCMatrix *	clone ()
	Return a matrix with the same dimension but uninitialized memory.

RefSCDimension	dim () const
	Return the dimension.

void	print (std::ostream &o=ExEnv::out0()) const
	Print out the matrix.

void	print (const char *title=0, std::ostream &out=ExEnv::out0(), int=10) const

Ref< MessageGrp >	messagegrp () const
	Returns the message group used by the matrix kit.

Public Member Functions inherited from sc::DescribedClass
	DescribedClass (const DescribedClass &)

DescribedClass &	operator= (const DescribedClass &)

ClassDesc *	class_desc () const MPQC__NOEXCEPT
	This returns the unique pointer to the ClassDesc corresponding to the given type_info object. More...

const char *	class_name () const
	Return the name of the object's exact type.

int	class_version () const
	Return the version of the class.

Ref< DescribedClass >	ref ()
	Return this object wrapped up in a Ref smart pointer. More...

Public Member Functions inherited from sc::RefCount
size_t	identifier () const
	Return the unique identifier for this object that can be compared for different objects of different types. More...

int	lock_ptr () const
	Lock this object.

int	unlock_ptr () const
	Unlock this object.

void	use_locks (bool inVal)
	start and stop using locks on this object

refcount_t	nreference () const
	Return the reference count.

refcount_t	reference ()
	Increment the reference count and return the new count.

refcount_t	dereference ()
	Decrement the reference count and return the new count.

int	managed () const

void	unmanage ()
	Turn off the reference counting mechanism for this object. More...

Protected Member Functions
size_t	compute_offset (int, int) const

void	init_blocklist ()

void	before_elemop ()

void	after_elemop ()

Protected Member Functions inherited from sc::RefCount
	RefCount (const RefCount &)

RefCount &	operator= (const RefCount &)

Protected Attributes
Ref< SCMatrixBlockList >	blocklist

double *	matrix

double **	rows

Protected Attributes inherited from sc::SymmSCMatrix
RefSCDimension	d

Ref< SCMatrixKit >	kit_

Friends
class	ReplSCMatrix

class	ReplDiagSCMatrix

class	ReplSCVector

Detailed Description

Replicated SymmSCMatrix.

Member Function Documentation

◆ eigensystem()

void sc::ReplSymmSCMatrix::eigensystem	(	SymmSCMatrix *	s,
		DiagSCMatrix *	d,
		SCMatrix *	m
	)

virtual

Solve generalized eigensystem for this with metric s, placing the eigenvalues in d and the eigenvectors in m.

Exceptions

AlgorithmException eigensystem could not be solved within available precision. Check condition number of s.

Implements sc::SymmSCMatrix.

◆ element_op()

void sc::ReplSymmSCMatrix::element_op ( const Ref< SCElementOp > & )

virtual

only applied to the unique elements of this.

For example, to compute the sum of all matrix elements you need to do the following:

SymmSCMatrix* A;   // presumed initialized elsewhere
A->element_op(new SCElementScaleDiagonal(0.5));   // scale the diagonal by 1/2
SCElementSum* sum_op = new SCElementSum;
A->element_op(sum_op);
std::cout << "Sum of element of matrix A = " << sum_op->result() * 2.0 << std::endl;

Implements sc::SymmSCMatrix.