[numerics] fix bug in extracting 5x5 diagonal block
@@ -1232,7 +1232,11 @@ void NM_extract_diag_block5(NumericsMatrix* M, int block_row_nb, double ** Block
{
case NM_DENSE:
{
- double* Mptr = M->matrix0 + (M->size0 + 1)*(block_row_nb + block_row_nb + block_row_nb);
+ double* Mptr = M->matrix0 + (M->size0 + 1)*(block_row_nb
+ + block_row_nb
+ + block_row_nb
+ + block_row_nb
+ + block_row_nb);
double* Bmat = *Block;
/* The part of MM which corresponds to the current block is copied into MLocal */
Bmat[0] = Mptr[0];
