How to improve the computing speed of a THM model on OGS6?

Dear @joergbuchwald

Thanks so much for your help. I finally made configuration and build work on Windows shown in this post How to use intel mkl in cmake-gui under Windows? - #3 by Rui_Feng. But I would like to inquire how to switch off scaling for PardisoLU since I encountered this error when I used PardisoLU as my linear solver:

PS D:\Software\OGS6\test_speed\PardisoLU> ogs cold_water_real_3600in_3600no_3600in.prj
info: This is OpenGeoSys-6 version NO_VERSION.
info: OGS started on 2022-02-05 16:08:25+0800.
debug: Project configuration from file 'cold_water_real_3600in_3600no_3600in.prj' read.
debug: Reading multiple meshes.
debug: Reading mesh file 'fault_playing_mesh_linear.vtu'.
debug: Reading mesh file 'fault_doing_geometry_left_boundary.vtu'.
debug: Reading mesh file 'fault_doing_geometry_bottom_boundary.vtu'.
debug: Reading mesh file 'fault_doing_geometry_right_boundary.vtu'.
debug: Reading mesh file 'fault_doing_geometry_top_boundary.vtu'.
debug: Reading mesh file 'fault_playing_production.vtu'.
debug: Reading mesh file 'fault_playing_injection.vtu'.
debug: Reading curves configuration.
debug: Reading parameters:
info: FunctionParameter: sigma0
info: CurveScaledParameter: thermal_expansivity_increasing
debug: Using curve thermal_expansivity_curve
debug: Using parameter thermal_expansivity_curve_base
info: GroupBasedParameter: thermal_expansivity_curve_base
debug: Using group_id_property MaterialIDs
info: ConstantParameter: E
debug: Using value 5e+10 for constant parameter.
info: ConstantParameter: nu
debug: Using value 0.3 for constant parameter.
info: ConstantParameter: T0
debug: Using value 423.15 for constant parameter.
info: ConstantParameter: temperature0
debug: Using value 423.15 for constant parameter.
info: ConstantParameter: temperature_left_boundary
debug: Using value 423.15 for constant parameter.
info: ConstantParameter: temperature_bottom_boundary
debug: Using value 423.15 for constant parameter.
info: ConstantParameter: temperature_right_boundary
debug: Using value 423.15 for constant parameter.
info: ConstantParameter: temperature_top_boundary
debug: Using value 423.15 for constant parameter.
info: FunctionParameter: P0
info: FunctionParameter: P_left_boundary
info: FunctionParameter: P_bottom_boundary
info: FunctionParameter: P_right_boundary
info: FunctionParameter: P_top_boundary
info: ConstantParameter: displacement0
debug: Using following values for the constant parameter:
debug:  0
debug:  0
info: ConstantParameter: dirichlet0
debug: Using following values for the constant parameter:
debug:  0
info: CurveScaledParameter: p_source_term_out_well
debug: Using curve p_source_term_out_well_curve
debug: Using parameter p_source_term_out_well_curve_base
info: ConstantParameter: p_source_term_out_well_curve_base
debug: Using value 1 for constant parameter.
info: CurveScaledParameter: p_source_term_in_well
debug: Using curve p_source_term_in_well_curve
debug: Using parameter p_source_term_in_well_curve_base
info: ConstantParameter: p_source_term_in_well_curve_base
debug: Using value 1 for constant parameter.
info: CurveScaledParameter: T_source_term_in_well
debug: Using curve T_source_term_in_well_curve
debug: Using parameter T_source_term_in_well_curve_base
info: ConstantParameter: T_source_term_in_well_curve_base
debug: Using value 1 for constant parameter.
debug: Found parameter `thermal_expansivity_curve_base'.
debug: Found parameter `p_source_term_out_well_curve_base'.
debug: Found parameter `p_source_term_in_well_curve_base'.
debug: Found parameter `T_source_term_in_well_curve_base'.
debug: Parse process variables:
debug: Found parameter `temperature0'.
debug: Constructing process variable temperature
debug: Found mesh 'fault_doing_geometry_left_boundary' with id 1.
debug: Found mesh 'fault_doing_geometry_bottom_boundary' with id 2.
debug: Found mesh 'fault_doing_geometry_right_boundary' with id 3.
debug: Found mesh 'fault_doing_geometry_top_boundary' with id 4.
debug: Found mesh 'fault_playing_injection' with id 6.
debug: Found parameter `P0'.
debug: Constructing process variable pressure
debug: Found mesh 'fault_doing_geometry_left_boundary' with id 1.
debug: Found mesh 'fault_doing_geometry_bottom_boundary' with id 2.
debug: Found mesh 'fault_doing_geometry_right_boundary' with id 3.
debug: Found mesh 'fault_doing_geometry_top_boundary' with id 4.
debug: Found mesh 'fault_playing_production' with id 5.
debug: Found mesh 'fault_playing_injection' with id 6.
debug: Found parameter `displacement0'.
debug: Constructing process variable displacement
debug: Found mesh 'fault_doing_geometry_left_boundary' with id 1.
debug: Found mesh 'fault_doing_geometry_bottom_boundary' with id 2.
debug: Found mesh 'fault_doing_geometry_right_boundary' with id 3.
info: No source terms for process variable 'displacement' found.
debug: Reading media:
debug: Create Constant property specific_heat_capacity.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property density.
debug: Create Constant property thermal_expansivity.
debug: Create Constant property viscosity.
debug: Create Constant property storage.
debug: Create Constant property density.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Parameter property thermal_expansivity.
debug: Found parameter `thermal_expansivity_increasing'.
debug: Create Constant property thermal_longitudinal_dispersivity.
debug: Create Constant property thermal_transversal_dispersivity.
debug: Create Constant property permeability.
debug: Create Constant property porosity.
debug: Create Constant property biot_coefficient.
debug: Create effective thermal_conductivity property from porosity mixing thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property density.
debug: Create Constant property thermal_expansivity.
debug: Create Constant property viscosity.
debug: Create Constant property storage.
debug: Create Constant property density.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Parameter property thermal_expansivity.
debug: Found parameter `thermal_expansivity_increasing'.
debug: Create Constant property thermal_longitudinal_dispersivity.
debug: Create Constant property thermal_transversal_dispersivity.
debug: Create Constant property permeability.
debug: Create Constant property porosity.
debug: Create Constant property biot_coefficient.
debug: Create effective thermal_conductivity property from porosity mixing thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property density.
debug: Create Constant property thermal_expansivity.
debug: Create Constant property viscosity.
debug: Create Constant property storage.
debug: Create Constant property density.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Parameter property thermal_expansivity.
debug: Found parameter `thermal_expansivity_increasing'.
debug: Create Constant property thermal_longitudinal_dispersivity.
debug: Create Constant property thermal_transversal_dispersivity.
debug: Create Constant property permeability.
debug: Create Constant property porosity.
debug: Create Constant property biot_coefficient.
debug: Create effective thermal_conductivity property from porosity mixing thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property density.
debug: Create Constant property thermal_expansivity.
debug: Create Constant property viscosity.
debug: Create Constant property storage.
debug: Create Constant property density.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Parameter property thermal_expansivity.
debug: Found parameter `thermal_expansivity_increasing'.
debug: Create Constant property thermal_longitudinal_dispersivity.
debug: Create Constant property thermal_transversal_dispersivity.
debug: Create Constant property permeability.
debug: Create Constant property porosity.
debug: Create Constant property biot_coefficient.
debug: Create effective thermal_conductivity property from porosity mixing thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property density.
debug: Create Constant property thermal_expansivity.
debug: Create Constant property viscosity.
debug: Create Constant property storage.
debug: Create Constant property density.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Parameter property thermal_expansivity.
debug: Found parameter `thermal_expansivity_increasing'.
debug: Create Constant property thermal_longitudinal_dispersivity.
debug: Create Constant property thermal_transversal_dispersivity.
debug: Create Constant property permeability.
debug: Create Constant property porosity.
debug: Create Constant property biot_coefficient.
debug: Create effective thermal_conductivity property from porosity mixing thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property density.
debug: Create Constant property thermal_expansivity.
debug: Create Constant property viscosity.
debug: Create Constant property storage.
debug: Create Constant property density.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Parameter property thermal_expansivity.
debug: Found parameter `thermal_expansivity_increasing'.
debug: Create Constant property thermal_longitudinal_dispersivity.
debug: Create Constant property thermal_transversal_dispersivity.
debug: Create Constant property permeability.
debug: Create Constant property porosity.
debug: Create Constant property biot_coefficient.
debug: Create effective thermal_conductivity property from porosity mixing thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property density.
debug: Create Constant property thermal_expansivity.
debug: Create Constant property viscosity.
debug: Create Constant property storage.
debug: Create Constant property density.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Parameter property thermal_expansivity.
debug: Found parameter `thermal_expansivity_increasing'.
debug: Create Constant property thermal_longitudinal_dispersivity.
debug: Create Constant property thermal_transversal_dispersivity.
debug: Create Constant property permeability.
debug: Create Constant property porosity.
debug: Create Constant property biot_coefficient.
debug: Create effective thermal_conductivity property from porosity mixing thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property density.
debug: Create Constant property thermal_expansivity.
debug: Create Constant property viscosity.
debug: Create Constant property storage.
debug: Create Constant property density.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Parameter property thermal_expansivity.
debug: Found parameter `thermal_expansivity_increasing'.
debug: Create Constant property thermal_longitudinal_dispersivity.
debug: Create Constant property thermal_transversal_dispersivity.
debug: Create Constant property permeability.
debug: Create Constant property porosity.
debug: Create Constant property biot_coefficient.
debug: Create effective thermal_conductivity property from porosity mixing thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property density.
debug: Create Constant property thermal_expansivity.
debug: Create Constant property viscosity.
debug: Create Constant property storage.
debug: Create Constant property density.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Parameter property thermal_expansivity.
debug: Found parameter `thermal_expansivity_increasing'.
debug: Create Constant property thermal_longitudinal_dispersivity.
debug: Create Constant property thermal_transversal_dispersivity.
debug: Create Constant property permeability.
debug: Create Constant property porosity.
debug: Create Constant property biot_coefficient.
debug: Create effective thermal_conductivity property from porosity mixing thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property density.
debug: Create Constant property thermal_expansivity.
debug: Create Constant property viscosity.
debug: Create Constant property storage.
debug: Create Constant property density.
debug: Create Constant property thermal_conductivity.
debug: Create Constant property specific_heat_capacity.
debug: Create Parameter property thermal_expansivity.
debug: Found parameter `thermal_expansivity_increasing'.
debug: Create Constant property thermal_longitudinal_dispersivity.
debug: Create Constant property thermal_transversal_dispersivity.
debug: Create Constant property permeability.
debug: Create Constant property porosity.
debug: Create Constant property biot_coefficient.
debug: Create effective thermal_conductivity property from porosity mixing thermal_conductivity.
debug: Reading linear solver configuration.
debug: Reading processes:
debug: Create ThermoHydroMechanicsProcess.
debug: Found process variable 'temperature' for config tag <temperature>.
debug: Found process variable 'pressure' for config tag <pressure>.
debug: Found process variable 'displacement' for config tag <displacement>.
debug: Associate displacement with process variable 'displacement'.
debug: Associate pressure with process variable 'pressure'.
debug: Associate temperature with process variable 'temperature'.
debug: Create LinearElasticIsotropic material
debug: Found parameter `E'.
debug: Use 'E' as youngs_modulus parameter.
debug: Found parameter `nu'.
debug: Use 'nu' as poissons_ratio parameter.
debug: Found 1 constitutive relations.
debug: Found parameter `sigma0'.
debug: Reading non-linear solver configuration.
debug: Reading time loop configuration.
debug: Parse output configuration:
debug: adding output variable `displacement'
debug: adding output variable `pressure'
debug: adding output variable `temperature'
debug: adding output variable `sigma'
debug: adding output variable `epsilon'
info: Initialize processes.
debug: Initialize process.
debug: Construct dof mappings.
debug: Compute sparsity pattern
debug: Initialize the extrapolator
debug: Create local assemblers.
debug: Calling local assembler builder for all mesh elements.
debug: Initialize boundary conditions.
debug: Constructing DirichletBoundaryCondition from config.
debug: Using parameter temperature_left_boundary
debug: Found parameter `temperature_left_boundary'.
debug: Found 21 nodes for Dirichlet BCs for the variable 0 and component 0
debug: Construct reduced local to global index map.
debug: Constructing DirichletBoundaryCondition from config.
debug: Using parameter temperature_bottom_boundary
debug: Found parameter `temperature_bottom_boundary'.
debug: Found 21 nodes for Dirichlet BCs for the variable 0 and component 0
debug: Construct reduced local to global index map.
debug: Constructing DirichletBoundaryCondition from config.
debug: Using parameter temperature_right_boundary
debug: Found parameter `temperature_right_boundary'.
debug: Found 21 nodes for Dirichlet BCs for the variable 0 and component 0
debug: Construct reduced local to global index map.
debug: Constructing DirichletBoundaryCondition from config.
debug: Using parameter temperature_top_boundary
debug: Found parameter `temperature_top_boundary'.
debug: Found 21 nodes for Dirichlet BCs for the variable 0 and component 0
debug: Construct reduced local to global index map.
debug: Constructing DirichletBoundaryCondition from config.
debug: Using parameter P_left_boundary
debug: Found parameter `P_left_boundary'.
debug: Found 21 nodes for Dirichlet BCs for the variable 1 and component 0
debug: Construct reduced local to global index map.
debug: Constructing DirichletBoundaryCondition from config.
debug: Using parameter P_bottom_boundary
debug: Found parameter `P_bottom_boundary'.
debug: Found 21 nodes for Dirichlet BCs for the variable 1 and component 0
debug: Construct reduced local to global index map.
debug: Constructing DirichletBoundaryCondition from config.
debug: Using parameter P_right_boundary
debug: Found parameter `P_right_boundary'.
debug: Found 21 nodes for Dirichlet BCs for the variable 1 and component 0
debug: Construct reduced local to global index map.
debug: Constructing DirichletBoundaryCondition from config.
debug: Using parameter P_top_boundary
debug: Found parameter `P_top_boundary'.
debug: Found 21 nodes for Dirichlet BCs for the variable 1 and component 0
debug: Construct reduced local to global index map.
debug: Constructing DirichletBoundaryCondition from config.
debug: Using parameter dirichlet0
debug: Found parameter `dirichlet0'.
debug: Found 21 nodes for Dirichlet BCs for the variable 2 and component 0
debug: Construct reduced local to global index map.
debug: Constructing DirichletBoundaryCondition from config.
debug: Using parameter dirichlet0
debug: Found parameter `dirichlet0'.
debug: Found 21 nodes for Dirichlet BCs for the variable 2 and component 1
debug: Construct reduced local to global index map.
debug: Constructing DirichletBoundaryCondition from config.
debug: Using parameter dirichlet0
debug: Found parameter `dirichlet0'.
debug: Found 21 nodes for Dirichlet BCs for the variable 2 and component 0
debug: Construct reduced local to global index map.
debug: Found 312 nodes for source term at mesh 'fault_playing_injection' for the variable 0 and component 0
debug: Construct reduced local to global index map.
debug: Constructing VolumetricSourceTerm from config.
debug: Found parameter `T_source_term_in_well'.
debug: Using 'T_source_term_in_well' as volumetric source term parameter.
debug: Create local assemblers.
debug: Create local assemblers.
debug: Calling local assembler builder for all mesh elements.
debug: Found 312 nodes for source term at mesh 'fault_playing_production' for the variable 1 and component 0
debug: Construct reduced local to global index map.
debug: Constructing VolumetricSourceTerm from config.
debug: Found parameter `p_source_term_out_well'.
debug: Using 'p_source_term_out_well' as volumetric source term parameter.
debug: Create local assemblers.
debug: Create local assemblers.
debug: Calling local assembler builder for all mesh elements.
debug: Found 312 nodes for source term at mesh 'fault_playing_injection' for the variable 1 and component 0
debug: Construct reduced local to global index map.
debug: Constructing VolumetricSourceTerm from config.
debug: Found parameter `p_source_term_in_well'.
debug: Using 'p_source_term_in_well' as volumetric source term parameter.
debug: Create local assemblers.
debug: Create local assemblers.
debug: Calling local assembler builder for all mesh elements.
info: Solve processes.
debug: Set the initial condition of variable temperature of process 0.
debug: Set the initial condition of variable pressure of process 0.
debug: Set the initial condition of variable displacement of process 0.
debug: PreTimestep ThermoHydroMechanicsProcess.
debug: Compute the secondary variables for ThermoHydroMechanicsProcess.
debug: Process output data.
debug:   process variable temperature
debug:   process variable pressure
debug:   process variable displacement
debug:   secondary variable epsilon
debug: Computing new singular value decomposition
debug:   secondary variable sigma
debug: Writing output to 'Initialization_ts_0_t_0.000000.vtu'.
info: [time] Output of timestep 0 took 6.23828 s.
info: === Time stepping at step #1 and time 100 with step size 100
info: Calculate non-equilibrium initial residuum.
debug: PreTimestep ThermoHydroMechanicsProcess.
debug: Assemble the Jacobian of ThermoHydroMechanics for the monolithic scheme.
debug: Assemble VolumetricSourceTerm.
debug: Assemble VolumetricSourceTerm.
debug: Assemble VolumetricSourceTerm.
info: [time] Assembly took 147.413 s.
info: [time] Applying Dirichlet BCs took 0.40803 s.
info: ------------------------------------------------------------------
info: *** Eigen solver computation
info: -> scale
info: -> solve with Eigen direct linear solver PardisoLU
error: Failed during Eigen linear solver initialization
info: ------------------------------------------------------------------
info: [time] Linear solver took 29.9208 s.
error: Newton: The linear solver failed.
info: [time] Solving process #0 took 177.785 s in time step #1
error: The nonlinear solver failed in time step #1 at t = 100 s for process #0.
debug: Process output data.
debug:   process variable temperature
debug:   process variable pressure
debug:   process variable displacement
debug:   secondary variable epsilon
debug:   secondary variable sigma
debug: Writing output to 'Initialization_ts_1_t_100.000000.vtu'.
info: [time] Output of timestep 1 took 6.31346 s.
critical: D:/Software/OGS6/ogs6_4_1_source_code_10/ogs-master/ProcessLib/TimeLoop.cpp:694 ProcessLib::TimeLoop::solveUncoupledEquationSystems()
PS D:\Software\OGS6\test_speed\PardisoLU>

Note of OS information:
ogs source code version: 6.4.1 (latest), from ogs / ogs ยท GitLab (by zip),
Operating system / hardware information: Windows 11 (I am not using WSL2, just Windows 11),Intel Core i7-1185G7 3.00GHz, 16 GB RAM),
CMake-gui version: 3.22.2
Intel MKL version: w_BaseKit_p_2022.1.2.154_offline.exe
my prj files:
test_PardisoLU.zip (1.0 MB)

Btw, I cannot find any BLAS variables in CMake-gui.

Best,
Rui