Hi,
I dont get what is going on here. One would think I would get an off diag.
let mut XI_matrix = DMatrix::from_iterator(xi_rows,xi_cols,xi_weights.iter().cloned());
println!("XI_matrix: {}",XI_matrix);
xi_weights: [1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0]
XI_matrix:
┌ ┐
│ 1 1 1 1 1 1 1 1 │
│ 0 0 0 0 0 0 0 0 │
│ 0 0 0 0 0 0 0 0 │
│ 0 0 0 0 0 0 0 0 │
│ 1 1 1 1 1 1 1 1 │
│ 0 0 0 0 0 0 0 0 │
│ 0 0 0 0 0 0 0 0 │
│ 0 0 0 0 0 0 0 0 │
└ ┘