initdata.c

LPARX is a C++ class library for for implementing portable scientific applications on distributed me

#include "redblack2d.h"

//
// Initialize the boundary region with random data and
// the interior with zero
//

void initData(Grid2(double)& U)
  {
   for_2(j, k, U)
      U(j,k) = drand48();
   end_for

   Region2 interior = grow(U.region(), -1);

   for_2(j, k, interior);
      U(j,k) = 0.0;
   end_for
  }