Changeset 758 in Sophya for trunk/SophyaLib/Samba/alm.h

Timestamp:

Mar 2, 2000, 5:25:35 PM (26 years ago)

Author:

ansari

Message:

Restructuration de Sophya en petits modules (creation de SkyMap en particulier)

Reza 2/3/2000

File:

• trunk/SophyaLib/Samba/alm.h (modified) (1 diff)

trunk/SophyaLib/Samba/alm.h

@@ -1 +1 @@
 #ifndef ALM_SEEN
 #define ALM_SEEN
+
 #include "nbrandom.h"
 #include "nbmath.h"
+#include "triangmtx.h"
 #include "tvector.h"
-#include "pexceptions.h"
-/*! Class for inferior triangular matrix (base class for the class Alm) */
-template <class T>
-class TriangularMatrix 
-  {
-    
-  public :
-
-TriangularMatrix()   {};
-  /* instanciate a triangular matrix from the number of rows */
-TriangularMatrix(int rowSize)  : long_diag_((uint_4)rowSize) {elem_.ReSize((uint_4) (rowSize*(rowSize+1)/2) ); };
-TriangularMatrix(const TriangularMatrix<T>& a,  bool share=false)  : elem_(a.elem_, share),  long_diag_(a.long_diag_) {;}
-/*! resize the matrix with a new number of rows */
-inline void ReSizeRow(int rowSize) 
-  {
-    long_diag_=(uint_4)rowSize;
-    elem_.ReSize(long_diag_*(long_diag_+1)/2);
-  }
- inline void SetTemp(bool temp=false) const {elem_.SetTemp(temp);}
-
-inline TriangularMatrix<T>& operator = (const TriangularMatrix<T>& a)
-{
-  elem_=a.elem_;
-  long_diag_ = a.long_diag_;
-  return *this;
-}
-inline T& operator()(int l, int m) 
-  {
-      return  elem_(adr_ij(l,m));
-  }
-inline T const& operator()(int l, int m) const 
-  {
-      return *(elem_.Begin()+ adr_ij(l,m));
-  }
- inline  int_4  rowNumber() const {return (int_4)long_diag_;}
-  private: 
- /*! compute the address of an element in the single array representing the matrix */
-inline uint_4 adr_ij(int i,int j) const 
-{
-  int adr= i*(i+1)/2+j;
-  if ( adr >= elem_.Size() || adr <0 )
-{
- cout << " attention depassement dans triangularMatrix " << endl;
- cout << " l= " << i << " m= " << j << " tableau reserve longueur " << elem_.Size() << endl; 
-}
-  return(i*(i+1)/2+j);
-}
-
- uint_4 long_diag_;
- NDataBlock<T> elem_;
-
-  };