Parallel geothermal numerical model with faults and multi-branch wells