eschnett/cactusamrex - Change JNXOHF6UIF4SD6K7KESWOLMHGAIN6GUEFVRFFI6LJP6EUYFIL47QC

TestArrayGroup: Added tests for arrays and DynamicData The thorn provides a parameter file for these tests. Arrays (distrib=const) are written in one function and read in another function. The DynamicData test calls the DynamicData function for a grid function, array, and scalar and compares the data to the correct values. These tests also serve as a basic test of the validity mechanisms. The initialize function writes all three types of variables everywhere. Then, the DynamicData function reads all of them everywhere.

Created by Samuel Cupp on October 9, 2020

JNXOHF6UIF4SD6K7KESWOLMHGAIN6GUEFVRFFI6LJP6EUYFIL47QC

In channels

main

Change contents

File addition: TestArrayGroup (drwxr-xr-x)
[1.0]

File addition: configuration.ccl (-rw-r--r--)

[0.17]

# Configuration definitions for thorn TestArrayGroup
REQUIRES CarpetX
PROVIDES TestArrayGroup
{
}

File addition: interface.ccl (-rw-r--r--)

[0.17]

# Interface definition for thorn CarpetX
IMPLEMENTS: TestArrayGroup
CCTK_REAL test_array[4] TYPE=array DIM=2 SIZE=5,6 DISTRIB=constant
{
  test1 test2 test3
} "Test arrays for verifying CarpetX implementation of non-distributed array data"
CCTK_REAL test_gf TYPE=gf "Test grid function for verifying CarpetX implementation of DynamicData"
CCTK_REAL test_scalar TYPE=scalar "Test scalar for verifying CarpetX implementation of DynamicData"

File addition: par (drwxr-xr-x)
[0.17]

File addition: test.par (-rw-r--r--)

[0.607]

ActiveThorns = "
    CarpetX
    IOUtil
    TestArrayGroup
"
 
$nlevels = 1
$ncells = 32
Cactus::cctk_show_schedule = no
Cactus::presync_mode = "mixed-error"
Cactus::cctk_itlast = 2
CarpetX::ncells_x = $ncells
CarpetX::ncells_y = $ncells
CarpetX::ncells_z = $ncells
CarpetX::max_num_levels = $nlevels
CarpetX::dtfac = 0.5
IO::out_dir = "arraytest"
IO::out_every = 1

File addition: param.ccl (-rw-r--r--)

[0.17]

# Parameter definitions for thorn TestArrayGroup

File addition: schedule.ccl (-rw-r--r--)

[0.17]

# Schedule definitions for thorn CarpetX
SCHEDULE TestArrayGroup_Initialize AT Initial
{
  LANG: C
  WRITES: test_scalar(everywhere)
  WRITES: test_array(everywhere)
  WRITES: test_gf(everywhere)
} "Initialize data in the test arrays"
SCHEDULE TestArrayGroup_Compare AT PostInitial
{
  LANG: C
  READS: test_array(everywhere)
} "Test data in the test arrays"
SCHEDULE TestArrayGroup_DynamicData AT PostInitial
{
  LANG: C
  READS: test_scalar(everywhere)
  READS: test_array(everywhere)
  READS: test_gf(everywhere)
} "Test DynamicData for grid functions, scalars, and distrib=const arrays"

File addition: src (drwxr-xr-x)
[0.17]

File addition: TestArray.cxx (-rw-r--r--)

[0.1673]

#include <cctk.h>
#include <cctk_Arguments.h>
#include <cctk_Parameters.h>
#include <util_Table.h>
#include <cassert>
#include <cmath>
#include <memory>
#include <iostream>
#include <string>
#include <utility>
#include <algorithm>
#include <array>
#include <cstdint>
#include <functional>
#include <tuple>
#include <type_traits>
#include <vector>
extern "C" void TestArrayGroup_Initialize(cGH *cctkGH) {
  DECLARE_CCTK_PARAMETERS;
  DECLARE_CCTK_ARGUMENTS_TestArrayGroup_Initialize;
  int imax = 5;
  int jmax = 6;
  int nmax = 4;
  // Initialize test array
  for(int n=0; n<nmax; n++) {
    for(int j=0; j<jmax; j++) {
      for(int i=0; i<imax; i++) {
        int index = i + j*imax + n*imax*jmax;
        test1[index] = 1+i*j*n;
        test2[index] = 1+7*i*j*n;
        test3[index] = 1+13*i*j*n;
      }
    }
  }
  // Initialize test grid function
  int gi = CCTK_GroupIndex("TestArrayGroup::test_gf");
  cGroup group;
  int ierr = CCTK_GroupData(gi, &group);
  assert(!ierr);
  std::vector<int> lsh(3);
  ierr = CCTK_GrouplshGI(cctkGH, group.dim, lsh.data(), gi);
  assert(!ierr);
 for(int k=0; k < lsh[2]; k++) {
   for(int j=0; j < lsh[1]; j++) {
     for(int i=0; i < lsh[0]; i++) {
       int index = CCTK_GFINDEX3D(cctkGH, i, j, k);
       test_gf[index] = i*j*k;
     }
   }
 }
  *test_scalar = 1;
 return;
}
extern "C" void TestArrayGroup_Compare(cGH *cctkGH) {
  DECLARE_CCTK_PARAMETERS;
  DECLARE_CCTK_ARGUMENTS_TestArrayGroup_Compare;
  int imax = 5;
  int jmax = 6;
  int nmax = 4;
  int error_count[3] = {0,0,0};
  for(int n=0; n<nmax; n++) {
    for(int j=0; j<jmax; j++) {
      for(int i=0; i<imax; i++) {
        int index = i + j*imax + n*imax*jmax;
//        std::string vname = "TestArrayGroup::test1[" + std::to_string(index) + "]";
//        CCTK_REAL *var1 = (CCTK_REAL *)CCTK_VarDataPtr(cctkGH,0,vname.c_str());
//
//        vname = "TestArrayGroup::test2[" + std::to_string(index) + "]";
//        CCTK_REAL *var2 = (CCTK_REAL *)CCTK_VarDataPtr(cctkGH,0,vname.c_str());
//
//        vname = "TestArrayGroup::test3[" + std::to_string(index) + "]";
//        CCTK_REAL *var3 = (CCTK_REAL *)CCTK_VarDataPtr(cctkGH,0,vname.c_str());
        if (test1[index] != 1+i*j*n) error_count[0] += 1;
        if (test2[index] != 1+7*i*j*n) error_count[1] += 1;
        if (test3[index] != 1+13*i*j*n) error_count[2] += 1;
        //if (*var1 != 1+i*j*n) error_count[0] += 1;
        //if (*var2 != 1+7*i*j*n) error_count[1] += 1;
        //if (*var3 != 1+13*i*j*n) error_count[2] += 1;
      }
    }
  }
  if (error_count[0] > 0) {
    int size = nmax*jmax*imax;
    std::string msg;
    msg = "TestArrayGroup: grid array test1 failed in " + std::to_string(error_count[0]) + 
          " of " + std::to_string(size) + " elements.";
    CCTK_ERROR(msg.c_str());
  }
  if (error_count[1] > 0) {
    int size = nmax*jmax*imax;
    std::string msg;
    msg = "TestArrayGroup: grid array test2 failed in " + std::to_string(error_count[1]) + 
          " of " + std::to_string(size) + " elements.";
    CCTK_ERROR(msg.c_str());
  }
  if (error_count[2] > 0) {
    int size = nmax*jmax*imax;
    std::string msg;
    msg = "TestArrayGroup: grid array test3 failed in " + std::to_string(error_count[2]) + 
          " of " + std::to_string(size) + " elements.";
    CCTK_ERROR(msg.c_str());
  }
 return;
}

File addition: TestDynamicData.cxx (-rw-r--r--)

[0.1673]

