Here are the two source-code files:

spongeserrmd.cpp - Includes the RMD algorithm

The input parameters are the following:

//
// Input for the program RMD part
//

int N = 8;                   // The size of the cube in 
power of 2
double Hurst_exponent = 0.1; // The input of the exponent for the sponge
double start_value = 1.;     // The inital values inside 
the grid

The output is saved in dma_data_*size*_*H*.dat where

*size* : is the size of the grid
*H* : is the value of H

spongeserdma.cpp - Includes the DMA algorithm

The input parameters are the following:

//
// Input for the program DMA part
//

int N = 8;			// The size of the cube in power of 2
double Hurst_exponent = 0.1;	// The input of the exponent 
for the sponge
double theta = 0.5;             // The values of theta in 
order to calculate the average
int n_max = 23;                 // Max values in the sum 
n_max <<N_max
int n_min = 3;                  // Starting value for the 
sum >= 3
int m = 0;                      // Default values for the 
m for the average case to use

The output saved in dma_data_*size*_*H*.dat where

*size* : is the size of the grid
*H* : is the value of H

is needed.