Guest User

Untitled

a guest
Aug 25th, 2019
71
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