This is supposed to be a mutating function, so I assume you meant du .= alpha*D*Q*u. Otherwise you’re creating a new vector called du instead of writing into the data of the buffer.
1 Like
This is supposed to be a mutating function, so I assume you meant du .= alpha*D*Q*u. Otherwise you’re creating a new vector called du instead of writing into the data of the buffer.