fn connection_per_phase(
info: &DomainInfo,
mass_balance: &mut PfrGlobalMassBalance,
phase_node: &[FluxType],
phase: PhaseCM,
) -> Result<RawDataFlux, CMError>fn connection_per_phase(
info: &DomainInfo,
mass_balance: &mut PfrGlobalMassBalance,
phase_node: &[FluxType],
phase: PhaseCM,
) -> Result<RawDataFlux, CMError>