Advertisement
Guest User

Untitled

a guest
Jun 26th, 2013
53
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.64 KB | None | 0 0
  1. index_t * C_dims = SG_MALLOC(index_t, 3);
  2. C_dims[0] = 10;
  3. C_dims[1] = 10;
  4. C_dims[2] = 30;
  5. SGNDArray< float64_t > C(C_dims, 3);
  6.  
  7. // Load Data (is transposed)
  8. SGMatrix<float64_t> C_flat(C_dims[1]*C_dims[2], C_dims[0]);
  9. CAsciiFile C_in("C.csv", 'r');
  10. C_flat.load(&C_in);
  11. C_in.close();
  12.  
  13. // Reshape Data
  14. Eigen::Map<EMatrix> EC_flat(C_flat.matrix, C_dims[1]*C_dims[2], C_dims[0]); // (still transposed)
  15. for(int i = 0; i < C_dims[2]; i++)
  16. {
  17. Eigen::Map<EMatrix> Ci(C.get_matrix(i), C_dims[0], C_dims[1]);
  18. EMatrix src = EC_flat.block(i*C_dims[1], 0, C_dims[1], C_dims[0]);
  19. Ci = src.transpose(); // un-transpose here
  20. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement