include/ds/SymmetricMatrix.h

#ifndef _SYMMETRICMATRIX_H_
#define _SYMMETRICMATRIX_H_
// ---------------------------------------------
// company   : 
// user      : 
// date      : 22 Mar 2015 17:21
// file      : SymmetricMatrix.h
// markers   : Full
// language  : C++
// generator : Metamill 6.0
// ---------------------------------------------

//#UBLK-BEG-HEADERH
//==============================================================
// Description : 
//--------------------------------------------------------------
// Author      : Gilberto Luis
// Created     : 16 Dez 2010
//--------------------------------------------------------------
// Change history : 
//   16 Dez 2010 (Gilberto Luis) Initial version generated
//
//==============================================================
//#UBLK-END-HEADERH

/*=================================================================================================
=                                                                                                 =
=       Copyright 2010-2015  Marco Lucio Bittencourt / Fabiano Fernandes Bargos                       =
=                        Gilberto Luis Valente / Jorge Luis Suzuki / Allan Patrick Cordeiro Dias  =
=                                                                                                 =
=       Licensed under the Apache License, Version 2.0 (the "License");                               =
=       you may not use this file except in compliance with the License.                              =
=       You may obtain a copy of the License at                                                       =
=                                                                                                 =
=               http://www.apache.org/licenses/LICENSE-2.0                                                =
=                                                                                                 =
=       Unless required by applicable law or agreed to in writing, software                           =
=       distributed under the License is distributed on an "AS IS" BASIS,                             =
=       WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.                      =
=       See the License for the specific language governing permissions and                           =
=       limitations under the License.                                                                =
=                                                                                                 =
=================================================================================================*/


class Matrix;


#define SizeAttributesSM 28
#define SizeExtraMemSM (ALIGNMENT - (SizeAttributesSM % ALIGNMENT))


#include <cstdio>
#include "acdp/matematc.h"
#include "ds/DS_defs.h"
#include "ds/Matrix.h"


class SymmetricMatrix
{
//#UBLK-BEG-CLASSDEF mm:73e74e16-c4f4-11df-9da6-88402942916d
//#UBLK-END-CLASSDEF

 public:

   // METHOD=mm:7053b843-c4f4-11df-9da6-88402942916d
   SymmetricMatrix();

   // METHOD=mm:7053b848-c4f4-11df-9da6-88402942916d
   SymmetricMatrix(long MatrixOrder, double InitValue = 0.0);

   // METHOD=mm:7053b849-c4f4-11df-9da6-88402942916d
   SymmetricMatrix(FILE* File);

   // METHOD=mm:7053b84a-c4f4-11df-9da6-88402942916d
   SymmetricMatrix(SymmetricMatrix& Instance);

   // METHOD=mm:7053b84b-c4f4-11df-9da6-88402942916d
   ~SymmetricMatrix();

   // METHOD=mm:7053b84c-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& operator=(SymmetricMatrix& Instance);

   // METHOD=mm:7053b84d-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& Equal(SymmetricMatrix& Instance);

   // METHOD=mm:7053b852-c4f4-11df-9da6-88402942916d
   operator double* const();

   // METHOD=mm:7053b853-c4f4-11df-9da6-88402942916d
   double& operator()(long Row, long Col);

   // METHOD=mm:7053b854-c4f4-11df-9da6-88402942916d
   double& GetElement(unsigned long Row, unsigned long Column);

   // METHOD=mm:7053b855-c4f4-11df-9da6-88402942916d
   void SetOrder(unsigned long MatrixOrder = 0, double InitValue = 0.0);

   // METHOD=mm:7053b856-c4f4-11df-9da6-88402942916d
   unsigned long GetOrder();

   // METHOD=mm:7053b857-c4f4-11df-9da6-88402942916d
   unsigned long GetNumberElements();

   // METHOD=mm:7053b85c-c4f4-11df-9da6-88402942916d
   long SizeOf();

   // METHOD=mm:7053b85d-c4f4-11df-9da6-88402942916d
   void Reset(double Value = 0.0);

   // METHOD=mm:a041e7b5-8cfc-48d5-905a-b919489cbfa4
   void Reset(unsigned long NumEls, double Value = 0.0);

   // METHOD=mm:7053b85e-c4f4-11df-9da6-88402942916d
   void Identity();

   // METHOD=mm:7053b85f-c4f4-11df-9da6-88402942916d
   void Insert(SymmetricMatrix& SubMatrix, long InitLine = 0, long InitColumn = 0);

   // METHOD=mm:7053b860-c4f4-11df-9da6-88402942916d
   double Maximum();

   // METHOD=mm:7053b861-c4f4-11df-9da6-88402942916d
   double Minimum();

   // METHOD=mm:7053b866-c4f4-11df-9da6-88402942916d
   double EuclideanNorm();

   // METHOD=mm:5b142dff-c9c2-11e0-952b-da8bfa8c82f2
   void Product(double* A, unsigned long N, double b = 1.0);

   // METHOD=mm:7053b867-c4f4-11df-9da6-88402942916d
   double OneNorm();

   // METHOD=mm:7053b868-c4f4-11df-9da6-88402942916d
   double InfinityNorm();

   // METHOD=mm:7053b869-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& Chs();

   // METHOD=mm:7053b86a-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& operator-();

   // METHOD=mm:7053b86b-c4f4-11df-9da6-88402942916d
   friend SymmetricMatrix operator+(SymmetricMatrix& A, SymmetricMatrix& B);

   // METHOD=mm:7053b870-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& Add(SymmetricMatrix& A, SymmetricMatrix& B);

   // METHOD=mm:7053b871-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& operator+=(SymmetricMatrix& B);

   // METHOD=mm:7053b872-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& Add(SymmetricMatrix& B);

   // METHOD=mm:e6c436de-f1f2-4cf1-ad6b-66a470ea900d
   SymmetricMatrix& Add(SymmetricMatrix& B, unsigned long NumEls);

   // METHOD=mm:69efb8d2-cc70-4abc-9ac4-5cee706bd3d2
   void AddMultiplicate(SymmetricMatrix& B, double scalar, unsigned long NumEls);

   // METHOD=mm:7053b873-c4f4-11df-9da6-88402942916d
   friend SymmetricMatrix operator-(SymmetricMatrix& A, SymmetricMatrix& B);

   // METHOD=mm:7053b874-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& Subtract(SymmetricMatrix& A, SymmetricMatrix& B);

   // METHOD=mm:7053b875-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& operator-=(SymmetricMatrix& B);

   // METHOD=mm:7053b87a-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& Subtract(SymmetricMatrix& B);

   // METHOD=mm:7053b87b-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& Add(double Scalar);

   // METHOD=mm:7053b87c-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& operator+=(double Scalar);

   // METHOD=mm:7053b87d-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& Subtract(double Scalar);

   // METHOD=mm:7053b87e-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& operator-=(double Scalar);

   // METHOD=mm:7053b87f-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& operator*=(double Scalar);

   // METHOD=mm:7053b884-c4f4-11df-9da6-88402942916d
   SymmetricMatrix& Multiplicate(double Scalar);

   // METHOD=mm:7053b885-c4f4-11df-9da6-88402942916d
   void Multiplicate(double* X, double* Y);

   // METHOD=mm:c934d080-86f6-41af-9deb-7978ecaebec6
   void Multiplicate(double* X, unsigned long Size, double* Y);

   // METHOD=mm:036711fa-3a57-478a-91d3-91c641f7ef28
   void Multiplicate(double* A, unsigned long nlA, unsigned long ncA, double* B, unsigned long nlB, unsigned long ncB);

   // METHOD=mm:7053b886-c4f4-11df-9da6-88402942916d
   void TripleProduct(Matrix& A, SymmetricMatrix& B);

   // METHOD=mm:4463039a-b9cd-4329-9350-ad32a51189a9
   void TripleProduct(double* A, unsigned short NlinesA, unsigned short NcolsA, SymmetricMatrix& B);

   // METHOD=mm:7053b887-c4f4-11df-9da6-88402942916d
   double ScalarProduct(SymmetricMatrix& B);

   // METHOD=mm:7053b888-c4f4-11df-9da6-88402942916d
   double operator^(SymmetricMatrix& B);

   // METHOD=mm:7053b889-c4f4-11df-9da6-88402942916d
   void Save(FILE* File);

   // METHOD=mm:7053b88e-c4f4-11df-9da6-88402942916d
   void Save(char* Tabname, int Version, char* Filename);

   // METHOD=mm:7053b88f-c4f4-11df-9da6-88402942916d
   void Restore(FILE* File);

   // METHOD=mm:7053b890-c4f4-11df-9da6-88402942916d
   void Restore(char* Tabname, int Version, char* Filename);

   // METHOD=mm:7053b891-c4f4-11df-9da6-88402942916d
   void Print(FILE* File, char* Title, char* Format);

   // METHOD=mm:7053b892-c4f4-11df-9da6-88402942916d
   void Print(FILE* File = stdout, const char* Message = "", int NumMaxCols = 5, int Mant = 14, int Dec = 7);

   // METHOD=mm:7053b893-c4f4-11df-9da6-88402942916d
   void Print(FILE* File, int Mant, int Dec);

   // METHOD=mm:7053b898-c4f4-11df-9da6-88402942916d
   void Free();

   // METHOD=mm:7053b899-c4f4-11df-9da6-88402942916d
   void GaussSolver(double* U, double* B, double Precision);

   // METHOD=mm:7053b89a-c4f4-11df-9da6-88402942916d
   void GaussSolver();

   // METHOD=mm:7053b89b-c4f4-11df-9da6-88402942916d
   void GaussSolver(double* U, double* B);

   // METHOD=mm:7053b89c-c4f4-11df-9da6-88402942916d
   long Jacobi(double* U, double* B, long NumberIterations, double& NormRes, double Precision = 1e-4, ConvCriterion_E Crit = WEIGHT_RESIDUE, NormDS_E Norm = EUCLID);

   // METHOD=mm:914ecc0d-4659-4746-adac-481fc2560bc2
   void EigJacobi(SymmetricMatrix& B, double tol, SymmetricMatrix& eigv, Matrix& eigvec);

   // METHOD=mm:7053b89d-c4f4-11df-9da6-88402942916d
   long GaussSeidel(double* U, double* B, long NumberIterations, double& NormRes, double Precision = 1e-4, ConvCriterion_E Crit = WEIGHT_RESIDUE, NormDS_E Norm = EUCLID);

   // METHOD=mm:7053b8a2-c4f4-11df-9da6-88402942916d
   long SOR(double* U, double* B, long NumberIterations, double& NormRes, double Precision = 1e-4, double Omega = 1.81, ConvCriterion_E Crit = WEIGHT_RESIDUE, NormDS_E Norm = EUCLID);

   // METHOD=mm:7053b8a3-c4f4-11df-9da6-88402942916d
   long CG(double* U, double* B, long MaxNumberIter, double& NormRes, double Precision = 1e-4, ConvCriterion_E Crit = WEIGHT_RESIDUE, NormDS_E Norm = EUCLID);

   // METHOD=mm:7053b8a4-c4f4-11df-9da6-88402942916d
   long CGDiagonal(double* U, double* B, long MaxNumberIter, double& NormRes, double Precision = 1e-4, ConvCriterion_E Crit = WEIGHT_RESIDUE, NormDS_E Norm = EUCLID);

   // METHOD=mm:7053b8a5-c4f4-11df-9da6-88402942916d
   long CGPreConditionated(double* U, double* B, long MaxNumberIter, double& NormRes, double Precision = 1e-4, double Omega = 1.81, PreMatrix_E PreMatrix = SGS, ConvCriterion_E Crit = WEIGHT_RESIDUE, NormDS_E Norm = EUCLID);


 protected:

   // METHOD=mm:7053b836-c4f4-11df-9da6-88402942916d
   void FreeAuxVectors();

   // METHOD=mm:7053b837-c4f4-11df-9da6-88402942916d
   void Alloc(double InitValue = 0.0);

