1. A method of increasing at least one of efficiency and speed in executing a matrix subroutine on a computer, said method comprising: for a cache on said computer having a working area size defined by dimensions k and m, k and m being integers greater than 1, with k>m, determining an integer a and, for data of a matrix to be processed by said matrix subroutine, a square block of matrix data having a size (a*m)×(a*m), such that k≧(a*m); and preliminarily storing data contiguously for a matrix subroutine call in a computer memory in an increment block size that is based on said cache working area size, a first dimension of said block size being (a*m) and a second dimension of said block size being m, thereby a of said increment blocks constitute said square block of matrix data of size (a*m)×(a*m), and each said increment block fits into said cache working area having said dimensions defined by k and m.