<div dir="ltr"><div>Hi Barry,</div><div><br></div><div>Thanks for popping in. This is my code:</div><div><br></div><div><i>                blockValues = this%getDiagonalBlockValues(iElement=iElement)<br>                call MatSetValuesBlocked(this%A, 4-bdim, mesh%cellGlobalAddr(iElement)-1, 4-bdim, mesh%cellGlobalAddr(iElement)-1, blockValues, INSERT_VALUES, ierr)</i><br></div><div><br></div><div><i>blockValues </i>is a 3x3 or 4x4 matrix and I am passing it straight away.<br></div></div>