   // METHOD=mm:7053b838-c4f4-11df-9da6-88402942916d
   friend void Multiplicate(Matrix& A, SymmetricMatrix& B, Matrix& C);

   // METHOD=mm:7053b839-c4f4-11df-9da6-88402942916d
   int ConvCriteria(double* U, double* Un, double* b, double Prec, double Normb, ConvCriterion_E Crit, NormDS_E Norm, double& NormRes);

   // METHOD=mm:7053b83e-c4f4-11df-9da6-88402942916d
   void praux(double* vetb, double* vetx);

   // METHOD=mm:7053b83f-c4f4-11df-9da6-88402942916d
   void fwdsubst(double* vetx, double* vetb, double omega);

   // METHOD=mm:7053b840-c4f4-11df-9da6-88402942916d
   void prodaux(double* U, double* Y, double cte);

   // METHOD=mm:7053b841-c4f4-11df-9da6-88402942916d
   void bcksubst(double* vetx, double* vetb, double omega);

   // METHOD=mm:7053b842-c4f4-11df-9da6-88402942916d
   void prodfx(double* z, double* x, double Omega);

   // METHOD=mm:67847220-da91-41fc-ab0d-3004d8a01271
   double norinfv(double* u, long dim);

   // METHOD=mm:6c42e25e-9904-4829-b270-7ba46c2cdd49
   double scalv(double* x, double* y, unsigned long dim);

   // METHOD=mm:cc7567f7-633c-49e0-a669-e53b0c64df43
   double normaev(double* u, long dim);

   // METHOD=mm:c38ba932-9b75-4ff4-befd-6b96746f7b59
   int mgausimt(double* mas, long neq, long* step);

   // METHOD=mm:f1ad9249-1397-4d06-b372-2f1acb048b5b
   void pmatbsa(double* mcs, double* ma, long nlina, long ncola, double* mbs);

   // METHOD=mm:8273a4c4-56da-4d00-9901-2a6a9e1cd3dc
   void promabs(double* mc, double* ma, double* mbs, long nlina, long ncola);

   // METHOD=mm:45a0369c-f976-4e4b-ad59-0da92419c04e
   void solgasim(double* vetx, double* mas, double* vetb, long int neq);


   // attr=mm:7053b82b-c4f4-11df-9da6-88402942916d
   unsigned long PolyOrder;

   // attr=mm:7053b82e-c4f4-11df-9da6-88402942916d
   double* d;

   // attr=mm:7053b82f-c4f4-11df-9da6-88402942916d
   double* f;

   // attr=mm:7053b834-c4f4-11df-9da6-88402942916d
   double* g;

   // attr=mm:7053b835-c4f4-11df-9da6-88402942916d
   double* h;

   // attr=mm:7053b82d-c4f4-11df-9da6-88402942916d
   double* SMatrix;

   // attr=mm:7053b82c-c4f4-11df-9da6-88402942916d
   DoublePtr_T* PSLines;

   // attr=mm:54b3d705-3866-11e0-893c-fc4410c867b1
   char ExtraMemory[SizeExtraMemSM];


 private:

   // METHOD=mm:639488c7-38f3-11e0-8edb-aca828ea8a14
   long SizeOfAttributes();


//#UBLK-BEG-CLASSMEMB mm:73e74e16-c4f4-11df-9da6-88402942916d
//#UBLK-END-CLASSMEMB
};

//#UBLK-BEG-GLOBALH
//*** Removed methods: 17 Jan 2012 15:08
//---
//uuid=mm:51a4cda9-4dc8-49d1-8ca4-090526094cc9
//   void prmatba(double* ma, double* mb, long nlina, long ncola, double* mc);

//*** Removed methods: 24 Jan 2012 17:02
//---
//uuid=mm:eaed1dd4-a20e-4bdd-be5f-da6ef24abb2a
//   friend Matrix& operator*(Matrix& A, SymmetricMatrix& B);

//#UBLK-END-GLOBALH


#endif // _SYMMETRICMATRIX_H_