1//--------------------------------------------------------------------- 2//--------------------------------------------------------------------- 3#include "header.h" 4 5void add() { 6 7//--------------------------------------------------------------------- 8//--------------------------------------------------------------------- 9 10//--------------------------------------------------------------------- 11// addition of update to the vector u 12//--------------------------------------------------------------------- 13 14 int c, i, j, k, m; 15 16 for (c = 1; c <= ncells; c++) { 17 for (k = start(3,c); k <= cell_size(3,c)-end(3,c)-1; k++) { 18 for (j = start(2,c); j <= cell_size(2,c)-end(2,c)-1; j++) { 19 for (i = start(1,c); i <= cell_size(1,c)-end(1,c)-1; i++) { 20 for (m = 1; m <= 5; m++) { 21 u(m,i,j,k,c) = u(m,i,j,k,c) + rhs(m,i,j,k,c); 22 } 23 } 24 } 25 } 26 } 27 28 return; 29} 30