/* heatSupportPrivate.c */

#include "heatSupport.h"
#include "heatSupportPrivate.h"

void updateInteriorTemperatures(double p[V][H]) {
  int i, j;
  double sourceElementTemperature;
  sourceElementTemperature = p[V/2][H/2];
  for (i = 1; i < V - 1; i = i + 1) {
    for (j = 1; j < H - 1; j = j + 1) {
      p[i][j] = (p[i-1][j] + p[i][j+1] + p[i+1][j] + p[i][j-1])/4;
    }
  }
  p[V/2][H/2] = sourceElementTemperature; /* restore source element */
}