// -*- C++ -*-
#define N 10
#define K 1

generic main()
{
  stream ref_t A_rif;
  stream long[N][N] B;

  genera_A   (                    output_stream A_rif);
  copia_one  (input_stream A_rif  output_stream B);
  stampa        (input_stream B);
} 


genera_A (output_stream ref_t A_rif) {
  proc_genera_A (output_stream A_rif);
}


proc proc_genera_A (output_stream ref_t A_rif)
$c{ 
    int i,j;
  int displ;
  int dim = N*N*sizeof(long);  // dimensione area allocata con i reference 
  int size = N*sizeof(long);
  REF_Reference_t ref_matrice;
  long riga_A[N];
  
  // messaggio di partenza x il regression test
  printf("STARTING gen\n");

  /* allocazione dell'area in memoria condivisa */ 
  REF_New (&ref_matrice, dim, REF_Default, 0);
  
  for (i=0; i<N; i++) {
    /* inizializzazione della riga i-esima */
    for(j=0; j<N; j++)  
      riga_A[j] = j*i;

    displ = i*N*sizeof(long);
    
    /* scrittura in memoria condivisa della riga i-esima */
    REF_D_Write (&ref_matrice, riga_A, displ, size); 
  }
  assist_out (A_rif, &ref_matrice);
}c$



/* parmod per la generazione dello stream di righe */
parmod copia_one (input_stream ref_t A_rif  output_stream long B[N][N]) {

  topology one Pv;

  do input_section {
    guard1: on , , A_rif {
      distribution A_rif broadcast to Pv; 
    }
  } while (true)
 
  virtual_processors {
    elab_one (in guard1 out B) {
      VP {
	proc_elab_one (in A_rif output_stream B);
      }
    }
  }
 
  output_section {
    collects B from ANY Pv;
  }  
}


proc proc_elab_one (in ref_t A_rif output_stream long B[N][N])
inc<"iostream">
$c++{
  long B_[N][N];
  int displ;
  int size = N * sizeof(long);
  
  for(int i=0; i<N; i++) {
    /* lettura della riga i-esima dalla memoria condivisa */
    displ = i * N * sizeof(long);
    REF_D_Read (&(A_rif), &(B_[i][0]), displ, size); 
  }  
  assist_out (B, B_);
  REF_Delete (&A_rif); 
}c++$


stampa (input_stream long B[N][N])
inc<"iostream">
$c++{
  int ok = 0;

  // verifica del risultato 
  for (int i=0; i<N; i++) 
    for(int j=0; j<N; j++)  
      if (B[i][j] != j*i) ok = -1;
  
  if (ok==0)
    std::cerr << "ENDING with SUCCESS ref_copia.ast" << std::endl;
  else std::cerr << "ENDING with FAILURE ref_copia.ast" << std::endl;
}c++$