---

elements.cc

Go to the documentation of this file.

// elements.cc
// SuperMix version 1.0  C++ source file
//
// Copyright (c) 1999 California Institute of Technology.
// All rights reserved.
//
// Redistribution and use in source and binary forms for noncommercial
// purposes are permitted provided that the above copyright notice and
// this paragraph are duplicated in all such forms and that any
// documentation and other materials related to such distribution and
// use acknowledge that the software was developed by California
// Institute of Technology. Redistribution and/or use in source or
// binary forms is not permitted for any commercial purpose. Use of
// this software does not include a permitted use of the Institute's
// name or trademark for any purpose.
//
// DISCLAIMER:
// THIS SOFTWARE AND/OR RELATED MATERIALS ARE PROVIDED "AS-IS" WITHOUT
// WARRANTY OF ANY KIND INCLUDING ANY WARRANTIES OF PERFORMANCE OR
// MERCHANTABILITY OR FITNESS FOR A PARTICULAR USE OR PURPOSE (AS SET
// FORTH IN UCC 23212-2313) OR FOR ANY PURPOSE WHATSOEVER, FOR THE
// LICENSED PRODUCT, HOWEVER USED.  IN NO EVENT SHALL CALTECH/JPL BE
// LIABLE FOR ANY DAMAGES AND/OR COSTS, INCLUDING BUT NOT LIMITED TO
// INCIDENTAL OR CONSEQUENTIAL DAMAGES OF ANY KIND, INCLUDING ECONOMIC
// DAMAGE OR INJURY TO PROPERTY AND LOST PROFITS, REGARDLESS OF
// WHETHER CALTECH/JPL SHALL BE ADVISED, HAVE REASON TO KNOW, OR IN
// FACT SHALL KNOW OF THE POSSIBILITY.  THE USER BEARS ALL RISK
// RELATING TO QUALITY AND PERFORMANCE OF THE SOFTWARE AND/OR RELATED
// MATERIALS.
//
// Change history:
// 11/16/99: changes to support new noise calculations
// 10/29/99: Ensured data.set_znorm() called in recalc()'s
// 10/22/99: Updated for parameter device::Z0
// 11/11/98: Changed table access to new syntax
// 9/15/98:  Fixed a nasty bug: wrong argument order in
//           passive noise calculations.
// 8/27/98:  changed branch code: improved recalc, etc.
// 8/24/98:  Added additional constructors; made recalc private
//           changed some void return values.
// 8/21/98:  Added 1-ports zterm and yterm
// 8/19/98:  Added admittance rep to spimp; many improvements
// 7/16/98:  changed branch to allow branches = 2
// 12/30/97: fixed spimp to handle Z = 0 case
// 12/29/97: modified by FR for new sdata (with znorm)
// 12/19/97: modified by FR for matmath classes

#include "elements.h"
#include "units.h"
#include "error.h"

void spimp::calcZ(complex Z)
{
  data.set_znorm(device::Z0);   // added 12/29/97

  if(is_series) {
    if (Z != 0.0) {
      complex zo = 2 * device::Z0;  // common coefficient      
      data.S[1][1] = data.S[2][2] = Z/(zo + Z);
      data.S[2][1] = data.S[1][2] = zo/(zo + Z);
    }
    else /* Z == 0.0 */ {
      data.S[1][1] = data.S[2][2] = 0.0;
      data.S[2][1] = data.S[1][2] = 1.0;
    }
  }
  else /* parallel */ {
    if (Z != 0.0) {
      complex z = 2 * Z;  // common coefficient
      data.S[1][1] = data.S[2][2] = -device::Z0/(z + device::Z0);
      data.S[2][1] = data.S[1][2] = z/(z + device::Z0);
    }
    else /* Z == 0.0 */ {
      data.S[1][1] = data.S[2][2] = -1.0;
      data.S[2][1] = data.S[1][2] = 0.0;
    }
  }
}

void spimp::calcY(complex Y)
{
  data.set_znorm(device::Z0);

  if(is_series) {
    if (Y != 0.0) {
      complex zoy = 2 * device::Z0 * Y;  // common coefficient      
      data.S[1][1] = data.S[2][2] = 1/(1 + zoy);
      data.S[2][1] = data.S[1][2] = zoy/(1 + zoy);
    }
    else /* Y == 0.0 */ {
      data.S[1][1] = data.S[2][2] = 1.0;
      data.S[2][1] = data.S[1][2] = 0.0;
    }
  }
  else /* parallel */ {
    if (Y != 0.0) {
      complex zoy = device::Z0 * Y;  // common coefficient      
      data.S[1][1] = data.S[2][2] = -zoy/(2 + zoy);
      data.S[2][1] = data.S[1][2] = 2/(2 + zoy);
    }
    else /* Y == 0.0 */ {
      data.S[1][1] = data.S[2][2] = 0.0;
      data.S[2][1] = data.S[1][2] = 1.0;
    }
  }
}

// **************************************************************

branch::branch(int b) : nport(b), branches(b)
{ 
  if(b < 2)
    error::fatal("Invalid number of branches in branch constructor.");
  info.noise = info.active = info.source = false;
  calc();
}

branch & branch::set_branches(int b)
{
  if(b < 2)
    error::fatal("Branches must have at least two ports.");
  branches = b;
  calc();
  return *this;
}

void branch::calc()
{
  data.set_znorm(0.0);   // since S is independent of device::Z0
  data.resize(branches);

  double diagonal = (2.0 - branches) / branches;
  double off_diagonal = 2.0 / branches;

  data.S.fill(off_diagonal);

  for(int i=1; i<=branches; i++)
    data.S[i][i] = diagonal;

  // Assume that data.C is already a zero matrix.
}

// **************************************************************

series_tee::series_tee() : nport(3)
{
  data.set_znorm(0.0);   // since S is independent of device::Z0

  data.S[1][1] = 1;
  data.S[2][2] = 1;
  data.S[3][3] = 1;

  data.S[2][1] = -2;
  data.S[1][2] = -2;

  data.S[3][1] = 2;
  data.S[3][2] = 2;
  data.S[1][3] = 2;
  data.S[2][3] = 2;

  data.S /= 3.0;
  info.noise = info.active = info.source = false;
}

// **************************************************************

void zterm::recalc()
{
  recalc_S();
  data.passive_noise(f, Temp);
}

void zterm::recalc_S()
{
  data.set_znorm(device::Z0);
  data.S[1][1] = (Z == 0.0) ? -1.0 : (Z - double(device::Z0))/(Z + double(device::Z0));
}

// **************************************************************

void yterm::recalc()
{
  recalc_S();
  data.passive_noise(f, Temp);
}

void yterm::recalc_S()
{
  data.set_znorm(device::Z0);
  complex y = double(device::Z0) * Y;
  data.S[1][1] = (y == 0.0) ? 1.0 : (1 - y)/(1 + y);
}

---

Please direct comments and corrections to supermix@submm.caltech.edu
Go to the supermix home page
Generated by 1.2.7