From: Ramachandra Kuravi on
1) I am trying to solve a composite beam response problem. The size of GLOBAL STIFFNESS MATRIX that I need should be of the order od 80thousand*80 thousand. I evaluated local stiffness matrices for each beam separately and assemgled them to GLOBAL MATRIX.

I could initialize GLOBAL MATRIX to the desired size using spalloc. I could obtain local stiffness matrices individually. But, when I try to assemble Local stiffness matrices into Global {one by one}, I am told that "memory has run out".
I need such high grid size to study convergence issues.

Schematic of the code:
Main function:
1)Initializes Global matrix and calls functions that evaluate Local stiffness matrices and assembles them. {Initializing done with SPALLOC}.

Subsidiary functions:
1)They, when called return local stiffness matrices to the Main function.