timeOrder = 2 else: timeIntegration = TimeIntegration.BackwardEuler_cfl stepController = StepControl.Min_dt_controller femSpaces = {0: ct.basis} elementQuadrature = ct.elementQuadrature elementBoundaryQuadrature = ct.elementBoundaryQuadrature massLumping = False numericalFluxType = Kappa.NumericalFlux conservativeFlux = None subgridError = Kappa.SubgridError(coefficients=physics.coefficients, nd=nd) shockCapturing = Kappa.ShockCapturing(coefficients=physics.coefficients, nd=nd, shockCapturingFactor=ct.kappa_shockCapturingFactor, lag=ct.kappa_lag_shockCapturing) fullNewtonFlag = True multilevelNonlinearSolver = NonlinearSolvers.Newton levelNonlinearSolver = NonlinearSolvers.Newton nonlinearSmoother = None linearSmoother = None #printNonlinearSolverInfo = True matrix = LinearAlgebraTools.SparseMatrix if not ct.useOldPETSc and not ct.useSuperlu: multilevelLinearSolver = LinearSolvers.KSP_petsc4py levelLinearSolver = LinearSolvers.KSP_petsc4py
timeIntegration = TimeIntegration.BackwardEuler_cfl stepController = StepControl.Min_dt_controller if timeIntegration == "VBDF": timeIntegration = TimeIntegration.VBDF if timeIntegration == "VBDF": timeOrder = 2 femSpaces = {0: basis} massLumping = False numericalFluxType = Kappa.NumericalFlux conservativeFlux = None subgridError = Kappa.SubgridError(coefficients=physics.coefficients, nd=nd) shockCapturing = Kappa.ShockCapturing( coefficients=physics.coefficients, nd=nd, #shockCapturingFactor=user_param.kappa_shockCapturingFactor, #lag=user_param.kappa_lag_shockCapturing) shockCapturingFactor=kappa_shockCapturingFactor, lag=kappa_lag_shockCapturing) fullNewtonFlag = True multilevelNonlinearSolver = NonlinearSolvers.Newton levelNonlinearSolver = NonlinearSolvers.Newton nonlinearSmoother = None linearSmoother = None #printNonlinearSolverInfo = True matrix = LinearAlgebraTools.SparseMatrix if not useOldPETSc and not useSuperlu: