(Vessel vessel,
int childIndex,
double [] a,
double [] b,
double [] c)
{result = new Matrix (2 * atm. size () + 1, 1);
int parentIndex = vessel. index;. set (2 * parentIndex, 0, a [parentIndex]);. set (2 * parentIndex + 1, 0, b [parentIndex]);. set (2 * childIndex, 0, - a [childIndex]);. set (2 * childIndex + 1, 0, - b [childIndex]);. set (2 * atm. size (), 0, c [parentIndex] - c [childIndex]);
return result;
}
private Matrix resolveSingleFluxCondition ( int i, int j, int k, double [] a, double [] b, double [] c)
{result = new Matrix (2 * atm. size () + 1, 1);. set (2 * i + 1, 0, 1);. set (2 * j + 1, 0, - 1);. set (2 * k + 1, 0, - 1);
return result;
}
private Matrix resolveSinglePressureEndCondition ( int vesselIndex, double [] a, double [] b, double [] c, double flux0, double fluxFin, double t)
{result = new Matrix (2 * atm. size () + 1, 1);. set (2 * vesselIndex, 0, a [vesselIndex]);. set (2 * vesselIndex + 1, 0, b [vesselIndex]);
if (vesselIndex == 0)
{. set (2 * atm. size (), 0, c [vesselIndex]);
}
else
{. set (2 * atm. size (), 0, c [vesselIndex] + 70 * (1 - Math. exp (- (t * t)/100)));
}
<...
