parallelPartitioningType = ct.parallelPartitioningType nLayersOfOverlapForParallel = ct.nLayersOfOverlapForParallel restrictFineSolutionToAllMeshes = ct.restrictFineSolutionToAllMeshes triangleOptions = ct.triangleOptions timeIntegration = TimeIntegration.BackwardEuler_cfl stepController = StepControl.Min_dt_cfl_controller femSpaces = {0: ct.basis} elementQuadrature = ct.elementQuadrature elementBoundaryQuadrature = ct.elementBoundaryQuadrature massLumping = False numericalFluxType = Dissipation.NumericalFlux conservativeFlux = None subgridError = Dissipation.SubgridError(coefficients=physics.coefficients, nd=ct.domain.nd) shockCapturing = Dissipation.ShockCapturing( coefficients=physics.coefficients, nd=ct.domain.nd, shockCapturingFactor=ct.dissipation_shockCapturingFactor, lag=ct.dissipation_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:
from user_param import * timeIntegration = TimeIntegration.BackwardEuler_cfl stepController = StepControl.Min_dt_controller if timeIntegration == "VBDF": timeIntegration = TimeIntegration.VBDF if timeIntegration == "VBDF": timeOrder = 2 femSpaces = {0: basis} elementQuadrature = elementQuadrature elementBoundaryQuadrature = elementBoundaryQuadrature massLumping = False numericalFluxType = Dissipation.NumericalFlux conservativeFlux = None subgridError = Dissipation.SubgridError(coefficients=physics.coefficients, nd=nd) shockCapturing = Dissipation.ShockCapturing( coefficients=physics.coefficients, nd=nd, shockCapturingFactor=dissipation_shockCapturingFactor, lag=dissipation_lag_shockCapturing) fullNewtonFlag = True multilevelNonlinearSolver = NonlinearSolvers.Newton levelNonlinearSolver = NonlinearSolvers.Newton nonlinearSmoother = None linearSmoother = None matrix = LinearAlgebraTools.SparseMatrix multilevelLinearSolver = LinearSolvers.LU