topic MATMUL causing stack overflow in Intel® Fortran Compiler

MATMUL causing stack overflow

Roman1 — Mon, 11 Dec 2017 23:37:46 GMT

I noticed that for large matrices, MATMUL is crashing with a stack overflow. I can fix this with /heap-arrays0. The program does not crash when calling dgemm from MKL. I ran some tests, and the results from MATMUL and dgemm are identical. However MATMUL needs a large stack, and dgemm doesn't. Is this the correct behavior, or is there a bug somewhere?

Roman

Most of the time, the MATMUL

Lorri_M_Intel — Tue, 12 Dec 2017 13:31:10 GMT

Most of the time, the MATMUL intrinsic is implemented by in-line code. That is, the actual instructions to multiply each element of the matrices is generated by the compiler (including the loops to go from element to element, etc). That's different than a call to dgemm, where it is a simple routine call.

That said, because using /heap-arrays resolves the stack overflow, I'm going to guess that the compiler cannot detect that there won't be overlap between the result variable and the two operands, and for safety reasons it creates a temporary array. By default, this is on the stack.

To answer your question, I would characterize this as "expected behavior", and not likely a bug.

--Lorri

I ran a test, and the

Roman1 — Wed, 13 Dec 2017 00:03:07 GMT

I ran a test, and the performance of dgemm is slightly better than matmul. Based on your reply, this might be because there is an extra step where values are copied from the temporary memory to the result variable.

In the case where the matmul

TimP — Wed, 13 Dec 2017 10:08:08 GMT

In the case where the matmul result is stored explicitly A = matmul(B,C) Some compilers are able to avoid allocation of a temporary. Usually, the difference in performance would come mostly from efficiency of cache usage. In the more general case, A = b*matmul(C,D) + E It seems unlikely for a compiler to optimize the temporary away.