blob: 5e5388d0de5a488720cd59ba95b7f5d796283667 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
|
/***************************************************************************
* Copyright (C) 2003-2004 by David Saxton *
* [email protected] *
* *
* This program is free software; you can redistribute it and/or modify *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation; either version 2 of the License, or *
* (at your option) any later version. *
***************************************************************************/
#include "currentsignal.h"
#include "element.h"
#include "matrix.h"
CurrentSignal::CurrentSignal( const double delta, const double current )
: Reactive::Reactive(delta)
{
m_current = current;
m_oldCurrent = m_newCurrent = 0.;
m_numCNodes = 2;
}
CurrentSignal::~CurrentSignal()
{
}
void CurrentSignal::setCurrent( const double i )
{
if ( m_oldCurrent != m_newCurrent ) add_initial_dc();
m_newCurrent *= i/m_current; // Instead of calling step again, we can just "adjust" what the current should be
m_current = i;
add_initial_dc();
}
void CurrentSignal::add_map()
{
// We don't need a map for current signal :-)
}
void CurrentSignal::add_initial_dc()
{
if ( !b_status )
return;
if ( m_newCurrent == m_oldCurrent )
return;
b_i( 0 ) -= m_newCurrent-m_oldCurrent;
b_i( 1 ) += m_newCurrent-m_oldCurrent;
m_oldCurrent = m_newCurrent;
}
void CurrentSignal::updateCurrents()
{
m_cnodeI[1] = m_newCurrent;
m_cnodeI[0] = -m_newCurrent;
}
void CurrentSignal::time_step()
{
add_initial_dc(); // Make sure our old and new are synced
m_newCurrent = m_current*advance();
add_initial_dc();
}
|
