Revert "Implementation of get_cellfield, get_derivative and get_derivatives functions for distributed transient cell fields"

This reverts commit 6ce000b.

Author: oriolcg

src/ODEs.jl

@@ -75,40 +75,13 @@ function ODEs.TransientCellField(f::DistributedCellField,derivatives::Tuple)
   DistributedCellField(fields,get_triangulation(f),f.metadata)
 end
 
-function ODEs.TransientCellField(f::DistributedCellField,derivatives::AbstractArray)
-  fields = map(local_views(f),local_views(derivatives)) do f, derivatives
-    ODEs.TransientCellField(f,derivatives)
-  end
-  DistributedCellField(fields,get_triangulation(f),f.metadata)
-end
-
 function ODEs.time_derivative(f::DistributedTransientSingleFieldCellField)
   fields = map(local_views(f)) do field
     ODEs.time_derivative(field)
   end
   DistributedCellField(fields,get_triangulation(f))
 end
 
-function ODEs.get_cellfield(f::DistributedTransientSingleFieldCellField)
-  cellfields = map(local_views(f)) do field
-    ODEs.get_cellfield(field)
-  end
-  DistributedCellField(cellfields,get_triangulation(f))
-end
-
-function ODEs.get_derivative(f::DistributedTransientSingleFieldCellField, k::Int)
-  derivatives = map(local_views(f)) do field
-    ODEs.get_derivative(field, k)
-  end
-  DistributedCellField(derivatives,get_triangulation(f))
-end
-
-function ODEs.get_derivatives(f::DistributedTransientSingleFieldCellField)
-  derivatives = map(local_views(f)) do field
-    ODEs.get_derivatives(field)
-  end
-end
-
 # MultiField FESpace
 
 function ODEs.has_transient(space::DistributedMultiFieldFESpace)
@@ -188,23 +161,3 @@ function ODEs.time_derivative(f::DistributedTransientMultiFieldCellField)
   part_fe_fun = map(ODEs.TransientMultiFieldCellField,fields)
   DistributedMultiFieldCellField(field_fe_fun,part_fe_fun)
 end
-
-function ODEs.get_cellfield(f::DistributedTransientMultiFieldCellField)
-  field_cellfield = map(ODEs.get_cellfield,f.field_fe_fun)
-  cellfields = to_parray_of_arrays(map(local_views,field_cellfield))
-  part_cellfields = map(MultiFieldCellField,cellfields)
-  DistributedMultiFieldCellField(field_cellfield,part_cellfields)
-end
-
-function ODEs.get_derivative(f::DistributedTransientMultiFieldCellField, k::Int)
-  field_derivative = map(df -> ODEs.get_derivative(df, k), f.field_fe_fun)
-  derivatives = to_parray_of_arrays(map(local_views,field_derivative))
-  part_derivatives = map(MultiFieldCellField,derivatives)
-  DistributedMultiFieldCellField(field_derivative,part_derivatives)
-end
-
-function ODEs.get_derivatives(f::DistributedTransientMultiFieldCellField)
-  derivatives = map(local_views(f)) do field
-    ODEs.get_derivatives(field)
-  end
-end

test/AutodiffTests.jl

@@ -6,7 +6,6 @@ using GridapDistributed
 using PartitionedArrays
 using SparseArrays
 using ForwardDiff
-using Gridap.ODEs: TransientCellField, get_jacs
 
 function main_sf(distribute,parts)
   ranks = distribute(LinearIndices((prod(parts),)))
@@ -207,54 +206,11 @@ function skeleton_mf_different_fespace_trians(distribute,parts)
   end
 end
 
-function main_transient_sf(distribute,parts)
-  ranks = distribute(LinearIndices((prod(parts),)))
-
-  domain = (0,4,0,4)
-  cells = (4,4)
-  model = CartesianDiscreteModel(ranks,parts,domain,cells)
-
-  u((x,y),t) = (x+y)^k + 2*t
-  u(t::Real) = x -> u(x,t)
-  σ(∇u) = (1.0+∇u⋅∇u)*∇u
-  dσ(∇du,∇u) = (2*∇u⋅∇du)*∇u + (1.0+∇u⋅∇u)*∇du
-  f(t) = x -> ∂t(u)(t)(x) - divergence(y->σ(∇(u(t),y)),x)
-
-  k = 1
-  reffe = ReferenceFE(lagrangian,Float64,k)
-  V = TestFESpace(model,reffe,dirichlet_tags="boundary")
-  U = TransientTrialFESpace(V,u)
-
-  Ω = Triangulation(model)
-  dΩ = Measure(Ω,2*k)
-  r(t,u,v) = ∫( ∂t(u)⋅v + ∇(v)⋅(σ∘∇(u)) - v*f(t) )dΩ
-  j_0(t,u,du,v) = ∫( ∇(v)⋅(dσ∘(∇(du),∇(u))) )dΩ
-  j_t(t,u,dut,v) = ∫( dut⋅v )dΩ
-
-  op = TransientFEOperator(r,(j_0,j_t),U,V)
-  op_AD = TransientFEOperator(r,U,V)
-
-  uh = interpolate(0.0,U(0.0))
-  ∂tuₕ = interpolate(0.0,U(0.0))
-  uhₜ = TransientCellField(uh,(∂tuₕ,))
-  
-  b = assemble_vector(v->r(1.0,uhₜ,v),V)
-  A_0 = assemble_matrix((du,v)->j_0(1.0,uhₜ,du,v),U(1.0),V)
-  A_t = assemble_matrix((dut,v)->j_t(1.0,uhₜ,dut,v),U(1.0),V)
-  jac_0_AD = get_jacs(op_AD)[1]
-  jac_t_AD = get_jacs(op_AD)[2]
-  A_0_AD = assemble_matrix((du,v)->jac_0_AD(1.0,uhₜ,du,v),U(1.0),V)
-  A_t_AD = assemble_matrix((dut,v)->jac_t_AD(1.0,uhₜ,dut,v),U(1.0),V)
-  @test reduce(&,map(≈,partition(A_0),partition(A_0_AD)))
-  @test reduce(&,map(≈,partition(A_t),partition(A_t_AD)))
-end
-
 function main(distribute,parts)
   main_sf(distribute,parts)
   main_mf(distribute,parts)
   mf_different_fespace_trians(distribute,parts)
   skeleton_mf_different_fespace_trians(distribute,parts)
-  main_transient_sf(distribute,parts)
 end
 
 end

test/sequential/runtests.jl

@@ -46,8 +46,4 @@ end
   include("AdaptivityTests.jl")
 end
 
-@time @testset "AutodiffTests" begin
-  include("AutodiffTests.jl")
-end
-
 end # module