ARRAY_OFFSET needed

This commit is contained in:
dwarning 2023-06-13 19:20:18 +02:00
parent c71171f6de
commit c1ab187e86
4 changed files with 6 additions and 3 deletions

View File

@ -1845,7 +1845,7 @@ register int I, Size = Matrix->Size;
register ElementPtr pElement; register ElementPtr pElement;
int NumberOfTies = 0; int NumberOfTies = 0;
long Product, MinMarkowitzProduct; long Product, MinMarkowitzProduct;
ElementPtr ChosenPivot, pLargestElement; ElementPtr ChosenPivot, pLargestElement = NULL;
RealNumber Magnitude, LargestElementMag, Ratio, RatioOfAccepted = 0, LargestInCol; RealNumber Magnitude, LargestElementMag, Ratio, RatioOfAccepted = 0, LargestInCol;
RealNumber FindLargestInCol(); RealNumber FindLargestInCol();

View File

@ -159,6 +159,8 @@ int *PrintOrdToIntRowMap, *PrintOrdToIntColMap;
SP_CALLOC( PrintOrdToIntColMap, int, Top + 1 ); SP_CALLOC( PrintOrdToIntColMap, int, Top + 1 );
if ( PrintOrdToIntRowMap == NULL OR PrintOrdToIntColMap == NULL) if ( PrintOrdToIntRowMap == NULL OR PrintOrdToIntColMap == NULL)
{ Matrix->Error = spNO_MEMORY; { Matrix->Error = spNO_MEMORY;
SP_FREE(PrintOrdToIntColMap);
SP_FREE(PrintOrdToIntRowMap);
return; return;
} }
for (I = 1; I <= Size; I++) for (I = 1; I <= Size; I++)

View File

@ -431,10 +431,10 @@ SMPcDProd(SMPmatrix *Matrix, SPcomplex *pMantissa, int *pExponent)
if (y < z) if (y < z)
y = z; y = z;
#ifdef debug_print
*pExponent = (int)(x + y); *pExponent = (int)(x + y);
x = scalbn(re, (int) -y); x = scalbn(re, (int) -y);
z = scalbn(im, (int) -y); z = scalbn(im, (int) -y);
#ifdef debug_print
printf(" ** values are: re %g, im %g, y %g, re' %g, im' %g\n", printf(" ** values are: re %g, im %g, y %g, re' %g, im' %g\n",
re, im, y, x, z); re, im, y, x, z);
#endif #endif

View File

@ -2192,7 +2192,8 @@ int Row, Col, Error;
} }
else ABORT(); else ABORT();
Matrix->Error = spOKAY; if (Matrix != NULL)
Matrix->Error = spOKAY;
return; return;
} }
#endif /* DOCUMENTATION */ #endif /* DOCUMENTATION */