Hi,<br><br>I'm a sophomore student at lanl, relatively new to fortran. I'm supposed find a way to square large sparse matrices in parallel. I found a way to do this with the dbcsr part of cp2k, I was wondering how you would go about setting the distribution of a dbcsr_matrix based on a square matrix you read in from a file (for example: REAL(8):: DATA(10, 10) ). <br><br>Any help is greatly appreciated,<br>Adam<br>