public CouplLinesInFreqRange(Params currentParams, double l, double fmin, double fmax, int nf, double z1In, double z2In, double z1Out, double z2Out)
{
L = l * Math.Pow(10, -3);
C = 0.2998;
Z1In = z1In;
Z2In = z2In;
Z1Out = z1Out;
Z2Out = z2Out;
Fmin = fmin;
Fmax = fmax;
Nf = nf;
S11 = Vector <double> .Build.Dense(Nf);
S12 = Vector <double> .Build.Dense(Nf);
S13 = Vector <double> .Build.Dense(Nf);
S14 = Vector <double> .Build.Dense(Nf);
S22 = Vector <double> .Build.Dense(Nf);
S24 = Vector <double> .Build.Dense(Nf);
S23 = Vector <double> .Build.Dense(Nf);
S33 = Vector <double> .Build.Dense(Nf);
S34 = Vector <double> .Build.Dense(Nf);
S44 = Vector <double> .Build.Dense(Nf);
S21 = Vector <double> .Build.Dense(Nf);
S31 = Vector <double> .Build.Dense(Nf);
S32 = Vector <double> .Build.Dense(Nf);
S41 = Vector <double> .Build.Dense(Nf);
S42 = Vector <double> .Build.Dense(Nf);
S43 = Vector <double> .Build.Dense(Nf);
F11 = Vector <double> .Build.Dense(Nf);
F12 = Vector <double> .Build.Dense(Nf);
F13 = Vector <double> .Build.Dense(Nf);
F14 = Vector <double> .Build.Dense(Nf);
F22 = Vector <double> .Build.Dense(Nf);
F24 = Vector <double> .Build.Dense(Nf);
F23 = Vector <double> .Build.Dense(Nf);
F33 = Vector <double> .Build.Dense(Nf);
F34 = Vector <double> .Build.Dense(Nf);
F44 = Vector <double> .Build.Dense(Nf);
F21 = Vector <double> .Build.Dense(Nf);
F31 = Vector <double> .Build.Dense(Nf);
F32 = Vector <double> .Build.Dense(Nf);
F41 = Vector <double> .Build.Dense(Nf);
F42 = Vector <double> .Build.Dense(Nf);
F43 = Vector <double> .Build.Dense(Nf);
E4 = Matrix <double> .Build.DenseIdentity(4);
Zl = Matrix <double> .Build.DenseIdentity(4);
Zl[0, 0] = Math.Sqrt(Z1In);
Zl[1, 1] = Math.Sqrt(Z2In);
Zl[2, 2] = Math.Sqrt(Z1Out);
Zl[3, 3] = Math.Sqrt(Z2Out);
ZlComplex = Zl.ToComplex();
E4Complex = E4.ToComplex();
Y = new Matrix <Complex> [Nf];
Fi = Vector <double> .Build.Dense(Nf, i => ((Fmax - Fmin) / (Nf - 1)) *(i + 1) + Fmin);
Wi = Fi.Multiply(2 * Math.PI);
Um = Matrix <double> .Build.Dense(2, 2);
Um[0, 0] = 1;
Um[0, 1] = 1;
Um[1, 0] = currentParams.Rc;
Um[1, 1] = currentParams.Rp;
Im = Matrix <double> .Build.Dense(2, 2);
Im[0, 0] = Math.Pow(currentParams.Zc1, -1);
Im[0, 1] = Math.Pow(currentParams.Zp1, -1);
Im[1, 0] = -Math.Pow(currentParams.Zc1 * currentParams.Rp, -1);
Im[1, 1] = -Math.Pow(currentParams.Zp1 * currentParams.Rc, -1);
ElectricalLengthForC = Vector <Complex> .Build.Dense(Nf, i => _j *Wi[i] *(Math.Sqrt(currentParams.Erc) / C) *L);
ElectricalLengthForP = Vector <Complex> .Build.Dense(Nf, i => _j *Wi[i] *(Math.Sqrt(currentParams.Erp) / C) *L);
var tet = new Vector <Complex> [Nf];
for (int i = 0; i < Nf; i++)
{
tet[i] = Vector <Complex> .Build.Dense(new[] { ElectricalLengthForC[i], ElectricalLengthForP[i] });
}
var coIntermediate = new Vector <Complex> [Nf];
var co = new Vector <Complex> [Nf];
for (int i = 0; i < Nf; i++)
{
var coIntermediate1 = Vector <Complex> .Build.Dense(2);
var coIntermediate2 = Vector <Complex> .Build.Dense(2);
coIntermediate[i] = tet[i];
coIntermediate[i].PointwiseTanh(coIntermediate1);
coIntermediate1.DivideByThis(1, coIntermediate2);
co[i] = coIntermediate2;
}
var scIntermediate = new Vector <Complex> [Nf];
var sc = new Vector <Complex> [Nf];
for (int i = 0; i < Nf; i++)
{
var scIntermediate1 = Vector <Complex> .Build.Dense(2);
var scIntermediate2 = Vector <Complex> .Build.Dense(2);
scIntermediate[i] = tet[i];
scIntermediate[i].PointwiseSinh(scIntermediate1);
scIntermediate1.DivideByThis(1, scIntermediate2);
sc[i] = scIntermediate2;
}
var diagCo = new Matrix <Complex> [Nf];
for (int i = 0; i < Nf; i++)
{
diagCo[i] = Matrix <Complex> .Build.DenseOfDiagonalVector(co[i]);
}
var diagSc = new Matrix <Complex> [Nf];
for (int i = 0; i < Nf; i++)
{
diagSc[i] = Matrix <Complex> .Build.DenseOfDiagonalVector(sc[i]);
}
var yaa = new Matrix <Complex> [Nf];
var imIntermediate = Matrix <Complex> .Build.Dense(2, 2);
var umIntermediate = Matrix <Complex> .Build.Dense(2, 2);
imIntermediate = Im.ToComplex();
umIntermediate = Um.Inverse().ToComplex();
var diagCoIntermediate = diagCo;
for (int i = 0; i < Nf; i++)
{
yaa[i] = imIntermediate.Multiply(diagCoIntermediate[i]).Multiply(umIntermediate);
}
var yab = new Matrix <Complex> [Nf];
var diagScIntermediate = diagSc;
for (int i = 0; i < Nf; i++)
{
yab[i] = -imIntermediate.Multiply(diagScIntermediate[i]).Multiply(umIntermediate);
}
var yaayab1 = new Matrix <Complex> [Nf];
var yaayab2 = new Matrix <Complex> [Nf];
var yo = new Matrix <Complex> [Nf];
for (int i = 0; i < Nf; i++)
{
yaayab1[i] = yaa[i].Append(yab[i]);
yaayab2[i] = yab[i].Append(yaa[i]);
yo[i] = yaayab1[i].Stack(yaayab2[i]);
}
if (Z1Out > Math.Pow(10, 9) && Z2Out > Math.Pow(10, 9))
{
var additionMatrix = Matrix <double> .Build.Dense(4, 4);
additionMatrix[3, 3] = Math.Pow(10, 10);
additionMatrix[2, 2] = Math.Pow(10, 10);
additionMatrix[3, 2] = -Math.Pow(10, 10);
additionMatrix[2, 3] = -Math.Pow(10, 10);
var complexAddMatrix = additionMatrix.ToComplex();
for (int i = 0; i < Nf; i++)
{
yo[i] = yo[i] + complexAddMatrix;
}
}
for (int i = 0; i < Nf; i++)
{
Y[i] = ZlComplex.Multiply(yo[i]).Multiply(ZlComplex);
}
var SS = new Matrix <Complex> [Nf];
for (int i = 0; i < Nf; i++)
{
SS[i] = 2 * (E4Complex + Y[i]).Inverse() - E4Complex;
}
var SS11 = Vector <Complex> .Build.Dense(Nf);
var SS12 = Vector <Complex> .Build.Dense(Nf);
var SS13 = Vector <Complex> .Build.Dense(Nf);
var SS14 = Vector <Complex> .Build.Dense(Nf);
var SS22 = Vector <Complex> .Build.Dense(Nf);
var SS24 = Vector <Complex> .Build.Dense(Nf);
var SS23 = Vector <Complex> .Build.Dense(Nf);
var SS33 = Vector <Complex> .Build.Dense(Nf);
var SS34 = Vector <Complex> .Build.Dense(Nf);
var SS44 = Vector <Complex> .Build.Dense(Nf);
var SS21 = Vector <Complex> .Build.Dense(Nf);
var SS31 = Vector <Complex> .Build.Dense(Nf);
var SS32 = Vector <Complex> .Build.Dense(Nf);
var SS41 = Vector <Complex> .Build.Dense(Nf);
var SS42 = Vector <Complex> .Build.Dense(Nf);
var SS43 = Vector <Complex> .Build.Dense(Nf);
for (int i = 0; i < Nf; i++)
{
SS11[i] = SS[i][0, 0];
SS12[i] = SS[i][0, 1];
SS13[i] = SS[i][0, 2];
SS14[i] = SS[i][0, 3];
SS22[i] = SS[i][1, 1];
SS24[i] = SS[i][1, 3];
SS23[i] = SS[i][1, 2];
SS33[i] = SS[i][2, 2];
SS34[i] = SS[i][2, 3];
SS44[i] = SS[i][3, 3];
SS21[i] = SS[i][1, 0];
SS31[i] = SS[i][2, 0];
SS32[i] = SS[i][2, 1];
SS41[i] = SS[i][3, 0];
SS42[i] = SS[i][3, 1];
SS43[i] = SS[i][3, 2];
}
for (int i = 0; i < Nf; i++)
{
S11[i] = 20 * Math.Log10(SS11.PointwiseAbs()[i].Real);
S12[i] = 20 * Math.Log10(SS12.PointwiseAbs()[i].Real);
S13[i] = 20 * Math.Log10(SS13.PointwiseAbs()[i].Real);
S14[i] = 20 * Math.Log10(SS14.PointwiseAbs()[i].Real);
S22[i] = 20 * Math.Log10(SS22.PointwiseAbs()[i].Real);
S24[i] = 20 * Math.Log10(SS24.PointwiseAbs()[i].Real);
S23[i] = 20 * Math.Log10(SS23.PointwiseAbs()[i].Real);
S33[i] = 20 * Math.Log10(SS33.PointwiseAbs()[i].Real);
S34[i] = 20 * Math.Log10(SS34.PointwiseAbs()[i].Real);
S44[i] = 20 * Math.Log10(SS44.PointwiseAbs()[i].Real);
S21[i] = 20 * Math.Log10(SS21.PointwiseAbs()[i].Real);
S31[i] = 20 * Math.Log10(SS31.PointwiseAbs()[i].Real);
S32[i] = 20 * Math.Log10(SS32.PointwiseAbs()[i].Real);
S41[i] = 20 * Math.Log10(SS41.PointwiseAbs()[i].Real);
S42[i] = 20 * Math.Log10(SS42.PointwiseAbs()[i].Real);
S43[i] = 20 * Math.Log10(SS43.PointwiseAbs()[i].Real);
F11[i] = 180 / Math.PI * SS11[i].Phase;
F12[i] = 180 / Math.PI * SS12[i].Phase;
F13[i] = 180 / Math.PI * SS13[i].Phase;
F14[i] = 180 / Math.PI * SS14[i].Phase;
F22[i] = 180 / Math.PI * SS22[i].Phase;
F24[i] = 180 / Math.PI * SS24[i].Phase;
F23[i] = 180 / Math.PI * SS23[i].Phase;
F33[i] = 180 / Math.PI * SS33[i].Phase;
F34[i] = 180 / Math.PI * SS34[i].Phase;
F44[i] = 180 / Math.PI * SS44[i].Phase;
F21[i] = 180 / Math.PI * SS21[i].Phase;
F31[i] = 180 / Math.PI * SS31[i].Phase;
F32[i] = 180 / Math.PI * SS32[i].Phase;
F41[i] = 180 / Math.PI * SS41[i].Phase;
F42[i] = 180 / Math.PI * SS42[i].Phase;
F43[i] = 180 / Math.PI * SS43[i].Phase;
}
}