Oriol Colomés
Oriol Colomés
I went ahead and started the implementation of `AdjointFEOperator`. I pushed a first version in cdbebc173b7267860256ad9bfd148a49e8320793. You can find the new definition of an `AffineFEOperator` for cases 1 and 2...
> op_opd stores 2 matrices and 2 vectors. When selfadjoint=true the two matrices are the same julia object. Else, one needs to assemble a matrix with this new bilinear form...
For serial `AffineOperators` the numerical setup can also be reused: `lu(A') = transpose(lu(A))`, right?
Ok. No, it is not causing any error. I just thought it might be doing things twice because it's overwriting the same entries as in line 318, right? Or the...
In case that this line is fine, shouldn't it be returning `_mat` instead of `mat`?
Actually, I'm only looping for the first and if not found it goes to second and third. So, keeping aside the time of getting the closest 3 vertices instead of...
Hi @ConnorMallon , @amartinhuertas , @fverdugo, I'm getting the following error when trying to use autodiff for a problem with mixed-dimensional PDEs (a variable defined in a 2D domain and...
I'll prepare one reproducer and post it here
Hi @amartinhuertas , yes I forgot about this issue. Here I attach a reproducer: ```julia module tmp using Gridap domain = (0.0, 1.0, 0.0, 1.0) partition = (3,3) model_Ω =...
> To check if removing of `FESolver` will cause trouble in `GridapODEs` cc @oriolcg It should be fine. In `GridapODEs` we have `TransientFESolver` that receives an `ODESolver`. Each `ODESolver` calls...