bgeltz

#include <stdio.h>
#include <unistd.h>
#include <string.h>
#include <mpi.h>

int MPI_Init(int *argc, char **argv[])
{
    printf("In MPI_Init wrapper.\n");
    int err = PMPI_Init(argc, argv);
    printf("Finished with PMPI_Init. err = %d\n", err);

bgeltz

#include <stdio.h>
#include <unistd.h>
#include <string.h>
#include <mpi.h>

int main(int argc, char **argv) {
    MPI_Init(&argc, &argv);

    int s = 10;
    char hostname[100];

bgeltz

#include <stdio.h>
#include <unistd.h>
#include <string.h>
#include <mpi.h>

int main(int argc, char **argv) {
    /* MPI_Init(&argc, &argv); */

    int s = 1;
    char hostname[100];

bgeltz

$ srun -N1 -n8 ./hacc_tpm indat cmbM000.tf m000 INIT ALL_TO_ALL -w -R -N 512 -t 2x2x2

Partition 3D: [2:2:2]
Allocated heap: 5.45 MB (large) + 5.32 MB (small) used, 0.13 MB unused
Done loading parameters.
Initializing redistribution using a pencil layout on 8 ranks.
distribution 1D: [8:1:1]
distribution 3D: [2:2:2]
distribution 2z: [4:2:1]
distribution 2x: [1:4:2]

bgeltz

#!/bin/bash

export HACC_PLATFORM="mr-fusion"
export HACC_OBJDIR="${HACC_PLATFORM}"

ompflg="-DRCB_UNTHREADED_BUILD"

export HACC_CC=icc
export HACC_CXX=icpc
export HACC_MPI_CC=mpicc

bgeltz

1           # skip, stride in particle data output
0           # jinit, 0=>internal ICs, 1=>read ICs, 2=>restart ICs
1           # icbc, 0=>enforce periodic BCs, otherwise not (for jinit=1)
1           # norm, 0=>COBE 4-yr B/W, 1=> sigma8, 2=>COBE Peebles
0           # trans, 0=>T(k)=1, 1=>T(k)=H/S, 2=>T(k)=K/H, 3=>T(k)=CMBFAST
1           # initp, 0=>particle ICs before BCs, 1=>after, 2=>density ICc, else no print
2           # norder, order of time integration (locked to 2 here)
9           # frm, output format, 1=> ascii, 2=> binary, 3=> ascii, 2d density
2           # pbchoice, 1=>box units, 2=> physical units
1           # entest, 0=>Layzer-Irvine test, else not