#include <cctk.h>
#include <cctk_Arguments.h>
#include <cctk_Parameters.h>
#include <util_Table.h>
#include <cassert>
#include <cmath>
#include <memory>
#include <iostream>
#include <string>
#include <utility>
#include <algorithm>
#include <array>
#include <cstdint>
#include <functional>
#include <tuple>
#include <type_traits>
#include <vector>
#include <map>
extern "C" void TestArrayGroup_DynamicData(cGH *cctkGH) {
  DECLARE_CCTK_PARAMETERS;
  DECLARE_CCTK_ARGUMENTS;
  std::map<std::string, int> error = {
                {"lsh", 0}, {"ash", 0}, {"gsh", 0},
                {"lbnd", 0}, {"ubnd", 0}, {"bbox", 0},
                {"nghostzones", 0} };
  CCTK_INT gf_gi, array_gi, scalar_gi;
  cGroupDynamicData gf_data, array_data, scalar_data;
  gf_gi = CCTK_GroupIndex("TestArrayGroup::test_gf");
  array_gi = CCTK_GroupIndex("TestArrayGroup::test_array");
  scalar_gi = CCTK_GroupIndex("TestArrayGroup::test_scalar");
  CCTK_GroupDynamicData(cctkGH, gf_gi, &gf_data);
  CCTK_GroupDynamicData(cctkGH, array_gi, &array_data);
  CCTK_GroupDynamicData(cctkGH, scalar_gi, &scalar_data);
  // Validate grid function dynamic data
  cGroup gf_group;
  int ierr = CCTK_GroupData(gf_gi, &gf_group);
  if (ierr) CCTK_ERROR("TestArrayGroup: error in GroupData for grid functions");
  if (gf_data.dim != gf_group.dim || gf_data.dim <= 0) CCTK_ERROR(
            "TestArrayGroup: incorrect dimension in grid function dynamic data");
  if (gf_data.activetimelevels != CCTK_ActiveTimeLevelsGI(cctkGH, gf_gi)) CCTK_ERROR(
            "TestArrayGroup: incorrect activetimelevels in grid function dynamic data");
  for (int i = 0; i < gf_data.dim; i++) {
    if (gf_data.lsh[i] != cctkGH->cctk_lsh[i]) error.at("lsh") += 1;
    if (gf_data.ash[i] != cctkGH->cctk_ash[i]) error.at("ash") += 1;
    if (gf_data.gsh[i] != cctkGH->cctk_gsh[i]) error.at("gsh") += 1;
    if (gf_data.lbnd[i] != cctkGH->cctk_lbnd[i]) error.at("lbnd") += 1;
    if (gf_data.ubnd[i] != cctkGH->cctk_ubnd[i]) error.at("lbnd") += 1;
    if (gf_data.bbox[2*i] != cctkGH->cctk_bbox[2*i]) error.at("bbox") += 1;
    if (gf_data.bbox[2*i+1] != cctkGH->cctk_bbox[2*i+1]) error.at("bbox") += 1;
    if (gf_data.nghostzones[i] != cctkGH->cctk_nghostzones[i]) error.at("nghostzones") += 1;
  }
  if (error.at("lsh")) CCTK_ERROR(
            "TestArrayGroup: incorrect lsh data in grid function dynamic data");
  if (error.at("ash")) CCTK_ERROR(
            "TestArrayGroup: incorrect ash data in grid function dynamic data");
  if (error.at("gsh")) CCTK_ERROR(
            "TestArrayGroup: incorrect gsh data in grid function dynamic data");
  if (error.at("lbnd")) CCTK_ERROR(
            "TestArrayGroup: incorrect lbnd data in grid function dynamic data");
  if (error.at("ubnd")) CCTK_ERROR(
            "TestArrayGroup: incorrect ubnd data in grid function dynamic data");
  if (error.at("bbox")) CCTK_ERROR(
            "TestArrayGroup: incorrect bbox data in grid function dynamic data");
  if (error.at("nghostzones")) CCTK_ERROR(
            "TestArrayGroup: incorrect nghostzones data in grid function dynamic data");
  // Validate grid scalar dynamic data
  for ( auto & [key, val] : error) {val = 0;}
  if (scalar_data.dim != 0) CCTK_ERROR(
            "TestArrayGroup: incorrect dimension in scalar dynamic data");
  if (scalar_data.activetimelevels != 1) CCTK_ERROR(
            "TestArrayGroup: incorrect activetimelevels in scalar dynamic data");
  for (int i = 0; i < gf_data.dim; i++) {
    if (scalar_data.lsh[i] != -1) error.at("lsh") += 1;
    if (scalar_data.ash[i] != -1) error.at("ash") += 1;
    if (scalar_data.gsh[i] != -1) error.at("gsh") += 1;
    if (scalar_data.lbnd[i] != -1) error.at("lbnd") += 1;
    if (scalar_data.ubnd[i] != -1) error.at("lbnd") += 1;
    if (scalar_data.bbox[2*i] != -1) error.at("bbox") += 1;
    if (scalar_data.bbox[2*i+1] != -1) error.at("bbox") += 1;
    if (scalar_data.nghostzones[i] != -1) error.at("nghostzones") += 1;
  }
  if (error.at("lsh")) CCTK_ERROR(
            "TestArrayGroup: incorrect lsh data in scalar dynamic data");
  if (error.at("ash")) CCTK_ERROR(
            "TestArrayGroup: incorrect ash data in scalar dynamic data");
  if (error.at("gsh")) CCTK_ERROR(
            "TestArrayGroup: incorrect gsh data in scalar dynamic data");
  if (error.at("lbnd")) CCTK_ERROR(
            "TestArrayGroup: incorrect lbnd data in scalar dynamic data");
  if (error.at("ubnd")) CCTK_ERROR(
            "TestArrayGroup: incorrect ubnd data in scalar dynamic data");
  if (error.at("bbox")) CCTK_ERROR(
            "TestArrayGroup: incorrect bbox data in scalar dynamic data");
  if (error.at("nghostzones")) CCTK_ERROR(
            "TestArrayGroup: incorrect nghostzones data in scalar dynamic data");
  // Validate grid array dynamic data
  for ( auto & [key, val] : error) {val = 0;}
  cGroup array_group;
  ierr = CCTK_GroupData(array_gi, &array_group);
  if (ierr) CCTK_ERROR("TestArrayGroup: error in GroupData for arrays");
  int sz[2] = {5,6};
  if (array_data.dim != array_group.dim || array_data.dim <= 0) CCTK_ERROR(
            "TestArrayGroup: incorrect dimension in array dynamic data");
  if (scalar_data.activetimelevels != 1) CCTK_ERROR(
            "TestArrayGroup: incorrect activetimelevels in array dynamic data");
  for (int i = 0; i < array_data.dim; i++) {
    if (array_data.lsh[i] != sz[i]) error.at("lsh") += 1;
    if (array_data.ash[i] != sz[i]) error.at("ash") += 1;
    if (array_data.gsh[i] != sz[i]) error.at("gsh") += 1;
    if (array_data.lbnd[i] != 0) error.at("lbnd") += 1;
    if (array_data.ubnd[i] != sz[i]-1) error.at("lbnd") += 1;
    if (array_data.bbox[2*i] != 1) error.at("bbox") += 1;
    if (array_data.bbox[2*i+1] != 1) error.at("bbox") += 1;
    if (array_data.nghostzones[i] != 0) error.at("nghostzones") += 1;
  }
  if (error.at("lsh")) CCTK_ERROR(
            "TestArrayGroup: incorrect lsh data in array dynamic data");
  if (error.at("ash")) CCTK_ERROR(
            "TestArrayGroup: incorrect ash data in array dynamic data");
  if (error.at("gsh")) CCTK_ERROR(
            "TestArrayGroup: incorrect gsh data in array dynamic data");
  if (error.at("lbnd")) CCTK_ERROR(
            "TestArrayGroup: incorrect lbnd data in array dynamic data");
  if (error.at("ubnd")) CCTK_ERROR(
            "TestArrayGroup: incorrect ubnd data in array dynamic data");
  if (error.at("bbox")) CCTK_ERROR(
            "TestArrayGroup: incorrect bbox data in array dynamic data");
  if (error.at("nghostzones")) CCTK_ERROR(
            "TestArrayGroup: incorrect nghostzones data in array dynamic data");
 return;
}

File addition: make.code.defn (-rw-r--r--)

[0.1673]

# Main make.code.defn file for thorn CarpetX
# Source files in this directory
SRCS =						\
	TestArray.cxx \
	TestDynamicData.cxx
# Subdirectories containing source files
SUBDIRS =

In channels

Change contents

File addition: TestArrayGroup (drwxr-xr-x)

File addition: configuration.ccl (-rw-r--r--)

File addition: interface.ccl (-rw-r--r--)

File addition: par (drwxr-xr-x)

File addition: test.par (-rw-r--r--)

File addition: param.ccl (-rw-r--r--)

File addition: schedule.ccl (-rw-r--r--)

File addition: src (drwxr-xr-x)

File addition: TestArray.cxx (-rw-r--r--)

File addition: TestDynamicData.cxx (-rw-r--r--)

File addition: make.code.defn (-rw-r--r--)