Skip to content
This repository has been archived by the owner on Feb 28, 2022. It is now read-only.

Commit

Permalink
Updated StokesAutodiff
Browse files Browse the repository at this point in the history
  • Loading branch information
oriolcg committed Oct 5, 2021
1 parent 2caea13 commit c866ee8
Showing 1 changed file with 6 additions and 6 deletions.
12 changes: 6 additions & 6 deletions test/TransientFEsTests/StokesEquationAutoDiffTests.jl
Original file line number Diff line number Diff line change
@@ -1,12 +1,12 @@
module StokesEquationTests
module StokesEquationAutoDiffTests

using Gridap
using ForwardDiff
using LinearAlgebra
using Test
using GridapODEs.ODETools
using GridapODEs.TransientFETools
using Gridap.FESpaces: get_algebraic_operator
using Gridap.FESpaces
using Gridap.Arrays: test_array

# using GridapODEs.ODETools: ThetaMethodLinear
Expand Down Expand Up @@ -76,11 +76,11 @@ dy = get_fe_basis(Y)
dx = get_trial_fe_basis(X₀)
xh = FEFunction(X₀,rand(num_free_dofs(X₀)))

cell_j = get_array(jac(0.5,xh,xh,dx,dy))
cell_j_t = get_array(jac_t(0.5,xh,xh,dx,dy))
cell_j = get_array(jac(0.5,(xh,xh),dx,dy))
cell_j_t = get_array(jac_t(0.5,(xh,xh),dx,dy))

cell_j_auto = get_array(jacobian(x->res(0.5,x,xh,dy),xh))
cell_j_t_auto = get_array(jacobian(x->res(0.5,xh,x,dy),xh))
cell_j_auto = get_array(jacobian(x->res(0.5,(x,xh),dy),xh))
cell_j_t_auto = get_array(jacobian(x->res(0.5,(xh,x),dy),xh))

test_array(cell_j_auto,cell_j,)
test_array(cell_j_t_auto,cell_j_t,)
Expand Down

0 comments on commit c866ee8

Please sign in to comment.