petsc-error deactivating temperature in subdomain in ThermoRichardsMechanics
Description
Solution with petsc returns error in first solution call when deactivating heat transport equation in ThermoRichardsMechanics. Deactivating other balance equations does not lead to that error.
Steps to Reproduce the Problem
- Deactivate Heat transport for Test Tests/Data/ThermoRichardsMechanics/TaskCDECOVALEX2023/Decovalex-0.prj in Material Group 0 by adding following snippet in
process_variables->process_variable[temperature]
<deactivated_subdomains>
<deactivated_subdomain>
<time_interval>
<start>0</start><end>1e99</end></time_interval>
<material_ids>0</material_ids>
</deactivated_subdomain>
</deactivated_subdomains>
- execute with
mpiexec -np 3
Expected behavior: Problem should run. Although physics is changed, the thermal isolation of the bentonite materials should not yield failure in first step.
Actual behavior:
[1]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[1]PETSC ERROR: Object is in wrong state
[1]PETSC ERROR: Matrix is missing diagonal entry in the zeroed row 0 ...
Specifications
- Version: 6.5.1-96-g23dd6660
- Platform: Linux - openSUSE
CC @michaelpitz