/* * Crummy way to Britishise. */ void ise(void) { struct suftab *tab; for (tab = suftab; tab->suf; tab++) { /* Assume that suffix will contain 'z' if a1 or d1 do */ if (strchr(tab->suf, 'z')) { tab->suf = estrdup(tab->suf); ztos(tab->suf); if (strchr(tab->d1, 'z')) { tab->d1 = estrdup(tab->d1); ztos(tab->d1); } if (strchr(tab->a1, 'z')) { tab->a1 = estrdup(tab->a1); ztos(tab->a1); } } } }
void msmbend::calcSP (nr_double_t frequency) { setMatrixS (ztos (calcMatrixZ (frequency))); }
matvec ztos (matvec z, nr_complex_t z0) { return ztos (z, vector (z.getCols (), z0)); }
// Convert impedance matrix vector scattering parameter matrix vector. matvec ztos (matvec z, vector z0) { assert (z.getCols () == z.getRows () && z.getCols () == z0.getSize ()); matvec res (z.getSize (), z.getCols (), z.getRows ()); for (int i = 0; i < z.getSize (); i++) res.set (ztos (z.get (i), z0), i); return res; }