OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Changed some ints to unsigned

Browse Source
This commit is contained in:
T.D. (Tongdong) Qiu 2020-06-25 12:29:20 +02:00
parent b692b66d3a
commit 4ad234bf55
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
opm/simulators/linalg/bda/BlockedMatrix.cpp
View File

@ -104,10 +104,10 @@ void sortBlockedRow(int *colIndices, double *data, int left, int right) {
template <unsigned int block_size> template <unsigned int block_size>
void blockMultSub(double *a, double *b, double *c) void blockMultSub(double *a, double *b, double *c)
{ {
for (int row = 0; row < block_size; row++) { for (unsigned int row = 0; row < block_size; row++) {
for (int col = 0; col < block_size; col++) { for (unsigned int col = 0; col < block_size; col++) {
double temp = 0.0; double temp = 0.0;
for (int k = 0; k < block_size; k++) { for (unsigned int k = 0; k < block_size; k++) {
temp += b[block_size * row + k] * c[block_size * k + col]; temp += b[block_size * row + k] * c[block_size * k + col];
} }
a[block_size * row + col] -= temp; a[block_size * row + col] -= temp;
@ -119,10 +119,10 @@ void blockMultSub(double *a, double *b, double *c)
template <unsigned int block_size> template <unsigned int block_size>
void blockMult(double *mat1, double *mat2, double *resMat) { void blockMult(double *mat1, double *mat2, double *resMat) {
for (int row = 0; row < block_size; row++) { for (unsigned int row = 0; row < block_size; row++) {
for (int col = 0; col < block_size; col++) { for (unsigned int col = 0; col < block_size; col++) {
double temp = 0; double temp = 0;
for (int k = 0; k < block_size; k++) { for (unsigned int k = 0; k < block_size; k++) {
temp += mat1[block_size * row + k] * mat2[block_size * k + col]; temp += mat1[block_size * row + k] * mat2[block_size * k + col];
} }
resMat[block_size * row + col] = temp; resMat[block_size * row + col] = temp;