order of projection = 2

Author: carlodev

src/Commons/Projection.jl

@@ -41,7 +41,7 @@ function create_coarse_spaces(params,simcase,order::Int64)
         @unpack  model = params
         simcase_coarse = deepcopy(simcase)
         # simcase_coarse.meshp.meshinfo.N .= ones(Int64, D) .* Int(ceil(N[1] / 2))
-        simcase_coarse.sprob.method.order = 1
+        simcase_coarse.sprob.method.order = 2
         boundary_conditions = create_boundary_conditions(simcase) 
         Vc, Uc, _, _ = creation_fe_spaces(simcase_coarse, model, boundary_conditions)
         merge!(params,Dict(:Uc=>Uc, :Vc=>Vc))
@@ -65,7 +65,7 @@ function project_solution(uh_fine, simcase::TaylorGreen{Periodic}, params::Dict{
     @sunpack D,N, order = simcase
     @unpack U,V,P,Q, Ω, degree,parts, dΩ = params
 
-    @assert order >1
+    @assert order >2
 
     Vc, Uc =   create_coarse_spaces(params,simcase,order)