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