deprecated register keyword removed from bnolib

bnolib/DDiagMatrix.cpp

@@ -27,7 +27,7 @@ DDiagMatrix& DDiagMatrix::operator+=(const DDiagMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Vector != NULL);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < Dimension; i++)
         Vector[i] += mt.Vector[i];
@@ -44,7 +44,7 @@ DDiagMatrix& DDiagMatrix::operator-=(const DDiagMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Vector != NULL);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < Dimension; i++)
         Vector[i] -= mt.Vector[i];
@@ -62,7 +62,7 @@ DDiagMatrix& DDiagMatrix::MulRight(const DDiagMatrix& mt)
     TAssert(Vector != NULL);
     TAssert(Vector != mt.Vector);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < Dimension; i++)
         Vector[i] *= mt.Vector[i];
@@ -77,7 +77,7 @@ bool DDiagMatrix::Expand(DUppMatrix& s, DDiagMatrix& d) const
 {
     TAssert(Vector != NULL);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < Dimension; i++)
         if (fabs(Vector[i]) < ZERO_ACCURACY)
@@ -100,7 +100,7 @@ bool DDiagMatrix::Inverse()
 {
     TAssert(Vector != NULL);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         if (fabs(Vector[i]) < ZERO_ACCURACY) {
@@ -145,7 +145,7 @@ void DDiagMatrix::GetArray(double* array) const
 
     memset(array, 0, Dimension * Dimension * sizeof(double));
 
-    register unsigned i;
+    unsigned i;
     for (i = 0; i < Dimension; i++) {
         *array = Vector[i];
         array += Dimension + 1;

bnolib/DLineMem.cpp

@@ -88,7 +88,7 @@ void DLineMem::Add(const CDMatrix& mt)
     TAssert(mt.Cols() == 1);
     TAssert(Dimension == mt.Rows());
 
-    register unsigned i;
+    unsigned i;
     for (i = 0; i < mt.Rows(); i++)
         Vector[i] += mt(i, 0);
 
@@ -104,7 +104,7 @@ void DLineMem::Sub(const CDMatrix& mt)
     TAssert(mt.Cols() == 1);
     TAssert(Dimension == mt.Rows());
 
-    register unsigned i;
+    unsigned i;
     for (i = 0; i < mt.Rows(); i++)
         Vector[i] -= mt(i, 0);
 
@@ -120,7 +120,7 @@ void DLineMem::AddVect(const DLineMem& mt)
     TAssert(mt.Vector != NULL);
     TAssert(Dimension == mt.Dimension);
 
-    register unsigned i;
+    unsigned i;
     for (i = 0; i < mt.Dimension; i++)
         Vector[i] += mt.Vector[i];
 
@@ -136,7 +136,7 @@ void DLineMem::SubVect(const DLineMem& mt)
     TAssert(mt.Vector != NULL);
     TAssert(Dimension == mt.Dimension);
 
-    register unsigned i;
+    unsigned i;
     for (i = 0; i < mt.Dimension; i++)
         Vector[i] -= mt.Vector[i];
 
@@ -163,7 +163,7 @@ double DLineMem::GetMaxElementAbs() const
     double max_value = fabs(Vector[0]);
     double value_abs;
 
-    for (register unsigned i = 1; i < Dimension; ++i) {
+    for (unsigned i = 1; i < Dimension; ++i) {
         value_abs = fabs(Vector[i]);
         if (value_abs > max_value)
             max_value = value_abs;
@@ -183,7 +183,7 @@ double DLineMem::GetMinElementAbs() const
     double min_value = fabs(Vector[0]);
     double value_abs;
 
-    for (register unsigned i = 1; i < Dimension; ++i) {
+    for (unsigned i = 1; i < Dimension; ++i) {
         value_abs = fabs(Vector[i]);
         if (value_abs < min_value)
             min_value = value_abs;
@@ -202,7 +202,7 @@ double DLineMem::GetMaxElement() const
 
     double max_value = Vector[0];
 
-    for (register unsigned i = 1; i < Dimension; ++i) {
+    for (unsigned i = 1; i < Dimension; ++i) {
         max_value = (Vector[i] > max_value) ? Vector[i] : max_value;
     }
 
@@ -219,7 +219,7 @@ double DLineMem::GetMinElement() const
 
     double min_value = Vector[0];
 
-    for (register unsigned i = 1; i < Dimension; ++i) {
+    for (unsigned i = 1; i < Dimension; ++i) {
         min_value = (Vector[i] < min_value) ? Vector[i] : min_value;
     }
 
@@ -233,7 +233,7 @@ void DLineMem::Init(double value)
 {
     TAssert(Vector != NULL);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < Dimension; i++)
         Vector[i] = value;
@@ -246,7 +246,7 @@ void DLineMem::Mult(double mltp)
 {
     TAssert(Vector != NULL);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < Dimension; i++)
         Vector[i] *= mltp;
@@ -267,7 +267,7 @@ void DLineMem::MulRect(const DRectMem& mt, const DLineMem& vc)
     if (Vector == NULL)
         return;
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < mt.NRows; i++) {
@@ -295,7 +295,7 @@ void DLineMem::MulLeftRect(const DRectMem& mt)
     if (Vector == NULL)
         return;
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < mt.NRows; i++) {

bnolib/DLowMatrix.cpp

@@ -13,7 +13,7 @@ void DLowMatrix::Init(const DLowMatrix& mt)
     TAssert(mt.Buffer != NULL);
     DTriagMem::Init(mt.Dim());
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = 0; j < i + 1; j++)
@@ -29,7 +29,7 @@ void DLowMatrix::Init(const DDiagMatrix& mt)
     DTriagMem::Init(mt.Dim());
 
     DTriagMem::Init(0.0);
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < Dimension; i++)
         Buffer[Index(i, i)] = mt.Vector[i];
@@ -44,7 +44,7 @@ DLowMatrix& DLowMatrix::operator+=(const DLowMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Buffer != NULL);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = 0; j < i + 1; j++)
@@ -63,7 +63,7 @@ DLowMatrix& DLowMatrix::operator-=(const DLowMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Buffer != NULL);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = 0; j < i + 1; j++)
@@ -81,7 +81,7 @@ DLowMatrix& DLowMatrix::operator+=(const DDiagMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Buffer != NULL);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < Dimension; i++)
         Buffer[Index(i, i)] += mt.Vector[i];
@@ -97,7 +97,7 @@ DLowMatrix& DLowMatrix::operator-=(const DDiagMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Buffer != NULL);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < Dimension; i++)
         Buffer[Index(i, i)] -= mt.Vector[i];
@@ -118,7 +118,7 @@ DLowMatrix& DLowMatrix::MulRight(const DLowMatrix& mt)
     DLowMatrix src(Dimension);
     XChng(src);
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
     double* p;
 
     for (i = 0; i < Dimension; i++) {
@@ -146,7 +146,7 @@ DLowMatrix& DLowMatrix::MulLeft(const DLowMatrix& mt)
     DLowMatrix src(Dimension);
     XChng(src);
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
     double* p;
 
     for (i = 0; i < Dimension; i++) {
@@ -173,7 +173,7 @@ DLowMatrix& DLowMatrix::MulRight(const DDiagMatrix& mt)
     DLowMatrix src(Dimension);
     XChng(src);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = 0; j < i + 1; j++)
@@ -195,7 +195,7 @@ DLowMatrix& DLowMatrix::MulLeft(const DDiagMatrix& mt)
     DLowMatrix src(Dimension);
     XChng(src);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = 0; j < i + 1; j++)
@@ -214,7 +214,7 @@ DLowMatrix& DLowMatrix::Transpose(const DUppMatrix& mt)
 
     DTriagMem::Init(mt.Dim());
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++)
         for (j = 0; j < i + 1; j++)
@@ -246,7 +246,7 @@ void DLowMatrix::GetLine(unsigned index, double* col) const
     TAssert(col != NULL);
     TAssert(index < Dimension);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i <= index; i++)
         col[i] = *(Buffer + Index(i, index));
@@ -262,7 +262,7 @@ void DLowMatrix::GetArray(double* array) const
     TAssert(Buffer != NULL);
     TAssert(array != NULL);
 
-    register unsigned i;
+    unsigned i;
     for (i = 0; i < Dimension; i++) {
         GetLine(i, array);
         array += Dimension;

bnolib/DRectMem.cpp

@@ -85,7 +85,7 @@ void DRectMem::SetLinesPtrs(double** ptr, unsigned n, unsigned m)
 {
     double* d_point = (double *)(ptr + n);
 
-    register unsigned i;
+    unsigned i;
     for (i = 0; i < n; ++i, d_point += m)
         ptr[i] = d_point;
 }
@@ -131,7 +131,7 @@ void DRectMem::Init(double value)
 {
     TAssert(Matrix != NULL);
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -150,7 +150,7 @@ void DRectMem::Init(const CDMatrix& mt)
     if (Matrix == NULL)
         return;
 
-    register unsigned i;
+    unsigned i;
     double* line;
 
     line = new double[NCols];
@@ -173,7 +173,7 @@ void DRectMem::InitRect(const DRectMem& mt)
     if (Matrix == NULL)
         return;
 
-    register unsigned i;
+    unsigned i;
 
     if (NRows == NRowsMem && NCols == NColsMem) {
         memcpy(Matrix[0], mt.Matrix[0], sizeof(double) * NRows * NCols);
@@ -193,7 +193,7 @@ void DRectMem::Init(const DSymmMatrix& mt)
     if (Matrix == NULL)
         return;
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -215,7 +215,7 @@ void DRectMem::Init(const DUppMatrix& mt)
     if (Matrix == NULL)
         return;
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -235,7 +235,7 @@ void DRectMem::Init(const DLowMatrix& mt)
     if (Matrix == NULL)
         return;
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -255,7 +255,7 @@ void DRectMem::Init(const DDiagMatrix& mt)
     if (Matrix == NULL)
         return;
 
-    register unsigned i;
+    unsigned i;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -272,7 +272,7 @@ void DRectMem::Mult(double mltp)
 {
     TAssert(Matrix != NULL);
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -291,7 +291,7 @@ void DRectMem::Add(const CDMatrix& mt)
     TAssert(NRows == mt.Rows());
     TAssert(NCols == mt.Cols());
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -311,7 +311,7 @@ void DRectMem::Sub(const CDMatrix& mt)
     TAssert(NRows == mt.Rows());
     TAssert(NCols == mt.Cols());
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -331,7 +331,7 @@ void DRectMem::AddRect(const DRectMem& mt)
     TAssert(NRows == mt.NRows);
     TAssert(NCols == mt.NCols);
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -351,7 +351,7 @@ void DRectMem::SubRect(const DRectMem& mt)
     TAssert(NRows == mt.NRows);
     TAssert(NCols == mt.NCols);
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -371,7 +371,7 @@ void DRectMem::Add(const DSymmMatrix& mt)
     TAssert(NRows == mt.Dim());
     TAssert(NCols == mt.Dim());
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -396,7 +396,7 @@ void DRectMem::Sub(const DSymmMatrix& mt)
     TAssert(NRows == mt.Dim());
     TAssert(NCols == mt.Dim());
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -421,7 +421,7 @@ void DRectMem::Add(const DUppMatrix& mt)
     TAssert(NRows == mt.Dim());
     TAssert(NCols == mt.Dim());
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -442,7 +442,7 @@ void DRectMem::Sub(const DUppMatrix& mt)
     TAssert(NRows == mt.Dim());
     TAssert(NCols == mt.Dim());
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -463,7 +463,7 @@ void DRectMem::Add(const DLowMatrix& mt)
     TAssert(NRows == mt.Dim());
     TAssert(NCols == mt.Dim());
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -484,7 +484,7 @@ void DRectMem::Sub(const DLowMatrix& mt)
     TAssert(NRows == mt.Dim());
     TAssert(NCols == mt.Dim());
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -505,7 +505,7 @@ void DRectMem::Add(const DDiagMatrix& mt)
     TAssert(NRows == mt.Dim());
     TAssert(NCols == mt.Dim());
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < NRows; i++)
         Matrix[i][i] += mt.Vector[i];
@@ -522,7 +522,7 @@ void DRectMem::Sub(const DDiagMatrix& mt)
     TAssert(NRows == mt.Dim());
     TAssert(NCols == mt.Dim());
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < NRows; i++)
         Matrix[i][i] -= mt.Vector[i];
@@ -545,7 +545,7 @@ void DRectMem::Mult(const DRectMem& m1, const DRectMem& m2)
     if (Matrix == NULL)
         return;
 
-    register unsigned i, j, k, f;
+    unsigned i, j, k, f;
     double *p, *line, *lines;
     f = m2.NRows;
 
@@ -618,7 +618,7 @@ void DRectMem::Mult(const DRectMem& m, const DUppMatrix& mt)
     if (Matrix == NULL)
         return;
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
     double *p, *line, *lines;
 
     for (i = 0; i < NRows; i++) {
@@ -666,7 +666,7 @@ void DRectMem::Mult(const DUppMatrix& mt, const DRectMem& m)
     if (Matrix == NULL)
         return;
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
     double *p, *line;
 
     for (i = 0; i < mt.Dimension; i++) {
@@ -712,7 +712,7 @@ void DRectMem::Mult(const DRectMem& m, const DLowMatrix& mt)
     if (Matrix == NULL)
         return;
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
     double *p, *line;
 
     for (i = 0; i < NRows; i++) {
@@ -759,7 +759,7 @@ void DRectMem::Mult(const DLowMatrix& mt, const DRectMem& m)
     if (Matrix == NULL)
         return;
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
     double *p, *line;
 
     for (i = 0; i < mt.Dimension; i++) {
@@ -806,7 +806,7 @@ void DRectMem::Mult(const DRectMem& m, const DDiagMatrix& mt)
     if (Matrix == NULL)
         return;
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NRows; i++) {
@@ -849,7 +849,7 @@ void DRectMem::Mult(const DDiagMatrix& mt, const DRectMem& m)
     if (Matrix == NULL)
         return;
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < mt.Dimension; i++) {

bnolib/DSqrMatrix.cpp

@@ -15,7 +15,7 @@ void DSqrMatrix::Transpose()
     TAssert(Matrix != NULL);
     TAssert(NRows == NCols);
 
-    register unsigned i, j;
+    unsigned i, j;
     double a;
     double* line;
 
@@ -123,7 +123,7 @@ void DSqrMatrix::GetArray(double* array) const
     TAssert(Matrix != NULL);
     TAssert(array != NULL);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < NRows; i++) {
         memcpy(array, Matrix[i], NCols * sizeof(double));
@@ -140,7 +140,7 @@ void DSqrMatrix::GetCol(unsigned index, double* col) const
     TAssert(index < NCols);
     TAssert(col != NULL);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < NRows; i++)
         col[i] = Matrix[i][index];

bnolib/DSymmMatrix.cpp

@@ -17,7 +17,7 @@ void DSymmMatrix::Init(const DDiagMatrix& mt)
     if (Buffer == NULL)
         return;
 
-    register unsigned i;
+    unsigned i;
 
     DTriagMem::Init(0.0);
     for (i = 0; i < Dimension; i++)
@@ -33,7 +33,7 @@ void DSymmMatrix::Init(const CDSymmMatrix& mt)
     if (Buffer == NULL)
         return;
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = i; j < Dimension; j++)
@@ -49,7 +49,7 @@ void DSymmMatrix::Init(const DSymmMatrix& mt)
     if (Buffer == NULL)
         return;
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = i; j < Dimension; j++)
@@ -64,7 +64,7 @@ DSymmMatrix& DSymmMatrix::operator+=(const CDSymmMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Buffer != NULL);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = i; j < Dimension; j++)
@@ -83,7 +83,7 @@ DSymmMatrix& DSymmMatrix::operator+=(const DSymmMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Buffer != NULL);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = i; j < Dimension; j++)
@@ -101,7 +101,7 @@ DSymmMatrix& DSymmMatrix::operator-=(const CDSymmMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Buffer != NULL);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = i; j < Dimension; j++)
@@ -120,7 +120,7 @@ DSymmMatrix& DSymmMatrix::operator-=(const DSymmMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Buffer != NULL);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = i; j < Dimension; j++)
@@ -138,7 +138,7 @@ DSymmMatrix& DSymmMatrix::operator+=(const DDiagMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Buffer != NULL);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < Dimension; i++)
         Buffer[Index(i, i)] += mt.Vector[i];
