@@ -2845,34 +2845,6 @@ end
28452845 # the sparsity detection will never depend on the hyperbolic part of the problem
28462846 @test jac_prototype_parabolic == jac_prototype_hyperbolic_parabolic
28472847
2848- # ##############################################################################
2849- # ## Compare sparsity pattern based on hyperbolic-parabolic ###
2850- # ## problem with sparsity pattern of parabolic-only problem ###
2851-
2852- # We construct a semidiscretization just as we did previously, but this time with advection_velocity=0
2853- equations_hyperbolic_zero_advection = LinearScalarAdvectionEquation1D (0 )
2854- equations_parabolic_zero_advection = LaplaceDiffusion1D (diffusivity (),
2855- equations_hyperbolic_zero_advection)
2856-
2857- semi_jac_type_zero_advection = SemidiscretizationHyperbolicParabolic (mesh,
2858- (equations_hyperbolic_zero_advection,
2859- equations_parabolic_zero_advection),
2860- initial_condition_convergence_test,
2861- solver,
2862- uEltype = jac_eltype)
2863-
2864- # Do sparsity detection on our semidiscretization with advection turned off
2865- rhs_hyp_para_wrapped! = (du_ode, u0_ode) -> rhs_hyperbolic_parabolic! (du_ode, u0_ode,
2866- semi_jac_type_zero_advection,
2867- tspan[1 ])
2868-
2869- jac_prototype_parabolic_zero_advection = jacobian_sparsity (rhs_hyp_para_wrapped!,
2870- du_ode, u0_ode,
2871- jac_detector)
2872-
2873- # Given that the stencil for parabolic solvers are always larger than those of hyperbolic solvers,
2874- # the sparsity detection will never depend on the hyperbolic part of the problem
2875- @test jac_prototype_parabolic == jac_prototype_parabolic_zero_advection
28762848end
28772849end
28782850
0 commit comments