Function ImmersX::Vascularization::solve_reference(const NetworkData&, const std::vector<TerminalPair>&, const std::vector<TerminalTriplet>&, const ReferenceConditions&, const ReferenceParameters&, MPI_Comm)

Function Documentation

ReferenceResult ImmersX::Vascularization::solve_reference(const NetworkData &network, const std::vector<TerminalPair> &pairs, const std::vector<TerminalTriplet> &triplets, const ReferenceConditions &conditions, const ReferenceParameters &parameters, MPI_Comm communicator)

Solve two- and three-terminal territories together.

Each triplet has one arterial and one portal supply plus one venous drain, all with the same tissue demand. Positive star-arm resistances determine the common microvascular pressure. Its algebraic row is eliminated locally, leaving the same distributed native graph operator and no artificial 1D vessel. Pair/triplet ids and terminal membership must be unique collectively.