@@ -154,7 +154,7 @@ DSymmMatrix& DSymmMatrix::operator-=(const DDiagMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Buffer != NULL);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < Dimension; i++)
         Buffer[Index(i, i)] -= mt.Vector[i];
@@ -181,7 +181,7 @@ bool DSymmMatrix::Expand(DUppMatrix& s, DDiagMatrix& d) const
         return true;
 
     double sum, element;
-    register unsigned i, j, k;
+    unsigned i, j, k;
 
     for (j = 1; j < Dimension; ++j) {
         if (fabs(Buffer[j] * ZERO_ACCURACY) > fabs(d1))
@@ -333,7 +333,7 @@ bool DSymmMatrix::Inverse()
 
     l.Transpose(u);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; ++i) {
         vc_unit.Init(0.0);
@@ -370,7 +370,7 @@ DSymmMatrix& DSymmMatrix::RectSymmTRecmt(const DRectMem& A, const DSymmMatrix& B
 
     double element;
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
 
     for (i = 0; i < A.NRows; i++) {
         for (j = 0; j < B.Dimension; j++) {
@@ -407,7 +407,7 @@ DSymmMatrix& DSymmMatrix::RectSymmTRecmt(const CDSqrMatrix& A, const DSymmMatrix
 
     double element;
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
 
     for (i = 0; i < A.Dim(); i++) {
         for (j = 0; j < B.Dimension; j++) {
@@ -445,7 +445,7 @@ DSymmMatrix& DSymmMatrix::TRectSymmRecmt(const DRectMem& A, const DSymmMatrix& B
 
     double element;
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
 
     for (i = 0; i < A.NCols; i++) {
         for (j = 0; j < B.Dimension; j++) {
@@ -482,7 +482,7 @@ DSymmMatrix& DSymmMatrix::TRectSymmRecmt(const CDSqrMatrix& A, const DSymmMatrix
 
     double element;
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
 
     for (i = 0; i < A.Dim(); i++) {
         for (j = 0; j < B.Dimension; j++) {
@@ -514,7 +514,7 @@ DSymmMatrix& DSymmMatrix::RectTRecmt(const DRectMem& A)
 
     DTriagMem::Init(A.NRows);
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
 
     double sum;
 
@@ -539,7 +539,7 @@ void DSymmMatrix::GetLine(unsigned index, double* line) const
     TAssert(line != NULL);
     TAssert(index < Dimension);
 
-    register unsigned i;
+    unsigned i;
     for (i = 0; i < index; i++)
         line[i] = *(Buffer + Index(i, index));
     memcpy(line + index, Buffer + Index(index, index), (Dimension - index) * sizeof(double));
@@ -553,7 +553,7 @@ void DSymmMatrix::GetArray(double* array) const
     TAssert(Buffer != NULL);
     TAssert(array != NULL);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         memcpy(array, Buffer + Index(i, i), (Dimension - i) * sizeof(double));
@@ -577,7 +577,7 @@ double TRecmtSymmRecmt(const DVector& A, const DSymmMatrix& B)
 
     double element;
 
-    register unsigned j, k;
+    unsigned j, k;
 
     for (j = 0; j < B.Dimension; j++) {
         element = 0.0;

bnolib/DTriagMem.cpp

@@ -80,7 +80,7 @@ void DTriagMem::Init(double value)
 {
     TAssert(Buffer != NULL);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++)
         for (j = i; j < Dimension; j++)
@@ -94,7 +94,7 @@ void DTriagMem::Mult(double mltp)
 {
     TAssert(Buffer != NULL);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++)
         for (j = i; j < Dimension; j++)

bnolib/DUppMatrix.cpp

@@ -13,7 +13,7 @@ void DUppMatrix::Init(const DUppMatrix& mt)
     TAssert(mt.Buffer != NULL);
     DTriagMem::Init(mt.Dim());
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = i; j < Dimension; j++)
@@ -29,7 +29,7 @@ void DUppMatrix::Init(const DDiagMatrix& mt)
     DTriagMem::Init(mt.Dim());
 
     DTriagMem::Init(0.0);
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < Dimension; i++)
         Buffer[Index(i, i)] = mt.Vector[i];
@@ -44,7 +44,7 @@ DUppMatrix& DUppMatrix::operator+=(const DUppMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Buffer != NULL);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = i; j < Dimension; j++)
@@ -63,7 +63,7 @@ DUppMatrix& DUppMatrix::operator-=(const DUppMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Buffer != NULL);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = i; j < Dimension; j++)
@@ -81,7 +81,7 @@ DUppMatrix& DUppMatrix::operator+=(const DDiagMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Buffer != NULL);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < Dimension; i++)
         Buffer[Index(i, i)] += mt.Vector[i];
@@ -97,7 +97,7 @@ DUppMatrix& DUppMatrix::operator-=(const DDiagMatrix& mt)
     TAssert(mt.Dim() == Dimension);
     TAssert(Buffer != NULL);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < Dimension; i++)
         Buffer[Index(i, i)] -= mt.Vector[i];
@@ -118,7 +118,7 @@ DUppMatrix& DUppMatrix::MulRight(const DUppMatrix& mt)
     DUppMatrix src(Dimension);
     XChng(src);
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
     double* p;
 
     for (i = 0; i < Dimension; i++) {
@@ -146,7 +146,7 @@ DUppMatrix& DUppMatrix::MulLeft(const DUppMatrix& mt)
     DUppMatrix src(Dimension);
     XChng(src);
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
     double* p;
 
     for (i = 0; i < Dimension; i++) {
@@ -173,7 +173,7 @@ DUppMatrix& DUppMatrix::MulRight(const DDiagMatrix& mt)
     DUppMatrix src(Dimension);
     XChng(src);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = i; j < Dimension; j++)
@@ -195,7 +195,7 @@ DUppMatrix& DUppMatrix::MulLeft(const DDiagMatrix& mt)
     DUppMatrix src(Dimension);
     XChng(src);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++) {
         for (j = i; j < Dimension; j++)
@@ -214,7 +214,7 @@ DUppMatrix& DUppMatrix::Transpose(const DLowMatrix& mt)
 
     DTriagMem::Init(mt.Dim());
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < Dimension; i++)
         for (j = i; j < Dimension; j++)
@@ -246,7 +246,7 @@ void DUppMatrix::GetCol(unsigned index, double* col) const
     TAssert(col != NULL);
     TAssert(index < Dimension);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i <= index; i++)
         col[i] = *(Buffer + Index(i, index));
@@ -262,7 +262,7 @@ void DUppMatrix::GetArray(double* array) const
     TAssert(Buffer != NULL);
     TAssert(array != NULL);
 
-    register unsigned i;
+    unsigned i;
     for (i = 0; i < Dimension; i++) {
         GetLine(i, array);
         array += Dimension;

bnolib/DVector.cpp

@@ -28,7 +28,7 @@ double DVector::Abs(void) const
     TAssert(Vector != NULL);
 
     double s = 0;
-    for (register unsigned i = 0; i < Dimension; i++)
+    for (unsigned i = 0; i < Dimension; i++)
         s += Vector[i] * Vector[i];
 
     return sqrt(s);
@@ -46,7 +46,7 @@ DVector& DVector::MulLeft(const DUppMatrix& mt)
     ptr.Init(mt.Dimension);
     XChng(ptr);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < mt.Dimension; i++) {
         Vector[i] = 0;
@@ -69,7 +69,7 @@ DVector& DVector::MulLeft(const DLowMatrix& mt)
     ptr.Init(mt.Dimension);
     XChng(ptr);
 
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 0; i < mt.Dimension; i++) {
         Vector[i] = 0;
@@ -88,7 +88,7 @@ DVector& DVector::MulLeft(const DDiagMatrix& mt)
     TAssert(mt.Vector != NULL);
     TAssert(Dimension == mt.Dim());
 
-    register unsigned i;
+    unsigned i;
     for (i = 0; i < Dimension; i++)
         Vector[i] *= mt.Vector[i];
     return *this;
@@ -104,7 +104,7 @@ DVector& DVector::TRectRectVector(const CDMatrix& a, const CDMatrix& p, const DV
 
     Init(a.Cols());
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
 
     for (i = 0; i < Dimension; i++) {
         Vector[i] = 0.0;

bnolib/IMinMnk.cpp

@@ -554,7 +554,7 @@ double DMinMnk::VtMSV(const DVector& v, const DSymmMatrix& s)
 
     double element;
 
-    register unsigned j, k;
+    unsigned j, k;
 
     for (j = 0; j < n; j++) {
         element = 0.0;

bnolib/MatrFunc.cpp

@@ -19,7 +19,7 @@ DSqrMatrix& Mult(DSqrMatrix& dst, const DLowMatrix& src1, const DUppMatrix& src2
 
     dst.Init(src1.Dimension);
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
     double* line;
 
     for (i = 0; i < src1.Dimension; i++) {
@@ -46,7 +46,7 @@ DSqrMatrix& Mult(DSqrMatrix& dst, const DUppMatrix& src1, const DLowMatrix& src2
 
     dst.Init(src1.Dimension);
 
-    register unsigned i, j, k;
+    unsigned i, j, k;
     double* line;
 
     for (i = 0; i < src1.Dim(); i++) {
@@ -122,7 +122,7 @@ bool TDMSolution(const DVector& B, const DDiagMatrix& A, DVector& X)
 
     X.Init(A.Dimension);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < A.Dimension; i++) {
         if (fabs(A.Vector[i]) < ZERO_ACCURACY) {
@@ -154,8 +154,8 @@ bool TDMSolution(const DVector& b, const DUppMatrix& a, DVector& x)
     x.Vector[n - 1] = b.Vector[n - 1] / a.Buffer[a.Index(n - 1, n - 1)];
 
     double sum;
-    register int i;
-    register unsigned j;
+    int i;
+    unsigned j;
 
     for (i = n - 2; i >= 0; --i) {
         sum = 0.0;
@@ -190,7 +190,7 @@ bool TDMSolution(const DVector& b, const DLowMatrix& a, DVector& x)
     x.Vector[0] = b.Vector[0] / a.Buffer[0];
 
     double sum;
-    register unsigned i, j;
+    unsigned i, j;
 
     for (i = 1; i < n; i++) {
         sum = 0.0;
@@ -251,7 +251,7 @@ double operator*(const DVector& a, const DVector& b)
     TAssert(b.Vector != NULL);
 
     double s = 0.0;
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < a.Dimension; i++)
         s += a.Vector[i] * b.Vector[i];

bnolib/Matrix.cpp

@@ -19,7 +19,7 @@ void DMatrix::Transpose()
         return;
     }
 
-    register unsigned i, j;
+    unsigned i, j;
     double* line;
 
     for (i = 0; i < NCols; i++) {
@@ -45,7 +45,7 @@ void DMatrix::GetArray(double* array) const
     TAssert(Matrix != NULL);
     TAssert(array != NULL);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < NRows; i++) {
         memcpy(array, Matrix[i], NCols * sizeof(double));
@@ -62,7 +62,7 @@ void DMatrix::GetCol(unsigned index, double* col) const
     TAssert(index < NCols);
     TAssert(col != NULL);
 
-    register unsigned i;
+    unsigned i;
 
     for (i = 0; i < NRows; i++)
         col[i] = Matrix[i][index];

bnolib/ephemeris.cpp

@@ -66,7 +66,7 @@ bool Ephemeride::ReadData(const TimeJD& start, const TimeJD& end)
     //  BodyOffset - Смещение коэффициентов тела относительно начала 32-дневной порции
     //  CfLen - Число коэффициентов тела в одной порции
 
-    register int i;
+    int i;
     D32Len = 2;
     for (i = 0; i < (int)(sizeof(Order) / sizeof(int)); i++) {
         long cfl = 3 * (Order[i] + 1) * (32 / DInterval[i]);