SHARE
TWEET

Untitled

a guest Aug 25th, 2019 65 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. void init_f_d_valentin(Model *m)
  2. {
  3.   int N = m->N;
  4.   int N_mat = m->N_mat;
  5.   int i;
  6.  
  7. //  m->d_mat = new crsMatrix * [N_mat];
  8. //  m->f_mat = new crsMatrix * [N_mat];
  9.  
  10. //  vector<pair<int , dcomplex> > * f_std;
  11. //  vector<pair<int , dcomplex> > * d_std;
  12.  
  13. //  Tensor_Coordinates *f_ijk = new Tensor_Coordinates;
  14. //  Tensor_Coordinates *d_ijk = new Tensor_Coordinates;
  15. //  Tensor_Coordinates_1 * f_1 = new Tensor_Coordinates_1[N_mat];
  16. //  Tensor_Coordinates_1 * d_1 = new Tensor_Coordinates_1[N_mat];
  17.  
  18. //  fijk_coord(f_ijk, N+1);
  19. //  sort_matrix(&f_ijk, f_1, N_mat);
  20. //  sort_matrix(f_ijk);
  21.  
  22. //  dijk_coord(d_ijk, N+1);
  23. //  sort_matrix(&d_ijk, d_1, N_mat);
  24. //  sort_matrix(d_ijk);
  25.  
  26. //  m->f_ijk = f_ijk;
  27. //  m->d_ijk = d_ijk;
  28.  
  29. //  int k, s, f;
  30. //
  31. //  k = 0;
  32. //  for( i = 0; i < N_mat; i++)
  33. //  {
  34. //    s = k;
  35. //    while (k < f_ijk->k)
  36. //    {
  37. //      if (f_ijk->coord1[k] != i) break;
  38. //      k++;
  39. //    }
  40. //    f = k;
  41. //    m->f_mat[i] = TensorToCrs(*f_ijk, s, f, N_mat);
  42. ////    m->f_mat[i] = TensorToCrs(f_1[i], 0, f_1[i].N, N_mat);
  43. //  }
  44. //
  45. //  k = 0;
  46. //  for( i = 0; i < N_mat; i++)
  47. //  {
  48. //    s = k;
  49. //    while (k < d_ijk->k)
  50. //    {
  51. //      if (d_ijk->coord1[k] != i) break;
  52. //      k++;
  53. //    }
  54. //    f = k;
  55. //    m->d_mat[i] = TensorToCrs(*d_ijk, s, f, N_mat);
  56. ////    m->d_mat[i] = TensorToCrs(d_1[i], 0, d_1[i].N, N_mat);
  57. //  }
  58.  
  59. //  for( i = 0; i < N_mat; i++)
  60. //  {
  61. //    delete[] f_1[i].coord2;
  62. //    delete[] f_1[i].coord3;
  63. //    delete[] f_1[i].data;
  64. //    delete[] d_1[i].coord2;
  65. //    delete[] d_1[i].coord3;
  66. //    delete[] d_1[i].data;
  67. //  }
  68. //  delete[] f_1;
  69. //  delete[] d_1;
  70.  
  71. //  free_matrix(&f_ijk);
  72. //  free_matrix(&d_ijk);
  73. }
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top