/// <summary>
/// Load the current-controlled switch for AC analysis
/// </summary>
/// <param name="ckt">The circuit</param>
public override void AcLoad(Circuit ckt)
{
CurrentSwitchModel model = Model as CurrentSwitchModel;
double current_state;
double g_now;
var state = ckt.State;
var cstate = state.Complex;
// Get the current state
current_state = state.States[0][CSWstate];
g_now = current_state > 0.0 ? model.CSWonConduct : model.CSWoffConduct;
// Load the Y-matrix
cstate.Matrix[CSWposNode, CSWposNode] += g_now;
cstate.Matrix[CSWposNode, CSWnegNode] -= g_now;
cstate.Matrix[CSWnegNode, CSWposNode] -= g_now;
cstate.Matrix[CSWnegNode, CSWnegNode] += g_now;
}
/// <summary>
/// Load the current-controlled switch
/// </summary>
/// <param name="ckt">The circuit</param>
public override void Load(Circuit ckt)
{
CurrentSwitchModel model = Model as CurrentSwitchModel;
double g_now;
double i_ctrl;
double previous_state;
double current_state = 0.0;
var state = ckt.State;
var rstate = state.Real;
// decide the state of the switch
if (state.Init == CircuitState.InitFlags.InitFix || state.Init == CircuitState.InitFlags.InitJct)
{
if (CSWzero_state)
{
// Switch specified "on"
state.States[0][CSWstate] = 1.0;
current_state = 1.0;
}
else
{
// Switch specified "off"
state.States[0][CSWstate] = 0.0;
current_state = 0.0;
}
}
else if (state.UseSmallSignal)
{
previous_state = state.States[0][CSWstate];
current_state = previous_state;
}
else if (state.UseDC)
{
// No time-dependence, so use current state instead
previous_state = state.States[0][CSWstate];
i_ctrl = rstate.OldSolution[CSWcontBranch];
if (i_ctrl > (model.CSWthresh + model.CSWhyst))
{
current_state = 1.0;
}
else if (i_ctrl < (model.CSWthresh - model.CSWhyst))
{
current_state = 0.0;
}
else
{
current_state = previous_state;
}
// Store the current state
if (current_state == 0)
{
state.States[0][CSWstate] = 0.0;
}
else
{
state.States[0][CSWstate] = 1.0;
}
// Ensure one more iteration
if (current_state != previous_state)
{
state.IsCon = false;
}
}
else
{
// Get the previous state
previous_state = state.States[1][CSWstate];
i_ctrl = rstate.OldSolution[CSWcontBranch];
// Calculate the current state
if (i_ctrl > (model.CSWthresh + model.CSWhyst))
{
current_state = 1;
}
else if (i_ctrl < (model.CSWthresh - model.CSWhyst))
{
current_state = 0;
}
else
{
current_state = previous_state;
}
// Store the current state
if (current_state == 0)
{
state.States[0][CSWstate] = 0.0;
}
else
{
state.States[0][CSWstate] = 1.0;
}
}
// Get the current conduction
g_now = current_state != 0.0 ? (model.CSWonConduct) : (model.CSWoffConduct);
CSWcond = g_now;
// Load the Y-matrix
rstate.Matrix[CSWposNode, CSWposNode] += g_now;
rstate.Matrix[CSWposNode, CSWnegNode] -= g_now;
rstate.Matrix[CSWnegNode, CSWposNode] -= g_now;
rstate.Matrix[CSWnegNode, CSWnegNode] += g_now;
}
/// <summary> /// Set the model for the current-controlled switch /// </summary> public void SetModel(CurrentSwitchModel model) => Model = (ICircuitObject)model;
/// <summary> /// Set the model for the current-controlled switch /// </summary> public void SetModel(CurrentSwitchModel model) => Model = model;
