void init_f_d_valentin(Model *m)
{
  int N = m->N;
  int N_mat = m->N_mat;
  int i;

//  m->d_mat = new crsMatrix * [N_mat];
//  m->f_mat = new crsMatrix * [N_mat];

//  vector<pair<int , dcomplex> > * f_std;
//  vector<pair<int , dcomplex> > * d_std;

//  Tensor_Coordinates *f_ijk = new Tensor_Coordinates;
//  Tensor_Coordinates *d_ijk = new Tensor_Coordinates;
//  Tensor_Coordinates_1 * f_1 = new Tensor_Coordinates_1[N_mat];
//  Tensor_Coordinates_1 * d_1 = new Tensor_Coordinates_1[N_mat];

//  fijk_coord(f_ijk, N+1);
//  sort_matrix(&f_ijk, f_1, N_mat);
//  sort_matrix(f_ijk);

//  dijk_coord(d_ijk, N+1);
//  sort_matrix(&d_ijk, d_1, N_mat);
//  sort_matrix(d_ijk);

//  m->f_ijk = f_ijk;
//  m->d_ijk = d_ijk;

//  int k, s, f;
//
//  k = 0;
//  for( i = 0; i < N_mat; i++)
//  {
//    s = k;
//    while (k < f_ijk->k)
//    {
//      if (f_ijk->coord1[k] != i) break;
//      k++;
//    }
//    f = k;
//    m->f_mat[i] = TensorToCrs(*f_ijk, s, f, N_mat);
////    m->f_mat[i] = TensorToCrs(f_1[i], 0, f_1[i].N, N_mat);
//  }
//
//  k = 0;
//  for( i = 0; i < N_mat; i++)
//  {
//    s = k;
//    while (k < d_ijk->k)
//    {
//      if (d_ijk->coord1[k] != i) break;
//      k++;
//    }
//    f = k;
//    m->d_mat[i] = TensorToCrs(*d_ijk, s, f, N_mat);
////    m->d_mat[i] = TensorToCrs(d_1[i], 0, d_1[i].N, N_mat);
//  }

//  for( i = 0; i < N_mat; i++)
//  {
//    delete[] f_1[i].coord2;
//    delete[] f_1[i].coord3;
//    delete[] f_1[i].data;
//    delete[] d_1[i].coord2;
//    delete[] d_1[i].coord3;
//    delete[] d_1[i].data;
//  }
//  delete[] f_1;
//  delete[] d_1;

//  free_matrix(&f_ijk);
//  free_matrix(&d_ijk);
}