er. resolveMatrix (); = leftResolver. getMatrix (leftMatrix); = leftResolver. getVector (leftMatrix);
for (j = 0; j
{
if (k2 == 0)
{[j] = V [j]. get (0, NumberOfSegmentation); [j] = V [j]. get (1, NumberOfSegmentation);
}
else
{[j] = U0 [j]. get (0, NumberOfSegmentation); [j] = U0 [j]. get (1, NumberOfSegmentation);
} [j] = vesselBeta [j]/(vesselAreas [j] * (Math. sqrt (An [j]) + Math. sqrt (vesselAreas [j]))); [j] = (rho * Qn [ j])/(2 * An [j] * An [j]); [j] = vesselBeta [j]/(Math. sqrt (An [j]) + Math. sqrt (vesselAreas [j]));
}
// resolve right matrix {8} rightResolver = new MatrixResolver (atm, a, b, c, < b align = "justify"> false , vesselAreas, flux1, tMin); rightMatrix = rightResolver. resolveMatrix (); = rightResolver. getMatrix (rightMatrix); = rightResolver. getVector (rightMatrix);
// experiment repeat marching was here. out . println ("k2" + k2);
// null matrixes Y_f, Z_f, beta, Z_0, Y_0, Ort, U_f. setMatrix (0, NumberOfVessels * 2 - 1, 0, NumberOfSegmentation, ZeroY_f); _f. setMatrix (0, NumberOfVessels * 2 - 1, 0, NumberOfSegmentation, ZeroY_f);. setMatrix (0, NumberOfVessels, 0, NumberOfSegmentation, ZeroY_f. getMatrix (0, NumberOfVessels, 0, NumberOfSegmentation));
for ( int i1 = 0; i1
{_0 [i1]. setMatrix (0, NumberOfVessels * 2 - 1, 0, NumberOfVessels - 1, Zero. getMatrix (0, NumberOfVessels * 2 - 1, 0, NumberOfVessels - 1)); _0 [i1]. setMatrix (0, NumberOfVessels * 2 - 1, 0, NumberOfVessels - 1, Zero. getMatrix (0, NumberOfVessels * 2 - 1, 0, NumberOfVessels - 1));
if (i1
{[i1]. setMatrix (0, NumberOfVessels, 0, NumberOfVessels, Zero
. getMatrix (0, NumberOfVessels, 0, NumberOfVessels));
}
}
for ( int i2 = 0; i2
