[markb@marsh ~/Gale]$ ./Gale-1_5_0 work/inputPol_Okaya.xml Finding module: "Underworld"... Finding module: "PICellerator"... Finding module: "StgFEM"... Finding module: "StgDomain"... Creating Stg_Components from the component-list Instantiating UnderworldContext as context Instantiating FeMesh as mesh-linear Instantiating CartesianGenerator as linearMesh-generator Instantiating MeshVariable as velocity Instantiating CompositeVC as velocityBCs Instantiating CompositeVC as velocityICs Instantiating DofLayout as velocityDofLayout Instantiating FeVariable as VelocityField Instantiating OperatorFeVariable as VelocityMagnitudeField Instantiating OperatorFeVariable as VelocityGradientsField Instantiating OperatorFeVariable as VelocityGradientsInvariantField Instantiating OperatorFeVariable as StrainRateField Instantiating OperatorFeVariable as VorticityField Instantiating OperatorFeVariable as StrainRateInvariantField Instantiating MeshVariable as pressure Instantiating DofLayout as pressureDofLayout Instantiating FeVariable as PressureField Instantiating StressField as StressField Instantiating ViscosityField as ViscosityField Instantiating SingleCellLayout as cellLayout Instantiating GaussParticleLayout as particleLayout Instantiating IntegrationPointsSwarm as gaussSwarm Instantiating GaussMapper as gaussMapper Instantiating BackgroundParticleLayout as backgroundLayout Instantiating ParticleMovementHandler as gaussMSwarmMovementHandler Instantiating MaterialPointsSwarm as gaussMaterialSwarm Instantiating TimeIntegrator as timeIntegrator Instantiating ElementCellLayout as elementCellLayout Instantiating PCDVC as weights Instantiating MappedParticleLayout as localLayout Instantiating IntegrationPointsSwarm as picIntegrationPoints Instantiating CoincidentMapper as mapper Instantiating MeshParticleLayout as materialSwarmParticleLayout Instantiating ParticleMovementHandler as pMovementHandler Instantiating ParticleShadowSync as pShadowSync Instantiating MaterialPointsSwarm as materialSwarm Instantiating SwarmAdvector as materialSwarmAdvector Instantiating SolutionVector as solutionVelocity Instantiating SolutionVector as solutionPressure Instantiating ForceVector as mom_force Instantiating ForceVector as cont_force Instantiating StiffnessMatrix as k_matrix Instantiating ConstitutiveMatrixCartesian as constitutiveMatrix Instantiating StiffnessMatrix as g_matrix Instantiating GradientStiffnessMatrixTerm as gradientStiffnessMatrixTerm Instantiating StiffnessMatrix as preconditioner Instantiating UzawaPreconditionerTerm as preconditionerTerm Instantiating Stokes_SLE_UzawaSolver as uzawa Instantiating Stokes_SLE as stokesEqn Instantiating StiffnessMatrix as c_matrix Instantiating MixedStabiliserTerm as mixedStabiliser Instantiating BuoyancyForceTerm as buoyancyForceTerm Instantiating Everywhere as background Instantiating EscapedRoutine as escapedRoutine Instantiating RegularRemesherCmpt as velocityRemesher Instantiating Box as boxShape Instantiating MaterialViscosity as backgroundViscosity Instantiating StoreVisc as storeViscosity Instantiating StoreStress as storeStress Instantiating RheologyMaterial as viscous Constructing Stg_Components from the live-component register Constructing context.. Constructing mesh-linear.. Stg_Component 'mesh-linear' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting parameter 'isElementMesh': Using default value = false Getting parameter 'elementType': Found - Value = "linear" Constructing linearMesh-generator.. Looking for essential component 'context': Found. Stg_Component 'linearMesh-generator' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for non-essential component 'mesh-linear': Found. Getting parameter 'dim': Found - Value = "dim" Key 'dim' points to key 'dim' in the root dictionary: Found - Value = "3" Getting parameter 'dims': Using default value = 3 Getting parameter 'maxDecomposedDims': Using default value = 0 Getting parameter 'contactDepth-left': Using default value = 0 Getting parameter 'contactDepth-right': Using default value = 0 Getting parameter 'contactDepth-bottom': Using default value = 0 Getting parameter 'contactDepth-top': Using default value = 0 Getting parameter 'contactDepth-back': Using default value = 0 Getting parameter 'contactDepth-front': Using default value = 0 Getting parameter 'contactGeometry-x': Using default value = 0 Getting parameter 'contactGeometry-y': Using default value = 0 Getting parameter 'contactGeometry-z': Using default value = 0 Getting unsigned int from root dictionary with key 'restartTimestep' and default value '0' Getting parameter 'shadowDepth': Found - Value = "shadowDepth" Key 'shadowDepth' points to key 'shadowDepth' in the root dictionary: Found - Value = "1" Getting parameter 'regular': Found - Value = "False" Getting parameter 'periodic_x': Using default value = false Getting parameter 'periodic_y': Using default value = false Getting parameter 'periodic_z': Using default value = false Constructing velocity.. Stg_Component 'velocity' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'mesh-linear': Found. Getting parameter 'topologicalDim': Using default value = 0 Getting parameter 'VectorComponentCount': Found - Value = "dim" Key 'VectorComponentCount' points to key 'dim' in the root dictionary: Found - Value = "3" Constructing velocityBCs.. Stg_Component 'velocityBCs' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting parameter 'vcName': Using default value = "velocityBCs" Looking for non-essential component 'mesh-linear': Found. Constructing velocityICs.. Stg_Component 'velocityICs' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting parameter 'vcName': Using default value = "velocityICs" Looking for non-essential component 'mesh-linear': Found. Constructing velocityDofLayout.. Stg_Component 'velocityDofLayout' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'mesh-linear': Found. Getting parameter 'BaseVariableCount': Found - Value = "dim" Key 'BaseVariableCount' points to key 'dim' in the root dictionary: Found - Value = "3" Looking for Variable 'vx' in Variable_Register Looking for Variable 'vy' in Variable_Register Looking for Variable 'vz' in Variable_Register Constructing VelocityField.. Stg_Component 'VelocityField' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'fieldComponentCount': Using default value = 0 Looking for essential component 'mesh-linear': Found. Stg_Component 'VelocityField' cannot find non-essential component with key 'GeometryMesh'. Looking for essential component 'velocityDofLayout': Found. Looking for non-essential component 'velocityICs': Found. Looking for non-essential component 'velocityBCs': Found. Looking for non-essential component 'velocityLinkedDofs': Not found. Getting parameter 'isReferenceSolution': Using default value = false Getting parameter 'loadReferenceEachTimestep': Using default value = false Getting parameter 'removeBCs': Using default value = true Constructing VelocityMagnitudeField.. Stg_Component 'VelocityMagnitudeField' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'fieldComponentCount': Using default value = 0 Getting parameter 'Operator': Found - Value = "Magnitude" Looking for FeVariable 'VelocityField' in fieldVariable_Register. Constructing VelocityGradientsField.. Stg_Component 'VelocityGradientsField' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'fieldComponentCount': Using default value = 0 Getting parameter 'Operator': Found - Value = "Gradient" Looking for FeVariable 'VelocityField' in fieldVariable_Register. Constructing VelocityGradientsInvariantField.. Stg_Component 'VelocityGradientsInvariantField' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'fieldComponentCount': Using default value = 0 Getting parameter 'Operator': Found - Value = "TensorInvariant" Looking for FeVariable 'VelocityGradientsField' in fieldVariable_Register. Constructing StrainRateField.. Stg_Component 'StrainRateField' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'fieldComponentCount': Using default value = 0 Getting parameter 'Operator': Found - Value = "TensorSymmetricPart" Looking for FeVariable 'VelocityGradientsField' in fieldVariable_Register. Constructing VorticityField.. Stg_Component 'VorticityField' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'fieldComponentCount': Using default value = 0 Getting parameter 'Operator': Found - Value = "TensorAntisymmetricPart" Looking for FeVariable 'VelocityGradientsField' in fieldVariable_Register. Constructing StrainRateInvariantField.. Stg_Component 'StrainRateInvariantField' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'fieldComponentCount': Using default value = 0 Getting parameter 'Operator': Found - Value = "SymmetricTensor_Invariant" Looking for FeVariable 'StrainRateField' in fieldVariable_Register. Constructing pressure.. Stg_Component 'pressure' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'mesh-linear': Found. Getting parameter 'topologicalDim': Using default value = 0 Constructing pressureDofLayout.. Stg_Component 'pressureDofLayout' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'mesh-linear': Found. Getting parameter 'BaseVariableCount': Using default value = 1 Looking for Variable 'pressure' in Variable_Register Constructing PressureField.. Stg_Component 'PressureField' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'fieldComponentCount': Using default value = 0 Looking for essential component 'mesh-linear': Found. Stg_Component 'PressureField' cannot find non-essential component with key 'GeometryMesh'. Looking for essential component 'pressureDofLayout': Found. Stg_Component 'PressureField' cannot find non-essential component with key 'IC'. Stg_Component 'PressureField' cannot find non-essential component with key 'BC'. Looking for non-essential component 'pressureLinkedDofs': Not found. Getting parameter 'isReferenceSolution': Using default value = false Getting parameter 'loadReferenceEachTimestep': Using default value = false Getting parameter 'removeBCs': Using default value = true Constructing StressField.. Looking for non-essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'fieldComponentCount': Using default value = 0 Looking for essential component 'picIntegrationPoints': Found. picIntegrationPoints has not been constructed yet. Constructing now. Constructing picIntegrationPoints.. Stg_Component 'picIntegrationPoints' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'ParticleType': Using default value = "IntegrationPoint" Looking for essential component 'elementCellLayout': Found. elementCellLayout has not been constructed yet. Constructing now. Constructing elementCellLayout.. Stg_Component 'elementCellLayout' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'mesh-linear': Found. Looking for essential component 'localLayout': Found. localLayout has not been constructed yet. Constructing now. Constructing localLayout.. Stg_Component 'localLayout' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting parameter 'cellParticleTblDelta': Using default value = 4 Getting parameter 'extraParticlesFactor': Using default value = 0.05 Looking for list 'ParticleCommHandlers': Not Found. Stg_Component 'picIntegrationPoints' cannot find non-essential component with key 'IC'. Looking for essential component 'mesh-linear': Found. Looking for essential component 'timeIntegrator': Found. timeIntegrator has not been constructed yet. Constructing now. Constructing timeIntegrator.. Getting parameter 'order': Found - Value = "1" Getting parameter 'simultaneous': Found - Value = "False" Looking for non-essential component 'context': Found. Looking for non-essential component 'weights': Found. weights has not been constructed yet. Constructing now. Constructing weights.. Stg_Component 'weights' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'resolution': Using default value = 10 Getting parameter 'resolutionX': Found - Value = "10" Getting parameter 'resolutionY': Found - Value = "10" Getting parameter 'resolutionZ': Found - Value = "10" Looking for essential component 'materialSwarm': Found. materialSwarm has not been constructed yet. Constructing now. Constructing materialSwarm.. Stg_Component 'materialSwarm' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'ParticleType': Using default value = "IntegrationPoint" Looking for essential component 'elementCellLayout': Found. Looking for essential component 'materialSwarmParticleLayout': Found. materialSwarmParticleLayout has not been constructed yet. Constructing now. Constructing materialSwarmParticleLayout.. Stg_Component 'materialSwarmParticleLayout' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting parameter 'cellParticleCount': Found - Value = "particlesPerCell" Key 'cellParticleCount' points to key 'particlesPerCell' in the root dictionary: Found - Value = "30" Getting parameter 'seed': Using default value = 13 Looking for essential component 'mesh-linear': Found. Getting parameter 'cellParticleTblDelta': Using default value = 4 Getting parameter 'extraParticlesFactor': Using default value = 0.05 Looking for list 'ParticleCommHandlers': Found. Looking for non-essential component 'pMovementHandler': Found. pMovementHandler has not been constructed yet. Constructing now. Constructing pMovementHandler.. Looking for non-essential component 'pShadowSync': Found. pShadowSync has not been constructed yet. Constructing now. Constructing pShadowSync.. Stg_Component 'materialSwarm' cannot find non-essential component with key 'IC'. Looking for essential component 'mesh-linear': Found. Looking for non-essential component 'escapedRoutine': Found. escapedRoutine has not been constructed yet. Constructing now. Constructing escapedRoutine.. Stg_Component 'escapedRoutine' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'particlesToRemoveDelta': Using default value = 20 Stg_Component 'materialSwarm' cannot find non-essential component with key 'Material'. Getting parameter 'upperT': Found - Value = "25" Getting parameter 'lowerT': Found - Value = "0.6" Getting parameter 'maxDeletions': Found - Value = "3" Getting parameter 'maxSplits': Found - Value = "3" Getting parameter 'splitInInterfaceCells': Using default value = false Getting parameter 'deleteInInterfaceCells': Using default value = false Getting parameter 'Inflow': Found - Value = "True" Getting parameter 'Threshold': Using default value = 0.8 Getting parameter 'CentPosRatio': Using default value = 0.01 Getting unsigned int from root dictionary with key 'particlesPerCell' and default value '25' Looking for essential component 'mapper': Found. mapper has not been constructed yet. Constructing now. Constructing mapper.. Stg_Component 'mapper' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'picIntegrationPoints': Found. Looking for essential component 'materialSwarm': Found. Getting parameter 'recalculateWeights': Using default value = true Looking for non-essential component 'context': Found. Looking for essential component 'mesh-linear': Found. Looking for essential component 'StrainRateField': Found. Looking for essential component 'constitutiveMatrix': Found. constitutiveMatrix has not been constructed yet. Constructing now. Constructing constitutiveMatrix.. Stg_Component 'constitutiveMatrix' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'k_matrix': Found. k_matrix has not been constructed yet. Constructing now. Constructing k_matrix.. Stg_Component 'k_matrix' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'VelocityField': Found. Looking for essential component 'VelocityField': Found. Looking for non-essential component 'mom_force': Found. mom_force has not been constructed yet. Constructing now. Constructing mom_force.. Stg_Component 'mom_force' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'VelocityField': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Stg_Component 'k_matrix' cannot find non-essential component with key 'ApplicationDepInfo'. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'isNonLinear': Using default value = false Getting parameter 'allowZeroElementContributions': Found - Value = "False" Stg_Component 'k_matrix' cannot find non-essential component with key 'transposeRHS'. Looking for essential component 'picIntegrationPoints': Found. Stg_Component 'constitutiveMatrix' cannot find non-essential component with key 'ExtraInfo'. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'storeConstitutiveMatrix': Using default value = false Getting parameter 'StressVariable': Using default value = "Stress" Stg_Component 'StressField' cannot find non-essential component with key 'SLE'. Constructing ViscosityField.. Looking for non-essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'fieldComponentCount': Using default value = 0 Looking for essential component 'picIntegrationPoints': Found. Looking for non-essential component 'context': Found. Looking for essential component 'mesh-linear': Found. Looking for essential component 'constitutiveMatrix': Found. Constructing cellLayout.. Stg_Component 'cellLayout' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'dimExistsI': Using default value = true Getting parameter 'dimExistsJ': Using default value = true Getting parameter 'dimExistsK': Using default value = true Getting parameter 'minX': Using default value = -1 Getting parameter 'minY': Using default value = -1 Getting parameter 'minZ': Using default value = -1 Getting parameter 'maxX': Using default value = 1 Getting parameter 'maxY': Using default value = 1 Getting parameter 'maxZ': Using default value = 1 Constructing particleLayout.. Stg_Component 'particleLayout' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'gaussParticles': Using default value = 2 Getting parameter 'gaussParticlesX': Using default value = 2 Getting parameter 'gaussParticlesY': Using default value = 2 Getting parameter 'gaussParticlesZ': Using default value = 2 Constructing gaussSwarm.. Stg_Component 'gaussSwarm' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'ParticleType': Using default value = "IntegrationPoint" Looking for essential component 'cellLayout': Found. Looking for essential component 'particleLayout': Found. Getting parameter 'cellParticleTblDelta': Using default value = 4 Getting parameter 'extraParticlesFactor': Using default value = 0.05 Looking for list 'ParticleCommHandlers': Not Found. Stg_Component 'gaussSwarm' cannot find non-essential component with key 'IC'. Looking for essential component 'mesh-linear': Found. Looking for essential component 'timeIntegrator': Found. Stg_Component 'gaussSwarm' cannot find non-essential component with key 'WeightsCalculator'. Looking for essential component 'gaussMapper': Found. gaussMapper has not been constructed yet. Constructing now. Constructing gaussMapper.. Stg_Component 'gaussMapper' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'gaussSwarm': Found. Looking for essential component 'gaussMaterialSwarm': Found. gaussMaterialSwarm has not been constructed yet. Constructing now. Constructing gaussMaterialSwarm.. Stg_Component 'gaussMaterialSwarm' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'ParticleType': Using default value = "IntegrationPoint" Looking for essential component 'cellLayout': Found. Looking for essential component 'backgroundLayout': Found. backgroundLayout has not been constructed yet. Constructing now. Constructing backgroundLayout.. Stg_Component 'backgroundLayout' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting parameter 'cellParticleTblDelta': Using default value = 4 Getting parameter 'extraParticlesFactor': Using default value = 0.05 Looking for list 'ParticleCommHandlers': Found. Looking for non-essential component 'gaussMSwarmMovementHandler': Found. gaussMSwarmMovementHandler has not been constructed yet. Constructing now. Constructing gaussMSwarmMovementHandler.. Stg_Component 'gaussMaterialSwarm' cannot find non-essential component with key 'IC'. Looking for essential component 'mesh-linear': Found. Stg_Component 'gaussMaterialSwarm' cannot find non-essential component with key 'EscapedRoutine'. Stg_Component 'gaussMaterialSwarm' cannot find non-essential component with key 'Material'. Getting parameter 'recalculateWeights': Using default value = true Constructing materialSwarmAdvector.. Stg_Component 'materialSwarmAdvector' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Stg_Component 'materialSwarmAdvector' cannot find non-essential component with key 'Variable'. Looking for essential component 'timeIntegrator': Found. Looking for list 'data': Not Found. Getting parameter 'allowFallbackToFirstOrder': Found - Value = "True" Looking for essential component 'VelocityField': Found. Looking for essential component 'materialSwarm': Found. Looking for non-essential component 'periodicBCsManager': Not found. Constructing solutionVelocity.. Stg_Component 'solutionVelocity' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'VelocityField': Found. Constructing solutionPressure.. Stg_Component 'solutionPressure' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'PressureField': Found. Constructing cont_force.. Stg_Component 'cont_force' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'PressureField': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Constructing g_matrix.. Stg_Component 'g_matrix' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'VelocityField': Found. Looking for essential component 'PressureField': Found. Looking for non-essential component 'mom_force': Found. Stg_Component 'g_matrix' cannot find non-essential component with key 'ApplicationDepInfo'. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'isNonLinear': Using default value = false Getting parameter 'allowZeroElementContributions': Found - Value = "False" Looking for non-essential component 'cont_force': Found. Constructing gradientStiffnessMatrixTerm.. Stg_Component 'gradientStiffnessMatrixTerm' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'g_matrix': Found. Looking for essential component 'gaussSwarm': Found. Stg_Component 'gradientStiffnessMatrixTerm' cannot find non-essential component with key 'ExtraInfo'. Constructing preconditioner.. Stg_Component 'preconditioner' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'PressureField': Found. Looking for essential component 'PressureField': Found. Looking for non-essential component 'cont_force': Found. Stg_Component 'preconditioner' cannot find non-essential component with key 'ApplicationDepInfo'. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'isNonLinear': Using default value = false Getting parameter 'allowZeroElementContributions': Found - Value = "True" Stg_Component 'preconditioner' cannot find non-essential component with key 'transposeRHS'. Constructing preconditionerTerm.. Stg_Component 'preconditionerTerm' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'preconditioner': Found. Looking for essential component 'picIntegrationPoints': Found. Stg_Component 'preconditionerTerm' cannot find non-essential component with key 'ExtraInfo'. Constructing uzawa.. Stg_Component 'uzawa' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting parameter 'statSolve': Using default value = false Getting parameter 'statReps': Using default value = 0 Getting parameter 'tolerance': Found - Value = "linearTolerance" Key 'tolerance' points to key 'linearTolerance' in the root dictionary: Using default value = 1e-05 Getting parameter 'maxIterations': Found - Value = "5000" Getting parameter 'minIterations': Found - Value = "1" Getting parameter 'useAbsoluteTolerance': Using default value = false Getting parameter 'monitor': Found - Value = "false" Looking for non-essential component 'preconditioner': Found. Constructing stokesEqn.. Looking for non-essential component 'uzawa': Found. Getting parameter 'makeConvergenceFile': Found - Value = "false" Getting parameter 'isNonLinear': Found - Value = "true" Getting parameter 'nonLinearTolerance': Found - Value = "nonLinearTolerance" Key 'nonLinearTolerance' points to key 'nonLinearTolerance' in the root dictionary: Using default value = 0.01 Getting parameter 'nonLinearMaxIterations': Found - Value = "nonLinearMaxIterations" Key 'nonLinearMaxIterations' points to key 'nonLinearMaxIterations' in the root dictionary: Using default value = 500 Getting parameter 'killNonConvergent': Using default value = true Getting parameter 'nonLinearMinIterations': Found - Value = "nonLinearMinIterations" Key 'nonLinearMinIterations' points to key 'nonLinearMinIterations' in the root dictionary: Using default value = 1 Getting parameter 'nonLinearSolutionType': Using default value = "default" Getting parameter 'optionsPrefix': Using default value = "" Getting parameter 'picard_FormFunctionType': Using default value = "PicardFormFunction_KSPResidual" Getting parameter 'picard_alpha': Using default value = 1 Getting parameter 'picard_rtol': Using default value = 1e-08 Getting parameter 'picard_atol': Using default value = 1e-50 Getting parameter 'picard_xtol': Using default value = 1e-08 Getting parameter 'picard_ActivateMonitor': Using default value = false Looking for non-essential component 'context': Found. Getting parameter 'linearSolveInitialGuess': Using default value = false Looking for essential component 'k_matrix': Found. Looking for essential component 'g_matrix': Found. Looking for non-essential component '': Not found. Looking for non-essential component 'c_matrix': Found. c_matrix has not been constructed yet. Constructing now. Constructing c_matrix.. Stg_Component 'c_matrix' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'PressureField': Found. Looking for essential component 'PressureField': Found. Looking for non-essential component 'cont_force': Found. Stg_Component 'c_matrix' cannot find non-essential component with key 'ApplicationDepInfo'. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'isNonLinear': Using default value = false Getting parameter 'allowZeroElementContributions': Found - Value = "True" Stg_Component 'c_matrix' cannot find non-essential component with key 'transposeRHS'. Looking for essential component 'solutionVelocity': Found. Looking for essential component 'solutionPressure': Found. Looking for essential component 'mom_force': Found. Looking for essential component 'cont_force': Found. Constructing mixedStabiliser.. Stg_Component 'mixedStabiliser' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'c_matrix': Found. Looking for essential component 'gaussSwarm': Found. Stg_Component 'mixedStabiliser' cannot find non-essential component with key 'ExtraInfo'. Looking for essential component 'picIntegrationPoints': Found. Looking for essential component 'storeViscosity': Found. storeViscosity has not been constructed yet. Constructing now. Constructing storeViscosity.. Stg_Component 'storeViscosity' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'materialSwarm': Found. Constructing buoyancyForceTerm.. Stg_Component 'buoyancyForceTerm' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'mom_force': Found. Looking for essential component 'picIntegrationPoints': Found. Stg_Component 'buoyancyForceTerm' cannot find non-essential component with key 'ExtraInfo'. Stg_Component 'buoyancyForceTerm' cannot find non-essential component with key 'TemperatureField'. Getting parameter 'gravity': Found - Value = "1" Getting parameter 'adjust': Using default value = false Stg_Component 'buoyancyForceTerm' cannot find non-essential component with key 'HydrostaticTerm'. Constructing background.. Stg_Component 'background' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'CentreX': Using default value = 0 Getting parameter 'CentreY': Using default value = 0 Getting parameter 'CentreZ': Using default value = 0 Getting parameter 'invert': Using default value = false Getting parameter 'alpha': Using default value = 0 Getting parameter 'beta': Using default value = 0 Getting parameter 'gamma': Using default value = 0 Constructing velocityRemesher.. Looking for essential component 'mesh-linear': Found. Getting parameter 'contactDepth': Using default value = 0 Getting parameter 'contactSize': Using default value = 0 Getting parameter 'diffuseCorners': Using default value = false Getting parameter 'diffuseSurface': Using default value = false Getting parameter 'diffusionCoef': Using default value = 1 Looking for essential component 'context': Found. Constructing boxShape.. Stg_Component 'boxShape' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting unsigned int from root dictionary with key 'dim' and default value '0' Getting parameter 'CentreX': Using default value = 0 Getting parameter 'CentreY': Using default value = 0 Getting parameter 'CentreZ': Using default value = 0 Getting parameter 'invert': Using default value = false Getting parameter 'alpha': Using default value = 0 Getting parameter 'beta': Using default value = 0 Getting parameter 'gamma': Using default value = 0 Getting parameter 'startX': Found - Value = "minX" Key 'startX' points to key 'minX' in the root dictionary: Found - Value = "225.0f" Getting parameter 'endX': Found - Value = "maxX" Key 'endX' points to key 'maxX' in the root dictionary: Found - Value = "900.0f" Getting parameter 'startY': Found - Value = "minY" Key 'startY' points to key 'minY' in the root dictionary: Found - Value = "0.0f" Getting parameter 'endY': Found - Value = "maxY" Key 'endY' points to key 'maxY' in the root dictionary: Found - Value = "750.0f" Getting parameter 'startZ': Found - Value = "minZ" Key 'startZ' points to key 'minZ' in the root dictionary: Found - Value = "-400.0f" Getting parameter 'endZ': Found - Value = "maxZ" Key 'endZ' points to key 'maxZ' in the root dictionary: Found - Value = "-100.0f" Constructing backgroundViscosity.. Stg_Component 'backgroundViscosity' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Getting parameter 'eta0': Found - Value = "1.0" Constructing storeStress.. Stg_Component 'storeStress' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'materialSwarm': Found. First Trying to find component by name 'StrainRateField': Found. Constructing viscous.. Stg_Component 'viscous' cannot find non-essential component with key 'Context'. Looking for essential component 'context': Found. Looking for essential component 'boxShape': Found. Looking for list 'Rheology': Found. Looking for essential component 'backgroundViscosity': Found. Looking for essential component 'storeViscosity': Found. Looking for essential component 'storeStress': Found. Stg_Component 'viscous' cannot find non-essential component with key 'Compressible'. Getting parameter 'isCompressible': Using default value = false Constructing EulerDeform.. Looking for essential component 'context': Found. Looking for essential component 'timeIntegrator': Found. Looking for non-essential component 'ArtDisplacementField': Not found. Constructing Underworld_VTKOutput.. Looking for essential component 'context': Found. Constructing StgFEM_StandardConditionFunctions.. Looking for essential component 'context': Found. Building Stg_Components from the live-component register Cartesian generator: 'linearMesh-generator' Target mesh: 'mesh-linear' Global element size: 90x100x40 Local offset of rank 0: 0x0x0 Local range of rank 0: 90x100x40 Generating vertices... ... done. Generating elements... ... done. Generating edges... ... done. Generating faces... ... done. Generating element-vertex incidence... ... done. Generating face-edge incidence... ... done. Generating edge-vertex incidence... ... done. Generating face-vertex incidence... ... done. Generating vertex neighbours... ... done. Generating geometry... ... done. Generating element types... ... element types are 'Mesh_HexType', ... mesh algorithm type is 'Mesh_Algorithms', ... done. Assigning FeMesh element types... ... FE element types are 'TrilinearElementType', ... done. Stiffness matrix: 'k_matrix' Calculating number of nonzero entries... Found 84022490 nonzero entries. Done. Stiffness matrix: 'g_matrix' Calculating number of nonzero entries... Found 28563927 nonzero entries. Done. Stiffness matrix: 'c_matrix' Calculating number of nonzero entries... Found 9870091 nonzero entries. Done. Stiffness matrix: 'preconditioner' Calculating number of nonzero entries... Found 9870091 nonzero entries. Done. Getting parameter 'density': Found - Value = "0" Getting parameter 'alpha': Using default value = 0 Looking for essential component 'velocityRemesher': Found. Looking for essential component 'mesh-linear': Found. Looking for essential component 'VelocityField': Found. Looking for essential component 'VelocityField': Found. Looking for essential component 'velocity': Found. Looking for essential component 'PressureField': Found. Looking for essential component 'pressure': Found. Initialising Stg_Components from the live-component register Laying out material 'viscous' within Box 'boxShape': In func Materials_Register_AssignParticleProperties(): for swarm "materialSwarm" Assigning initial particle properties to the 10800000 global particles done 10% (1080000 particles)... done 20% (2160000 particles)... done 30% (3240000 particles)... done 40% (4320000 particles)... done 50% (5400000 particles)... done 60% (6480000 particles)... done 70% (7560000 particles)... done 80% (8640000 particles)... done 90% (9720000 particles)... done 100% (10800000 particles)... Materials_Register_AssignParticleProperties(): finished setup of material properties for swarm "materialSwarm" took 18.8151 [min] / 18.8151 [max] secs In func WeightsCalculator_CalculateAll(): for swarm "picIntegrationPoints" Calculating weights for the particles in the 360000 global cells done 10% (36000 cells)... done 20% (72000 cells)... done 30% (108000 cells)... done 40% (144000 cells)... done 50% (180000 cells)... done 60% (216000 cells)... done 70% (252000 cells)... done 80% (288000 cells)... done 90% (324000 cells)... done 100% (360000 cells)... WeightsCalculator_CalculateAll(): finished update of weights for swarm "picIntegrationPoints" Laying out material 'viscous' within Box 'boxShape': In func Materials_Register_AssignParticleProperties(): for swarm "gaussMaterialSwarm" Assigning initial particle properties to the 1 global particles done 100% (1 particles)... Materials_Register_AssignParticleProperties(): finished setup of material properties for swarm "gaussMaterialSwarm" took 3.91006e-05 [min] / 3.91006e-05 [max] secs TimeStep = 0, Time = 0 In SystemLinearEquations_NonLinearExecute Non linear solver - iteration 0 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.254598 Loop = 0, absResidual = 2.97830962e-01, relResidual = 5.04271259e-01 Loop = 1, absResidual = 1.90359720e-01, relResidual = 3.22306771e-01 Loop = 2, absResidual = 1.17113343e-01, relResidual = 1.98289972e-01 Loop = 3, absResidual = 8.10360573e-02, relResidual = 1.37205865e-01 Loop = 4, absResidual = 5.75611375e-02, relResidual = 9.74594014e-02 Loop = 5, absResidual = 3.66711889e-02, relResidual = 6.20896714e-02 Loop = 6, absResidual = 2.42660178e-02, relResidual = 4.10859075e-02 Loop = 7, absResidual = 1.64104921e-02, relResidual = 2.77853566e-02 Loop = 8, absResidual = 1.01503306e-02, relResidual = 1.71859901e-02 Loop = 9, absResidual = 6.19248857e-03, relResidual = 1.04847864e-02 Loop = 10, absResidual = 3.78971057e-03, relResidual = 6.41653275e-03 Loop = 11, absResidual = 2.14621014e-03, relResidual = 3.63384680e-03 Loop = 12, absResidual = 1.19793683e-03, relResidual = 2.02828178e-03 Loop = 13, absResidual = 7.74890663e-04, relResidual = 1.31200291e-03 Loop = 14, absResidual = 5.18223279e-04, relResidual = 8.77427597e-04 Loop = 15, absResidual = 3.13012352e-04, relResidual = 5.29975567e-04 Loop = 16, absResidual = 1.84595175e-04, relResidual = 3.12546555e-04 Loop = 17, absResidual = 1.15391893e-04, relResidual = 1.95375305e-04 Loop = 18, absResidual = 7.10036393e-05, relResidual = 1.20219518e-04 Loop = 19, absResidual = 4.33989557e-05, relResidual = 7.34807619e-05 Loop = 20, absResidual = 2.78915762e-05, relResidual = 4.72245066e-05 Loop = 21, absResidual = 1.85346855e-05, relResidual = 3.13819260e-05 Loop = 22, absResidual = 1.24018158e-05, relResidual = 2.09980831e-05 Loop = 23, absResidual = 8.56949235e-06, relResidual = 1.45094005e-05 Loop = 24, absResidual = 6.06759154e-06, relResidual = 1.02733175e-05 Loop = 25, absResidual = 4.25428732e-06, relResidual = 7.20312894e-06 Summary: Uzawa its. = 0026 , Uzawa residual = 7.2031289446049e-06 |G^T u|/|u| = 5.09487658e+00 |f - K u - G p|/|f| = 3.36952336e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252001e+00 , u_rms = 4.70271704e-01 |p|_{\infty} = 4.04644557e-01 , p_rms = 8.18002514e-03 min/max(u) = -1.00252001e+00 [63637] / 1.58662087e-02 [12] min/max(p) = -6.15045708e-02 [376196] / 4.04644557e-01 [371371] \sum_i p_i = 1.06763582e+03 Linear solver (stokesEqn-execute), solution time 2.483332e+03 (secs) Non linear solver - iteration 1 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094855 Loop = 0, absResidual = 6.57134316e-06, relResidual = 1.11262424e-05 Loop = 1, absResidual = 5.02048116e-06, relResidual = 8.50040688e-06 Summary: Uzawa its. = 0002 , Uzawa residual = 8.5004068762152e-06 |G^T u|/|u| = 5.09485970e+00 |f - K u - G p|/|f| = 3.48695435e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252051e+00 , u_rms = 4.70273347e-01 |p|_{\infty} = 4.04644615e-01 , p_rms = 8.18002237e-03 min/max(u) = -1.00252051e+00 [63637] / 1.58641471e-02 [12] min/max(p) = -6.15045993e-02 [376196] / 4.04644615e-01 [371371] \sum_i p_i = 1.06763412e+03 Linear solver (stokesEqn-execute), solution time 3.587157e+02 (secs) In func SystemLinearEquations_NonLinearExecute: Iteration 1 of 500 - Residual 6.8313e-06 - Tolerance = 0.01 Non linear solver - Residual 6.83125060e-06; Tolerance 1.0000e-02 - Converged - 3.287704e+03 (secs) In func SystemLinearEquations_NonLinearExecute: Converged after 1 iterations. Time Integration Time Integration EP: EulerDeform_IntegrationSetup - 0.2569 [min] / 0.2569 [max] (secs) EP: SwarmAdvector_AdvectionSetup - 0.0176 [min] / 0.0176 [max] (secs) 1st order: materialSwarmAdvector - 265.8846 [min] / 265.8846 [max] (secs) 1st order: EulerDeform_Velocity - 0.0740 [min] / 0.0740 [max] (secs) EP: EulerDeform_Execute - 19.5150 (secs) EP: SwarmAdvector_AdvectionFinish - 246.6321 (secs) EP: MaterialPointsSwarm_Update - 8.6440 (secs) In func WeightsCalculator_CalculateAll(): for swarm "picIntegrationPoints" Calculating weights for the particles in the 360000 global cells done 10% (36000 cells)... done 20% (72000 cells)... done 30% (108000 cells)... done 40% (144000 cells)... done 50% (180000 cells)... done 60% (216000 cells)... done 70% (252000 cells)... done 80% (288000 cells)... done 90% (324000 cells)... done 100% (360000 cells)... WeightsCalculator_CalculateAll(): finished update of weights for swarm "picIntegrationPoints" EP: IntegrationPointsSwarm_Update - 184.7581 (secs) EP: IntegrationPointsSwarm_Update - 0.0508 (secs) Time Integration - 725.848 [min] / 725.848 [max] (secs) Time Integration - 725.9 [min] / 725.9 [max] (secs) TimeStep = 1, Time = 3.74057 In SystemLinearEquations_NonLinearExecute Non linear solver - iteration 0 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094877 Loop = 0, absResidual = 5.24365005e-04, relResidual = 8.87906782e-04 Loop = 1, absResidual = 3.47013380e-04, relResidual = 5.87597438e-04 Loop = 2, absResidual = 2.44703397e-04, relResidual = 4.14356038e-04 Loop = 3, absResidual = 1.69808799e-04, relResidual = 2.87537083e-04 Loop = 4, absResidual = 1.14363930e-04, relResidual = 1.93652338e-04 Loop = 5, absResidual = 7.75192662e-05, relResidual = 1.31263302e-04 Loop = 6, absResidual = 5.22919232e-05, relResidual = 8.85458656e-05 Loop = 7, absResidual = 3.52158413e-05, relResidual = 5.96309518e-05 Loop = 8, absResidual = 2.37242331e-05, relResidual = 4.01722222e-05 Loop = 9, absResidual = 1.59535769e-05, relResidual = 2.70141772e-05 Loop = 10, absResidual = 1.07221948e-05, relResidual = 1.81558827e-05 Loop = 11, absResidual = 7.13611673e-06, relResidual = 1.20835799e-05 Loop = 12, absResidual = 4.67090156e-06, relResidual = 7.90923333e-06 Summary: Uzawa its. = 0013 , Uzawa residual = 7.9092333282479e-06 |G^T u|/|u| = 5.09486495e+00 |f - K u - G p|/|f| = 3.48422877e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252146e+00 , u_rms = 4.70273265e-01 |p|_{\infty} = 4.04580241e-01 , p_rms = 8.17557600e-03 min/max(u) = -1.00252146e+00 [61501] / 1.58569680e-02 [13546] min/max(p) = -6.15384817e-02 [370554] / 4.04580241e-01 [370188] \sum_i p_i = 1.06755947e+03 Linear solver (stokesEqn-execute), solution time 1.094831e+03 (secs) Non linear solver - iteration 1 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094865 Loop = 0, absResidual = 3.50350646e-06, relResidual = 5.93248428e-06 Summary: Uzawa its. = 0001 , Uzawa residual = 5.9324842830361e-06 |G^T u|/|u| = 5.09486496e+00 |f - K u - G p|/|f| = 3.48429528e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252146e+00 , u_rms = 4.70273263e-01 |p|_{\infty} = 4.04580422e-01 , p_rms = 8.17557490e-03 min/max(u) = -1.00252146e+00 [61501] / 1.58569919e-02 [13546] min/max(p) = -6.15386502e-02 [370554] / 4.04580422e-01 [370188] \sum_i p_i = 1.06755907e+03 Linear solver (stokesEqn-execute), solution time 2.431509e+02 (secs) In func SystemLinearEquations_NonLinearExecute: Iteration 1 of 500 - Residual 1.1341e-07 - Tolerance = 0.01 Non linear solver - Residual 1.13407737e-07; Tolerance 1.0000e-02 - Converged - 1.789368e+03 (secs) In func SystemLinearEquations_NonLinearExecute: Converged after 1 iterations. Time Integration EP: EulerDeform_IntegrationSetup - 0.0811 [min] / 0.0811 [max] (secs) EP: SwarmAdvector_AdvectionSetup - 0.0000 [min] / 0.0000 [max] (secs) 1st order: materialSwarmAdvector - 268.5346 [min] / 268.5346 [max] (secs) 1st order: EulerDeform_Velocity - 0.0439 [min] / 0.0439 [max] (secs) EP: EulerDeform_Execute - 20.1018 (secs) EP: SwarmAdvector_AdvectionFinish - 247.3663 (secs) EP: MaterialPointsSwarm_Update - 8.8303 (secs) In func WeightsCalculator_CalculateAll(): for swarm "picIntegrationPoints" Calculating weights for the particles in the 360000 global cells done 10% (36000 cells)... done 20% (72000 cells)... done 30% (108000 cells)... done 40% (144000 cells)... done 50% (180000 cells)... done 60% (216000 cells)... done 70% (252000 cells)... done 80% (288000 cells)... done 90% (324000 cells)... done 100% (360000 cells)... WeightsCalculator_CalculateAll(): finished update of weights for swarm "picIntegrationPoints" EP: IntegrationPointsSwarm_Update - 182.5954 (secs) EP: IntegrationPointsSwarm_Update - 0.0000 (secs) Time Integration - 727.6 [min] / 727.6 [max] (secs) TimeStep = 2, Time = 7.48114 In SystemLinearEquations_NonLinearExecute Non linear solver - iteration 0 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094860 Loop = 0, absResidual = 5.24679665e-04, relResidual = 8.88406894e-04 Loop = 1, absResidual = 3.42156137e-04, relResidual = 5.79351348e-04 Loop = 2, absResidual = 2.38828350e-04, relResidual = 4.04392941e-04 Loop = 3, absResidual = 1.62191493e-04, relResidual = 2.74628597e-04 Loop = 4, absResidual = 1.07605481e-04, relResidual = 1.82201555e-04 Loop = 5, absResidual = 7.20343572e-05, relResidual = 1.21971221e-04 Loop = 6, absResidual = 4.80177695e-05, relResidual = 8.13054522e-05 Loop = 7, absResidual = 3.21153592e-05, relResidual = 5.43789067e-05 Loop = 8, absResidual = 2.14235072e-05, relResidual = 3.62750699e-05 Loop = 9, absResidual = 1.40696917e-05, relResidual = 2.38233191e-05 Loop = 10, absResidual = 9.19402007e-06, relResidual = 1.55676527e-05 Loop = 11, absResidual = 5.96702819e-06, relResidual = 1.01035915e-05 Loop = 12, absResidual = 3.83714567e-06, relResidual = 6.49719609e-06 Summary: Uzawa its. = 0013 , Uzawa residual = 6.4971960883126e-06 |G^T u|/|u| = 5.09486274e+00 |f - K u - G p|/|f| = 3.48351237e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252345e+00 , u_rms = 4.70273368e-01 |p|_{\infty} = 4.04488461e-01 , p_rms = 8.17661571e-03 min/max(u) = -1.00252345e+00 [68710] / 1.58548945e-02 [25294] min/max(p) = -6.15473862e-02 [367915] / 4.04488461e-01 [375011] \sum_i p_i = 1.06757395e+03 Linear solver (stokesEqn-execute), solution time 1.009253e+03 (secs) Non linear solver - iteration 1 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094863 Loop = 0, absResidual = 2.85067458e-06, relResidual = 4.82686697e-06 Summary: Uzawa its. = 0001 , Uzawa residual = 4.8268669724353e-06 |G^T u|/|u| = 5.09486274e+00 |f - K u - G p|/|f| = 3.48347267e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252346e+00 , u_rms = 4.70273369e-01 |p|_{\infty} = 4.04488669e-01 , p_rms = 8.17661618e-03 min/max(u) = -1.00252346e+00 [68710] / 1.58548920e-02 [25294] min/max(p) = -6.15476669e-02 [367915] / 4.04488669e-01 [375011] \sum_i p_i = 1.06757407e+03 Linear solver (stokesEqn-execute), solution time 2.365018e+02 (secs) In func SystemLinearEquations_NonLinearExecute: Iteration 1 of 500 - Residual 6.0194e-08 - Tolerance = 0.01 Non linear solver - Residual 6.01936812e-08; Tolerance 1.0000e-02 - Converged - 1.699186e+03 (secs) In func SystemLinearEquations_NonLinearExecute: Converged after 1 iterations. Time Integration EP: EulerDeform_IntegrationSetup - 0.0449 [min] / 0.0449 [max] (secs) EP: SwarmAdvector_AdvectionSetup - 0.0000 [min] / 0.0000 [max] (secs) 1st order: materialSwarmAdvector - 269.7257 [min] / 269.7257 [max] (secs) 1st order: EulerDeform_Velocity - 0.0444 [min] / 0.0444 [max] (secs) EP: EulerDeform_Execute - 20.1944 (secs) EP: SwarmAdvector_AdvectionFinish - 248.5865 (secs) EP: MaterialPointsSwarm_Update - 8.8264 (secs) In func WeightsCalculator_CalculateAll(): for swarm "picIntegrationPoints" Calculating weights for the particles in the 360000 global cells done 10% (36000 cells)... done 20% (72000 cells)... done 30% (108000 cells)... done 40% (144000 cells)... done 50% (180000 cells)... done 60% (216000 cells)... done 70% (252000 cells)... done 80% (288000 cells)... done 90% (324000 cells)... done 100% (360000 cells)... WeightsCalculator_CalculateAll(): finished update of weights for swarm "picIntegrationPoints" EP: IntegrationPointsSwarm_Update - 178.5731 (secs) EP: IntegrationPointsSwarm_Update - 0.0000 (secs) Time Integration - 726 [min] / 726 [max] (secs) TimeStep = 3, Time = 11.2217 In SystemLinearEquations_NonLinearExecute Non linear solver - iteration 0 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094871 Loop = 0, absResidual = 5.00733125e-04, relResidual = 8.47911019e-04 Loop = 1, absResidual = 3.12140783e-04, relResidual = 5.28560217e-04 Loop = 2, absResidual = 2.10305056e-04, relResidual = 3.56117791e-04 Loop = 3, absResidual = 1.40721137e-04, relResidual = 2.38288614e-04 Loop = 4, absResidual = 9.46210982e-05, relResidual = 1.60225613e-04 Loop = 5, absResidual = 6.41771227e-05, relResidual = 1.08673636e-04 Loop = 6, absResidual = 4.28366716e-05, relResidual = 7.25370143e-05 Loop = 7, absResidual = 2.85821825e-05, relResidual = 4.83993294e-05 Loop = 8, absResidual = 1.91927734e-05, relResidual = 3.24998752e-05 Loop = 9, absResidual = 1.29815619e-05, relResidual = 2.19821874e-05 Loop = 10, absResidual = 8.71103951e-06, relResidual = 1.47507445e-05 Loop = 11, absResidual = 5.75111469e-06, relResidual = 9.73858782e-06 Summary: Uzawa its. = 0012 , Uzawa residual = 9.7385878230786e-06 |G^T u|/|u| = 5.09486553e+00 |f - K u - G p|/|f| = 3.48212751e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252201e+00 , u_rms = 4.70273304e-01 |p|_{\infty} = 4.04386554e-01 , p_rms = 8.17392146e-03 min/max(u) = -1.00252201e+00 [72181] / 1.58542071e-02 [2065] min/max(p) = -6.16548388e-02 [368916] / 4.04386554e-01 [371644] \sum_i p_i = 1.06753994e+03 Linear solver (stokesEqn-execute), solution time 9.734394e+02 (secs) Non linear solver - iteration 1 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094866 Loop = 0, absResidual = 4.30526906e-06, relResidual = 7.29028077e-06 Summary: Uzawa its. = 0001 , Uzawa residual = 7.2902807716469e-06 |G^T u|/|u| = 5.09486553e+00 |f - K u - G p|/|f| = 3.48220197e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252201e+00 , u_rms = 4.70273304e-01 |p|_{\infty} = 4.04386226e-01 , p_rms = 8.17392035e-03 min/max(u) = -1.00252201e+00 [72181] / 1.58542102e-02 [2065] min/max(p) = -6.16556356e-02 [368916] / 4.04386226e-01 [371644] \sum_i p_i = 1.06753963e+03 Linear solver (stokesEqn-execute), solution time 2.409279e+02 (secs) In func SystemLinearEquations_NonLinearExecute: Iteration 1 of 500 - Residual 1.0671e-07 - Tolerance = 0.01 Non linear solver - Residual 1.06711618e-07; Tolerance 1.0000e-02 - Converged - 1.670002e+03 (secs) In func SystemLinearEquations_NonLinearExecute: Converged after 1 iterations. Time Integration EP: EulerDeform_IntegrationSetup - 0.0448 [min] / 0.0448 [max] (secs) EP: SwarmAdvector_AdvectionSetup - 0.0000 [min] / 0.0000 [max] (secs) 1st order: materialSwarmAdvector - 270.4601 [min] / 270.4601 [max] (secs) 1st order: EulerDeform_Velocity - 0.0437 [min] / 0.0437 [max] (secs) EP: EulerDeform_Execute - 20.2127 (secs) EP: SwarmAdvector_AdvectionFinish - 249.4426 (secs) EP: MaterialPointsSwarm_Update - 8.8470 (secs) In func WeightsCalculator_CalculateAll(): for swarm "picIntegrationPoints" Calculating weights for the particles in the 360000 global cells done 10% (36000 cells)... done 20% (72000 cells)... done 30% (108000 cells)... done 40% (144000 cells)... done 50% (180000 cells)... done 60% (216000 cells)... done 70% (252000 cells)... done 80% (288000 cells)... done 90% (324000 cells)... done 100% (360000 cells)... WeightsCalculator_CalculateAll(): finished update of weights for swarm "picIntegrationPoints" EP: IntegrationPointsSwarm_Update - 176.7632 (secs) EP: IntegrationPointsSwarm_Update - 0.0000 (secs) Time Integration - 725.8 [min] / 725.8 [max] (secs) TimeStep = 4, Time = 14.9623 In SystemLinearEquations_NonLinearExecute Non linear solver - iteration 0 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094866 Loop = 0, absResidual = 4.55994498e-04, relResidual = 7.72149630e-04 Loop = 1, absResidual = 2.77538928e-04, relResidual = 4.69965276e-04 Loop = 2, absResidual = 1.86096176e-04, relResidual = 3.15122429e-04 Loop = 3, absResidual = 1.25773002e-04, relResidual = 2.12975326e-04 Loop = 4, absResidual = 8.44564686e-05, relResidual = 1.43012758e-04 Loop = 5, absResidual = 5.64164177e-05, relResidual = 9.55316702e-05 Loop = 6, absResidual = 3.68117288e-05, relResidual = 6.23344424e-05 Loop = 7, absResidual = 2.38068588e-05, relResidual = 4.03128925e-05 Loop = 8, absResidual = 1.55073740e-05, relResidual = 2.62591175e-05 Loop = 9, absResidual = 1.01732995e-05, relResidual = 1.72267636e-05 Loop = 10, absResidual = 6.67565176e-06, relResidual = 1.13040882e-05 Loop = 11, absResidual = 4.34305977e-06, relResidual = 7.35423784e-06 Summary: Uzawa its. = 0012 , Uzawa residual = 7.3542378355806e-06 |G^T u|/|u| = 5.09486470e+00 |f - K u - G p|/|f| = 3.47768198e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252038e+00 , u_rms = 4.70273416e-01 |p|_{\infty} = 4.04393644e-01 , p_rms = 8.17403656e-03 min/max(u) = -1.00252038e+00 [72181] / 1.58533068e-02 [7405] min/max(p) = -6.16386518e-02 [369553] / 4.04393644e-01 [376740] \sum_i p_i = 1.06753719e+03 Linear solver (stokesEqn-execute), solution time 9.382749e+02 (secs) Non linear solver - iteration 1 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094865 Loop = 0, absResidual = 3.24581429e-06, relResidual = 5.49623802e-06 Summary: Uzawa its. = 0001 , Uzawa residual = 5.4962380222878e-06 |G^T u|/|u| = 5.09486469e+00 |f - K u - G p|/|f| = 3.47760482e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252038e+00 , u_rms = 4.70273416e-01 |p|_{\infty} = 4.04392975e-01 , p_rms = 8.17403632e-03 min/max(u) = -1.00252038e+00 [72181] / 1.58532794e-02 [7405] min/max(p) = -6.16383376e-02 [369553] / 4.04392975e-01 [376740] \sum_i p_i = 1.06753715e+03 Linear solver (stokesEqn-execute), solution time 2.373994e+02 (secs) In func SystemLinearEquations_NonLinearExecute: Iteration 1 of 500 - Residual 6.0391e-08 - Tolerance = 0.01 Non linear solver - Residual 6.03905865e-08; Tolerance 1.0000e-02 - Converged - 1.636663e+03 (secs) In func SystemLinearEquations_NonLinearExecute: Converged after 1 iterations. Time Integration EP: EulerDeform_IntegrationSetup - 0.0451 [min] / 0.0451 [max] (secs) EP: SwarmAdvector_AdvectionSetup - 0.0000 [min] / 0.0000 [max] (secs) 1st order: materialSwarmAdvector - 271.2931 [min] / 271.2931 [max] (secs) 1st order: EulerDeform_Velocity - 0.0448 [min] / 0.0448 [max] (secs) EP: EulerDeform_Execute - 20.1336 (secs) EP: SwarmAdvector_AdvectionFinish - 250.4052 (secs) EP: MaterialPointsSwarm_Update - 8.7950 (secs) In func WeightsCalculator_CalculateAll(): for swarm "picIntegrationPoints" Calculating weights for the particles in the 360000 global cells done 10% (36000 cells)... done 20% (72000 cells)... done 30% (108000 cells)... done 40% (144000 cells)... done 50% (180000 cells)... done 60% (216000 cells)... done 70% (252000 cells)... done 80% (288000 cells)... done 90% (324000 cells)... done 100% (360000 cells)... WeightsCalculator_CalculateAll(): finished update of weights for swarm "picIntegrationPoints" EP: IntegrationPointsSwarm_Update - 175.4320 (secs) EP: IntegrationPointsSwarm_Update - 0.0000 (secs) Time Integration - 726.1 [min] / 726.1 [max] (secs) TimeStep = 5, Time = 18.7028 In SystemLinearEquations_NonLinearExecute Non linear solver - iteration 0 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094863 Loop = 0, absResidual = 4.94740830e-04, relResidual = 8.37748497e-04 Loop = 1, absResidual = 3.14035288e-04, relResidual = 5.31758397e-04 Loop = 2, absResidual = 2.12070339e-04, relResidual = 3.59100355e-04 Loop = 3, absResidual = 1.40244173e-04, relResidual = 2.37476549e-04 Loop = 4, absResidual = 9.10331430e-05, relResidual = 1.54147129e-04 Loop = 5, absResidual = 5.91546207e-05, relResidual = 1.00166979e-04 Loop = 6, absResidual = 3.88467109e-05, relResidual = 6.57794377e-05 Loop = 7, absResidual = 2.58047665e-05, relResidual = 4.36954118e-05 Loop = 8, absResidual = 1.73483825e-05, relResidual = 2.93761510e-05 Loop = 9, absResidual = 1.15936171e-05, relResidual = 1.96315621e-05 Loop = 10, absResidual = 7.56415736e-06, relResidual = 1.28084465e-05 Loop = 11, absResidual = 4.80519050e-06, relResidual = 8.13666646e-06 Summary: Uzawa its. = 0012 , Uzawa residual = 8.1366664614475e-06 |G^T u|/|u| = 5.09486582e+00 |f - K u - G p|/|f| = 3.48293413e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252139e+00 , u_rms = 4.70273214e-01 |p|_{\infty} = 4.04386219e-01 , p_rms = 8.17465330e-03 min/max(u) = -1.00252139e+00 [61768] / 1.58582875e-02 [12] min/max(p) = -6.16113300e-02 [372647] / 4.04386219e-01 [369824] \sum_i p_i = 1.06754339e+03 Linear solver (stokesEqn-execute), solution time 9.052485e+02 (secs) Non linear solver - iteration 1 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094866 Loop = 0, absResidual = 3.51195520e-06, relResidual = 5.94682106e-06 Summary: Uzawa its. = 0001 , Uzawa residual = 5.9468210559967e-06 |G^T u|/|u| = 5.09486583e+00 |f - K u - G p|/|f| = 3.48303583e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252140e+00 , u_rms = 4.70273213e-01 |p|_{\infty} = 4.04386672e-01 , p_rms = 8.17465394e-03 min/max(u) = -1.00252140e+00 [61768] / 1.58583222e-02 [12] min/max(p) = -6.16116326e-02 [372647] / 4.04386672e-01 [369824] \sum_i p_i = 1.06754349e+03 Linear solver (stokesEqn-execute), solution time 2.351723e+02 (secs) In func SystemLinearEquations_NonLinearExecute: Iteration 1 of 500 - Residual 6.3224e-08 - Tolerance = 0.01 Non linear solver - Residual 6.32242775e-08; Tolerance 1.0000e-02 - Converged - 1.599141e+03 (secs) In func SystemLinearEquations_NonLinearExecute: Converged after 1 iterations. Time Integration EP: EulerDeform_IntegrationSetup - 0.0451 [min] / 0.0451 [max] (secs) EP: SwarmAdvector_AdvectionSetup - 0.0000 [min] / 0.0000 [max] (secs) 1st order: materialSwarmAdvector - 272.7115 [min] / 272.7115 [max] (secs) 1st order: EulerDeform_Velocity - 0.0438 [min] / 0.0438 [max] (secs) EP: EulerDeform_Execute - 20.2229 (secs) EP: SwarmAdvector_AdvectionFinish - 251.3280 (secs) EP: MaterialPointsSwarm_Update - 8.8777 (secs) In func WeightsCalculator_CalculateAll(): for swarm "picIntegrationPoints" Calculating weights for the particles in the 360000 global cells done 10% (36000 cells)... done 20% (72000 cells)... done 30% (108000 cells)... done 40% (144000 cells)... done 50% (180000 cells)... done 60% (216000 cells)... done 70% (252000 cells)... done 80% (288000 cells)... done 90% (324000 cells)... done 100% (360000 cells)... WeightsCalculator_CalculateAll(): finished update of weights for swarm "picIntegrationPoints" EP: IntegrationPointsSwarm_Update - 174.8932 (secs) EP: IntegrationPointsSwarm_Update - 0.0000 (secs) Time Integration - 728.1 [min] / 728.1 [max] (secs) TimeStep = 6, Time = 22.4434 In SystemLinearEquations_NonLinearExecute Non linear solver - iteration 0 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094866 Loop = 0, absResidual = 4.80592157e-04, relResidual = 8.13787837e-04 Loop = 1, absResidual = 3.08697492e-04, relResidual = 5.22718194e-04 Loop = 2, absResidual = 2.13852910e-04, relResidual = 3.62117639e-04 Loop = 3, absResidual = 1.45123480e-04, relResidual = 2.45737932e-04 Loop = 4, absResidual = 9.70590472e-05, relResidual = 1.64350314e-04 Loop = 5, absResidual = 6.43583408e-05, relResidual = 1.08978131e-04 Loop = 6, absResidual = 4.25295959e-05, relResidual = 7.20154654e-05 Loop = 7, absResidual = 2.82588137e-05, relResidual = 4.78507160e-05 Loop = 8, absResidual = 1.89032124e-05, relResidual = 3.20088542e-05 Loop = 9, absResidual = 1.26606502e-05, relResidual = 2.14383089e-05 Loop = 10, absResidual = 8.40299505e-06, relResidual = 1.42288114e-05 Loop = 11, absResidual = 5.44181153e-06, relResidual = 9.21463234e-06 Summary: Uzawa its. = 0012 , Uzawa residual = 9.2146323439682e-06 |G^T u|/|u| = 5.09486526e+00 |f - K u - G p|/|f| = 3.47930330e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252335e+00 , u_rms = 4.70273269e-01 |p|_{\infty} = 4.04772343e-01 , p_rms = 8.17477963e-03 min/max(u) = -1.00252335e+00 [53579] / 1.58540853e-02 [8206] min/max(p) = -6.16314451e-02 [370736] / 4.04772343e-01 [376740] \sum_i p_i = 1.06754388e+03 Linear solver (stokesEqn-execute), solution time 8.887495e+02 (secs) Non linear solver - iteration 1 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094865 Loop = 0, absResidual = 4.00215742e-06, relResidual = 6.77686262e-06 Summary: Uzawa its. = 0001 , Uzawa residual = 6.7768626227652e-06 |G^T u|/|u| = 5.09486526e+00 |f - K u - G p|/|f| = 3.47916440e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252336e+00 , u_rms = 4.70273270e-01 |p|_{\infty} = 4.04772257e-01 , p_rms = 8.17477931e-03 min/max(u) = -1.00252336e+00 [53579] / 1.58540786e-02 [8206] min/max(p) = -6.16321090e-02 [370736] / 4.04772257e-01 [376740] \sum_i p_i = 1.06754386e+03 Linear solver (stokesEqn-execute), solution time 2.331179e+02 (secs) In func SystemLinearEquations_NonLinearExecute: Iteration 1 of 500 - Residual 7.54e-08 - Tolerance = 0.01 Non linear solver - Residual 7.54001652e-08; Tolerance 1.0000e-02 - Converged - 1.579833e+03 (secs) In func SystemLinearEquations_NonLinearExecute: Converged after 1 iterations. Time Integration EP: EulerDeform_IntegrationSetup - 0.0450 [min] / 0.0450 [max] (secs) EP: SwarmAdvector_AdvectionSetup - 0.0000 [min] / 0.0000 [max] (secs) 1st order: materialSwarmAdvector - 273.5980 [min] / 273.5980 [max] (secs) 1st order: EulerDeform_Velocity - 0.0439 [min] / 0.0439 [max] (secs) EP: EulerDeform_Execute - 20.1852 (secs) EP: SwarmAdvector_AdvectionFinish - 252.0948 (secs) EP: MaterialPointsSwarm_Update - 8.7881 (secs) In func WeightsCalculator_CalculateAll(): for swarm "picIntegrationPoints" Calculating weights for the particles in the 360000 global cells done 10% (36000 cells)... done 20% (72000 cells)... done 30% (108000 cells)... done 40% (144000 cells)... done 50% (180000 cells)... done 60% (216000 cells)... done 70% (252000 cells)... done 80% (288000 cells)... done 90% (324000 cells)... done 100% (360000 cells)... WeightsCalculator_CalculateAll(): finished update of weights for swarm "picIntegrationPoints" EP: IntegrationPointsSwarm_Update - 174.7945 (secs) EP: IntegrationPointsSwarm_Update - 0.0000 (secs) Time Integration - 729.6 [min] / 729.6 [max] (secs) TimeStep = 7, Time = 26.184 In SystemLinearEquations_NonLinearExecute Non linear solver - iteration 0 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094860 Loop = 0, absResidual = 5.41321759e-04, relResidual = 9.16592639e-04 Loop = 1, absResidual = 3.57724172e-04, relResidual = 6.05716171e-04 Loop = 2, absResidual = 2.51496339e-04, relResidual = 4.25845976e-04 Loop = 3, absResidual = 1.71175881e-04, relResidual = 2.89843426e-04 Loop = 4, absResidual = 1.13189031e-04, relResidual = 1.91657237e-04 Loop = 5, absResidual = 7.41325232e-05, relResidual = 1.25524836e-04 Loop = 6, absResidual = 4.81717433e-05, relResidual = 8.15667662e-05 Loop = 7, absResidual = 3.15485548e-05, relResidual = 5.34195653e-05 Loop = 8, absResidual = 2.09705060e-05, relResidual = 3.55082927e-05 Loop = 9, absResidual = 1.41699843e-05, relResidual = 2.39933147e-05 Loop = 10, absResidual = 9.52783911e-06, relResidual = 1.61330060e-05 Loop = 11, absResidual = 6.31188888e-06, relResidual = 1.06876008e-05 Loop = 12, absResidual = 4.10390776e-06, relResidual = 6.94893856e-06 Summary: Uzawa its. = 0013 , Uzawa residual = 6.9489385551243e-06 |G^T u|/|u| = 5.09486434e+00 |f - K u - G p|/|f| = 3.47969205e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252043e+00 , u_rms = 4.70273230e-01 |p|_{\infty} = 4.04611633e-01 , p_rms = 8.17665033e-03 min/max(u) = -1.00252043e+00 [65239] / 1.58616032e-02 [17551] min/max(p) = -6.16417106e-02 [371919] / 4.04611633e-01 [371462] \sum_i p_i = 1.06757369e+03 Linear solver (stokesEqn-execute), solution time 1.053741e+03 (secs) Non linear solver - iteration 1 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094864 Loop = 0, absResidual = 3.05480780e-06, relResidual = 5.17255088e-06 Summary: Uzawa its. = 0001 , Uzawa residual = 5.1725508794043e-06 |G^T u|/|u| = 5.09486433e+00 |f - K u - G p|/|f| = 3.47956879e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252043e+00 , u_rms = 4.70273231e-01 |p|_{\infty} = 4.04611700e-01 , p_rms = 8.17665121e-03 min/max(u) = -1.00252043e+00 [65239] / 1.58615917e-02 [17551] min/max(p) = -6.16416566e-02 [371919] / 4.04611700e-01 [371462] \sum_i p_i = 1.06757388e+03 Linear solver (stokesEqn-execute), solution time 2.422166e+02 (secs) In func SystemLinearEquations_NonLinearExecute: Iteration 1 of 500 - Residual 7.2333e-08 - Tolerance = 0.01 Non linear solver - Residual 7.23332665e-08; Tolerance 1.0000e-02 - Converged - 1.756101e+03 (secs) In func SystemLinearEquations_NonLinearExecute: Converged after 1 iterations. Time Integration EP: EulerDeform_IntegrationSetup - 0.0451 [min] / 0.0451 [max] (secs) EP: SwarmAdvector_AdvectionSetup - 0.0000 [min] / 0.0000 [max] (secs) 1st order: materialSwarmAdvector - 274.1201 [min] / 274.1201 [max] (secs) 1st order: EulerDeform_Velocity - 0.0436 [min] / 0.0436 [max] (secs) EP: EulerDeform_Execute - 20.0812 (secs) EP: SwarmAdvector_AdvectionFinish - 252.7205 (secs) EP: MaterialPointsSwarm_Update - 8.7365 (secs) In func WeightsCalculator_CalculateAll(): for swarm "picIntegrationPoints" Calculating weights for the particles in the 360000 global cells done 10% (36000 cells)... done 20% (72000 cells)... done 30% (108000 cells)... done 40% (144000 cells)... done 50% (180000 cells)... done 60% (216000 cells)... done 70% (252000 cells)... done 80% (288000 cells)... done 90% (324000 cells)... done 100% (360000 cells)... WeightsCalculator_CalculateAll(): finished update of weights for swarm "picIntegrationPoints" EP: IntegrationPointsSwarm_Update - 173.1889 (secs) EP: IntegrationPointsSwarm_Update - 0.0000 (secs) Time Integration - 728.9 [min] / 728.9 [max] (secs) TimeStep = 8, Time = 29.9245 In SystemLinearEquations_NonLinearExecute Non linear solver - iteration 0 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094872 Loop = 0, absResidual = 5.01956175e-04, relResidual = 8.49971776e-04 Loop = 1, absResidual = 3.27249129e-04, relResidual = 5.54137069e-04 Loop = 2, absResidual = 2.31144423e-04, relResidual = 3.91401174e-04 Loop = 3, absResidual = 1.57854226e-04, relResidual = 2.67297513e-04 Loop = 4, absResidual = 1.05287548e-04, relResidual = 1.78285373e-04 Loop = 5, absResidual = 7.09369566e-05, relResidual = 1.20118875e-04 Loop = 6, absResidual = 4.76104961e-05, relResidual = 8.06197432e-05 Loop = 7, absResidual = 3.19498680e-05, relResidual = 5.41013088e-05 Loop = 8, absResidual = 2.15542789e-05, relResidual = 3.64982634e-05 Loop = 9, absResidual = 1.44736571e-05, relResidual = 2.45085142e-05 Loop = 10, absResidual = 9.58042833e-06, relResidual = 1.62227184e-05 Loop = 11, absResidual = 6.28223971e-06, relResidual = 1.06378340e-05 Loop = 12, absResidual = 4.07066235e-06, relResidual = 6.89292866e-06 Summary: Uzawa its. = 0013 , Uzawa residual = 6.8929286606817e-06 |G^T u|/|u| = 5.09486558e+00 |f - K u - G p|/|f| = 3.48100386e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252272e+00 , u_rms = 4.70273293e-01 |p|_{\infty} = 4.04516043e-01 , p_rms = 8.17425685e-03 min/max(u) = -1.00252272e+00 [69778] / 1.58535497e-02 [3400] min/max(p) = -6.18343568e-02 [376742] / 4.04516043e-01 [370643] \sum_i p_i = 1.06753943e+03 Linear solver (stokesEqn-execute), solution time 1.083000e+03 (secs) Non linear solver - iteration 1 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094866 Loop = 0, absResidual = 3.04554854e-06, relResidual = 5.15708428e-06 Summary: Uzawa its. = 0001 , Uzawa residual = 5.1570842803466e-06 |G^T u|/|u| = 5.09486559e+00 |f - K u - G p|/|f| = 3.48106845e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00252271e+00 , u_rms = 4.70273292e-01 |p|_{\infty} = 4.04515981e-01 , p_rms = 8.17425626e-03 min/max(u) = -1.00252271e+00 [69778] / 1.58535732e-02 [3400] min/max(p) = -6.18342110e-02 [376742] / 4.04515981e-01 [370643] \sum_i p_i = 1.06753921e+03 Linear solver (stokesEqn-execute), solution time 2.409225e+02 (secs) In func SystemLinearEquations_NonLinearExecute: Iteration 1 of 500 - Residual 8.884e-08 - Tolerance = 0.01 Non linear solver - Residual 8.88396405e-08; Tolerance 1.0000e-02 - Converged - 1.785409e+03 (secs) In func SystemLinearEquations_NonLinearExecute: Converged after 1 iterations. Time Integration EP: EulerDeform_IntegrationSetup - 0.1338 [min] / 0.1338 [max] (secs) EP: SwarmAdvector_AdvectionSetup - 0.0000 [min] / 0.0000 [max] (secs) 1st order: materialSwarmAdvector - 274.6229 [min] / 274.6229 [max] (secs) 1st order: EulerDeform_Velocity - 0.0437 [min] / 0.0437 [max] (secs) EP: EulerDeform_Execute - 20.1463 (secs) EP: SwarmAdvector_AdvectionFinish - 254.2913 (secs) EP: MaterialPointsSwarm_Update - 8.8022 (secs) In func WeightsCalculator_CalculateAll(): for swarm "picIntegrationPoints" Calculating weights for the particles in the 360000 global cells done 10% (36000 cells)... done 20% (72000 cells)... done 30% (108000 cells)... done 40% (144000 cells)... done 50% (180000 cells)... done 60% (216000 cells)... done 70% (252000 cells)... done 80% (288000 cells)... done 90% (324000 cells)... done 100% (360000 cells)... WeightsCalculator_CalculateAll(): finished update of weights for swarm "picIntegrationPoints" EP: IntegrationPointsSwarm_Update - 172.4401 (secs) EP: IntegrationPointsSwarm_Update - 0.0000 (secs) Time Integration - 730.5 [min] / 730.5 [max] (secs) TimeStep = 9, Time = 33.6651 In SystemLinearEquations_NonLinearExecute Non linear solver - iteration 0 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094861 Loop = 0, absResidual = 5.42388866e-04, relResidual = 9.18413859e-04 Loop = 1, absResidual = 3.54968975e-04, relResidual = 6.01060322e-04 Loop = 2, absResidual = 2.50335810e-04, relResidual = 4.23887531e-04 Loop = 3, absResidual = 1.70196021e-04, relResidual = 2.88188779e-04 Loop = 4, absResidual = 1.13290195e-04, relResidual = 1.91831529e-04 Loop = 5, absResidual = 7.50833379e-05, relResidual = 1.27136788e-04 Loop = 6, absResidual = 4.90270629e-05, relResidual = 8.30163319e-05 Loop = 7, absResidual = 3.22481501e-05, relResidual = 5.46050073e-05 Loop = 8, absResidual = 2.13837428e-05, relResidual = 3.62085710e-05 Loop = 9, absResidual = 1.41485805e-05, relResidual = 2.39574469e-05 Loop = 10, absResidual = 9.30283563e-06, relResidual = 1.57522650e-05 Loop = 11, absResidual = 6.13101797e-06, relResidual = 1.03815034e-05 Loop = 12, absResidual = 4.05903433e-06, relResidual = 6.87306399e-06 Summary: Uzawa its. = 0013 , Uzawa residual = 6.8730639932630e-06 |G^T u|/|u| = 5.09486404e+00 |f - K u - G p|/|f| = 3.48025732e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00251954e+00 , u_rms = 4.70273290e-01 |p|_{\infty} = 4.04538044e-01 , p_rms = 8.17574837e-03 min/max(u) = -1.00251954e+00 [64705] / 1.58511972e-02 [6871] min/max(p) = -6.17478587e-02 [369826] / 4.04538044e-01 [370734] \sum_i p_i = 1.06755810e+03 Linear solver (stokesEqn-execute), solution time 1.054940e+03 (secs) Non linear solver - iteration 1 Linear solver (stokesEqn-execute) Relative residual < 1e-05 for Uzawa stopping condition Initial l2Norm( Div u ) / l2Norm( u ) = 5.094864 Loop = 0, absResidual = 3.05894316e-06, relResidual = 5.17963397e-06 Summary: Uzawa its. = 0001 , Uzawa residual = 5.1796339710273e-06 |G^T u|/|u| = 5.09486403e+00 |f - K u - G p|/|f| = 3.48020365e-07 |f - K u - G p|_w/|f|_w = 1.00000000e+00 |u|_{\infty} = 1.00251953e+00 , u_rms = 4.70273290e-01 |p|_{\infty} = 4.04537897e-01 , p_rms = 8.17574898e-03 min/max(u) = -1.00251953e+00 [64705] / 1.58511752e-02 [6871] min/max(p) = -6.17478255e-02 [369826] / 4.04537897e-01 [370734] \sum_i p_i = 1.06755824e+03 Linear solver (stokesEqn-execute), solution time 2.410544e+02 (secs) In func SystemLinearEquations_NonLinearExecute: Iteration 1 of 500 - Residual 7.4206e-08 - Tolerance = 0.01 Non linear solver - Residual 7.42058786e-08; Tolerance 1.0000e-02 - Converged - 1.762327e+03 (secs) In func SystemLinearEquations_NonLinearExecute: Converged after 1 iterations. Time Integration EP: EulerDeform_IntegrationSetup - 0.0621 [min] / 0.0621 [max] (secs) EP: SwarmAdvector_AdvectionSetup - 0.0000 [min] / 0.0000 [max] (secs) 1st order: materialSwarmAdvector - 275.6823 [min] / 275.6823 [max] (secs) 1st order: EulerDeform_Velocity - 0.0443 [min] / 0.0443 [max] (secs) EP: EulerDeform_Execute - 20.1607 (secs) EP: SwarmAdvector_AdvectionFinish - 254.6320 (secs) EP: MaterialPointsSwarm_Update - 8.8646 (secs) In func WeightsCalculator_CalculateAll(): for swarm "picIntegrationPoints" Calculating weights for the particles in the 360000 global cells done 10% (36000 cells)... done 20% (72000 cells)... done 30% (108000 cells)... Error in _PCDVC_Calculate3D: Problem has an under resolved cell (Cell Id = 125459), add more particles to your model *** glibc detected *** ./Gale-1_5_0: double free or corruption (out): 0x00000000a035d320 *** ======= Backtrace: ========= /lib64/libc.so.6[0x33b4471ce2] /lib64/libc.so.6(cfree+0x8c)[0x33b447590c] ./Gale-1_5_0[0x4ddd71] ./Gale-1_5_0[0x4d922a] ./Gale-1_5_0(_PCDVC_Calculate+0x2a)[0x4da130] ./Gale-1_5_0[0x4c7509] ./Gale-1_5_0[0x4c7f4e] ./Gale-1_5_0[0x4c0e38] ./Gale-1_5_0(_IntegrationPointsSwarm_UpdateHook+0x28)[0x4c0f3b] ./Gale-1_5_0[0x579023] ./Gale-1_5_0(_TimeIntegrator_ExecuteEuler+0x156)[0x579443] ./Gale-1_5_0(TimeIntegrator_UpdateClass+0x39)[0x579257] ./Gale-1_5_0(_EntryPoint_Run_Class_VoidPtr+0x33)[0x5c6650] ./Gale-1_5_0(_AbstractContext_Step+0x103)[0x5cb569] ./Gale-1_5_0(_ContextEntryPoint_Run_Step+0x47)[0x5ceb20] ./Gale-1_5_0[0x5cae90] ./Gale-1_5_0(_AbstractContext_Execute_Hook+0x1b3)[0x5cb8b8] ./Gale-1_5_0(_EntryPoint_Run_VoidPtr+0x49)[0x5c72a3] ./Gale-1_5_0(_AbstractContext_Execute+0x5f)[0x5cb2c3] ./Gale-1_5_0[0x5c36b1] ./Gale-1_5_0[0x5d446c] ./Gale-1_5_0[0x4592b4] /lib64/libc.so.6(__libc_start_main+0xf4)[0x33b441d974] ./Gale-1_5_0(realloc+0x229)[0x458cb9] ======= Memory map: ======== 00400000-00c2e000 r-xp 00000000 00:18 175833140 /home/markb/Gale/Gale-1_5_0 00e2d000-00e3f000 rw-p 0082d000 00:18 175833140 /home/markb/Gale/Gale-1_5_0 00e3f000-00f79000 rw-p 00e3f000 00:00 0 1fe93000-abe9d000 rw-p 1fe93000 00:00 0 [heap] 33b3400000-33b341c000 r-xp 00000000 08:02 30671045 /lib64/ld-2.5.so 33b361b000-33b361c000 r--p 0001b000 08:02 30671045 /lib64/ld-2.5.so 33b361c000-33b361d000 rw-p 0001c000 08:02 30671045 /lib64/ld-2.5.so 33b4400000-33b454c000 r-xp 00000000 08:02 30671178 /lib64/libc-2.5.so 33b454c000-33b474c000 ---p 0014c000 08:02 30671178 /lib64/libc-2.5.so 33b474c000-33b4750000 r--p 0014c000 08:02 30671178 /lib64/libc-2.5.so 33b4750000-33b4751000 rw-p 00150000 08:02 30671178 /lib64/libc-2.5.so 33b4751000-33b4756000 rw-p 33b4751000 00:00 0 33b4800000-33b4882000 r-xp 00000000 08:02 30670883 /lib64/libm-2.5.so 33b4882000-33b4a81000 ---p 00082000 08:02 30670883 /lib64/libm-2.5.so 33b4a81000-33b4a82000 r--p 00081000 08:02 30670883 /lib64/libm-2.5.so 33b4a82000-33b4a83000 rw-p 00082000 08:02 30670883 /lib64/libm-2.5.so 33b4c00000-33b4c02000 r-xp 00000000 08:02 30671179 /lib64/libdl-2.5.so 33b4c02000-33b4e02000 ---p 00002000 08:02 30671179 /lib64/libdl-2.5.so 33b4e02000-33b4e03000 r--p 00002000 08:02 30671179 /lib64/libdl-2.5.so 33b4e03000-33b4e04000 rw-p 00003000 08:02 30671179 /lib64/libdl-2.5.so 33b5000000-33b5016000 r-xp 00000000 08:02 30671180 /lib64/libpthread-2.5.so 33b5016000-33b5215000 ---p 00016000 08:02 30671180 /lib64/libpthread-2.5.so 33b5215000-33b5216000 r--p 00015000 08:02 30671180 /lib64/libpthread-2.5.so 33b5216000-33b5217000 rw-p 00016000 08:02 30671180 /lib64/libpthread-2.5.so 33b5217000-33b521b000 rw-p 33b5217000 00:00 0 33b6000000-33b600d000 r-xp 00000000 08:02 30671196 /lib64/libgcc_s-4.1.2-20080825.so.1 33b600d000-33b620d000 ---p 0000d000 08:02 30671196 /lib64/libgcc_s-4.1.2-20080825.so.1 33b620d000-33b620e000 rw-p 0000d000 08:02 30671196 /lib64/libgcc_s-4.1.2-20080825.so.1 33b7400000-33b7407000 r-xp 00000000 08:02 30671198 /lib64/librt-2.5.so 33b7407000-33b7607000 ---p 00007000 08:02 30671198 /lib64/librt-2.5.so 33b7607000-33b7608000 r--p 00007000 08:02 30671198 /lib64/librt-2.5.so 33b7608000-33b7609000 rw-p 00008000 08:02 30671198 /lib64/librt-2.5.so 2b32e83a9000-2b32e83aa000 rw-p 2b32e83a9000 00:00 0 2b32e83d1000-2b32e83d5000 rw-p 2b32e83d1000 00:00 0 2b32e83d5000-2b32e83df000 r-xp 00000000 08:02 30670915 /lib64/libnss_files-2.5.so 2b32e83df000-2b32e85de000 ---p 0000a000 08:02 30670915 /lib64/libnss_files-2.5.so 2b32e85de000-2b32e85df000 r--p 00009000 08:02 30670915 /lib64/libnss_files-2.5.so 2b32e85df000-2b32e85e0000 rw-p 0000a000 08:02 30670915 /lib64/libnss_files-2.5.so 2b32e8751000-2b32e88c2000 rw-p 2b32e8751000 00:00 0 2b32e8a33000-2b338c6e8000 rw-p 2b32e8a33000 00:00 0 2b339394d000-2b33d7159000 rw-p 2b339394d000 00:00 0 2b33d8000000-2b33d8021000 rw-p 2b33d8000000 00:00 0 2b33d8021000-2b33dc000000 ---p 2b33d8021000 00:00 0 2b33f1cfd000-2b3410961000 rw-p 2b338f858000 00:00 0 7fff73ed9000-7fff73eee000 rw-p 7ffffffea000 00:00 0 [stack] ffffffffff600000-ffffffffffe00000 ---p 00000000 00:00 0 [vdso] Abort [markb@marsh ~/Gale]$ auto-logout Connection to marsh.geology.yale.edu closed. wedge:~ markbrandon$