kicad/eeschema/sim/sim_plot_panel.cpp

651 lines
16 KiB
C++
Raw Normal View History

2016-08-11 12:41:07 +00:00
/*
* This program source code file is part of KiCad, a free EDA CAD application.
*
* Copyright (C) 2016 CERN
* @author Tomasz Wlostowski <tomasz.wlostowski@cern.ch>
* @author Maciej Suminski <maciej.suminski@cern.ch>
2016-08-11 12:41:07 +00:00
*
* 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 3
2016-08-11 12:41:07 +00:00
* of the License, or (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program; if not, you may find one here:
* https://www.gnu.org/licenses/gpl-3.0.html
* or you may search the http://www.gnu.org website for the version 3 license,
2016-08-11 12:41:07 +00:00
* or you may write to the Free Software Foundation, Inc.,
* 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA
*/
2016-08-11 12:41:01 +00:00
#include "sim_plot_panel.h"
#include <algorithm>
2016-08-11 12:41:14 +00:00
#include <limits>
2016-08-11 12:42:17 +00:00
static wxString formatFloat( double x, int nDigits )
{
wxString rv, fmt;
2016-08-11 12:42:17 +00:00
if( nDigits )
{
2016-08-11 12:42:17 +00:00
fmt = wxT( "%.0Nf" );
fmt[3] = '0' + nDigits;
2016-08-11 12:42:17 +00:00
}
else
{
fmt = wxT( "%.0f" );
}
2016-08-11 12:42:17 +00:00
rv.Printf( fmt, x );
return rv;
}
2016-08-11 12:42:17 +00:00
static void getSISuffix( double x, const wxString& unit, int& power, wxString& suffix )
{
const int n_powers = 11;
2016-08-11 12:42:17 +00:00
const struct
{
double exponent;
char suffix;
2016-08-11 12:42:17 +00:00
} powers[] =
{
{ -18, 'a' },
{ -15, 'f' },
{ -12, 'p' },
{ -9, 'n' },
{ -6, 'u' },
{ -3, 'm' },
{ 0, 0 },
{ 3, 'k' },
{ 6, 'M' },
{ 9, 'G' },
{ 12, 'T' },
{ 14, 'P' }
};
power = 0;
suffix = unit;
2016-08-11 12:42:17 +00:00
if( x == 0.0 )
return;
2016-08-11 12:42:17 +00:00
for( int i = 0; i < n_powers - 1; i++ )
{
2016-08-11 12:42:17 +00:00
double r_cur = pow( 10, powers[i].exponent );
2016-08-11 12:42:17 +00:00
if( fabs( x ) >= r_cur && fabs( x ) < r_cur * 1000.0 )
{
power = powers[i].exponent;
2016-08-11 12:42:17 +00:00
if( powers[i].suffix )
suffix = wxString( powers[i].suffix ) + unit;
else
suffix = unit;
2016-08-11 12:42:17 +00:00
return;
}
}
}
2016-08-11 12:42:17 +00:00
static int countDecimalDigits( double x, int maxDigits )
{
2016-08-11 12:42:17 +00:00
int64_t k = (int)( ( x - floor( x ) ) * pow( 10.0, (double) maxDigits ) );
int n = 0;
2016-08-11 12:42:17 +00:00
while( k && ( ( k % 10LL ) == 0LL || ( k % 10LL ) == 9LL ) )
{
k /= 10LL;
}
n = 0;
2016-08-11 12:42:17 +00:00
while( k != 0LL )
{
n++;
k /= 10LL;
}
return n;
}
2016-08-11 12:42:17 +00:00
static void formatSILabels( mpScaleBase* scale, const wxString& aUnit, int nDigits )
{
double maxVis = scale->AbsVisibleMaxValue();
wxString suffix;
int power, digits = 0;
2016-08-11 12:42:17 +00:00
getSISuffix( maxVis, aUnit, power, suffix );
2016-08-11 12:42:17 +00:00
double sf = pow( 10.0, power );
2016-08-11 12:42:17 +00:00
for( auto &l : scale->TickLabels() )
{
int k = countDecimalDigits( l.pos / sf, nDigits );
2016-08-11 12:42:17 +00:00
digits = std::max( digits, k );
}
2016-08-11 12:42:17 +00:00
for( auto &l : scale->TickLabels() )
{
l.label = formatFloat ( l.pos / sf, digits ) + suffix;
l.visible = true;
}
}
2016-08-11 12:42:17 +00:00
class FREQUENCY_LOG_SCALE : public mpScaleXLog
{
public:
2016-08-11 12:42:17 +00:00
FREQUENCY_LOG_SCALE( wxString name, int flags ) :
mpScaleXLog( name, flags ) {};
void formatLabels() override
{
2016-08-11 12:42:17 +00:00
const wxString unit = wxT( "Hz" );
wxString suffix;
int power;
2016-08-11 12:42:17 +00:00
for( auto &l : TickLabels() )
{
2016-08-11 12:42:17 +00:00
getSISuffix( l.pos, unit, power, suffix );
double sf = pow( 10.0, power );
int k = countDecimalDigits( l.pos / sf, 3 );
2016-08-11 12:42:17 +00:00
l.label = formatFloat( l.pos / sf, k ) + suffix;
l.visible = true;
}
}
};
2016-08-11 12:42:17 +00:00
class FREQUENCY_LIN_SCALE : public mpScaleX
{
public:
2016-08-11 12:42:17 +00:00
FREQUENCY_LIN_SCALE( wxString name, int flags ) :
mpScaleX( name, flags, false , 0 ) {};
void formatLabels() override
{
2016-08-11 12:42:17 +00:00
formatSILabels( this, wxT( "Hz" ), 3 );
}
};
class TIME_SCALE : public mpScaleX
{
public:
2016-08-11 12:42:17 +00:00
TIME_SCALE( wxString name, int flags ) :
mpScaleX( name, flags, false, 0 ) {};
void formatLabels() override
{
2016-08-11 12:42:17 +00:00
formatSILabels( this, wxT( "s" ), 3 );
}
};
2016-08-11 12:42:17 +00:00
2016-08-11 12:41:53 +00:00
class VOLTAGE_SCALE_X : public mpScaleX
{
public:
2016-08-11 12:42:17 +00:00
VOLTAGE_SCALE_X( wxString name, int flags ) :
mpScaleX( name, flags, false, 0 ) {};
2016-08-11 12:41:53 +00:00
void formatLabels() override
2016-08-11 12:41:53 +00:00
{
2016-08-11 12:42:17 +00:00
formatSILabels( this, wxT( "V" ), 3 );
2016-08-11 12:41:53 +00:00
}
};
2016-08-11 12:42:17 +00:00
class GAIN_SCALE : public mpScaleY
{
public:
GAIN_SCALE( wxString name, int flags ) :
2016-08-11 12:42:17 +00:00
mpScaleY( name, flags, false ) {};
void formatLabels() override
{
2016-08-11 12:42:17 +00:00
formatSILabels( this, wxT( "dB" ), 3 );
}
};
2016-08-11 12:42:17 +00:00
class PHASE_SCALE : public mpScaleY
{
public:
2016-08-11 12:42:17 +00:00
PHASE_SCALE( wxString name, int flags ) :
mpScaleY( name, flags, false ) {};
void formatLabels() override
{
2016-08-11 12:42:17 +00:00
formatSILabels( this, wxT( "\u00B0" ), 3 ); // degree sign
}
};
2016-08-11 12:42:17 +00:00
2016-08-11 12:41:53 +00:00
class VOLTAGE_SCALE_Y : public mpScaleY
{
public:
2016-08-11 12:42:17 +00:00
VOLTAGE_SCALE_Y( wxString name, int flags ) :
mpScaleY( name, flags, false ) {};
void formatLabels() override
{
2016-08-11 12:42:17 +00:00
formatSILabels( this, wxT( "V" ), 3 );
}
};
2016-08-11 12:42:17 +00:00
class CURRENT_SCALE : public mpScaleY
{
public:
2016-08-11 12:42:17 +00:00
CURRENT_SCALE( wxString name, int flags ) :
mpScaleY( name, flags, false ) {};
void formatLabels() override
{
2016-08-11 12:42:17 +00:00
formatSILabels( this, wxT( "A" ), 3 );
}
};
2016-08-11 12:41:22 +00:00
void CURSOR::Plot( wxDC& aDC, mpWindow& aWindow )
{
if( !m_window )
m_window = &aWindow;
2016-08-11 12:41:22 +00:00
if( !m_visible )
return;
const auto& dataX = m_trace->GetDataX();
const auto& dataY = m_trace->GetDataY();
if( dataX.size() <= 1 )
return;
if( m_updateRequired )
2016-08-11 12:41:22 +00:00
{
2016-08-11 12:41:55 +00:00
m_coords.x = m_trace->s2x( aWindow.p2x( m_dim.x ) );
2016-08-11 12:41:22 +00:00
// Find the closest point coordinates
auto maxXIt = std::upper_bound( dataX.begin(), dataX.end(), m_coords.x );
int maxIdx = maxXIt - dataX.begin();
int minIdx = maxIdx - 1;
// Out of bounds checks
if( minIdx < 0 )
{
minIdx = 0;
maxIdx = 1;
m_coords.x = dataX[0];
}
else if( maxIdx >= (int) dataX.size() )
{
maxIdx = dataX.size() - 1;
minIdx = maxIdx - 1;
m_coords.x = dataX[maxIdx];
}
const double leftX = dataX[minIdx];
const double rightX = dataX[maxIdx];
const double leftY = dataY[minIdx];
const double rightY = dataY[maxIdx];
2016-08-11 12:41:55 +00:00
// Linear interpolation
2016-08-11 12:41:22 +00:00
m_coords.y = leftY + ( rightY - leftY ) / ( rightX - leftX ) * ( m_coords.x - leftX );
m_updateRequired = false;
// Notify the parent window about the changes
wxQueueEvent( aWindow.GetParent(), new wxCommandEvent( EVT_SIM_CURSOR_UPDATE ) );
2016-08-11 12:41:22 +00:00
}
else
2016-08-11 12:42:14 +00:00
{
m_updateRef = true;
}
if( m_updateRef )
2016-08-11 12:41:22 +00:00
{
UpdateReference();
2016-08-11 12:42:14 +00:00
m_updateRef = false;
2016-08-11 12:41:22 +00:00
}
// Line length in horizontal and vertical dimensions
2016-08-11 12:41:55 +00:00
const wxPoint cursorPos( aWindow.x2p( m_trace->x2s( m_coords.x ) ),
aWindow.y2p( m_trace->y2s( m_coords.y ) ) );
2016-08-11 12:41:55 +00:00
wxCoord leftPx = m_drawOutsideMargins ? 0 : aWindow.GetMarginLeft();
wxCoord rightPx = m_drawOutsideMargins ? aWindow.GetScrX() : aWindow.GetScrX() - aWindow.GetMarginRight();
wxCoord topPx = m_drawOutsideMargins ? 0 : aWindow.GetMarginTop();
wxCoord bottomPx = m_drawOutsideMargins ? aWindow.GetScrY() : aWindow.GetScrY() - aWindow.GetMarginBottom();
aDC.SetPen( wxPen( *wxWHITE, 1, m_continuous ? wxPENSTYLE_SOLID : wxPENSTYLE_LONG_DASH ) );
2016-08-11 12:41:55 +00:00
if( topPx < cursorPos.y && cursorPos.y < bottomPx )
aDC.DrawLine( leftPx, cursorPos.y, rightPx, cursorPos.y );
if( leftPx < cursorPos.x && cursorPos.x < rightPx )
aDC.DrawLine( cursorPos.x, topPx, cursorPos.x, bottomPx );
2016-08-11 12:41:22 +00:00
}
2016-08-11 12:41:55 +00:00
bool CURSOR::Inside( wxPoint& aPoint )
{
if( !m_window )
return false;
return ( std::abs( (double) aPoint.x - m_window->x2p( m_trace->x2s( m_coords.x ) ) ) <= DRAG_MARGIN )
|| ( std::abs( (double) aPoint.y - m_window->y2p( m_trace->y2s( m_coords.y ) ) ) <= DRAG_MARGIN );
2016-08-11 12:41:55 +00:00
}
void CURSOR::UpdateReference()
{
if( !m_window )
return;
m_reference.x = m_window->x2p( m_trace->x2s( m_coords.x ) );
m_reference.y = m_window->y2p( m_trace->y2s( m_coords.y ) );
}
SIM_PLOT_PANEL::SIM_PLOT_PANEL( SIM_TYPE aType, wxWindow* parent, wxWindowID id, const wxPoint& pos,
const wxSize& size, long style, const wxString& name )
: mpWindow( parent, id, pos, size, style ), m_colorIdx( 0 ),
m_axis_x( nullptr ), m_axis_y1( nullptr ), m_axis_y2( nullptr ), m_type( aType )
2016-08-11 12:41:01 +00:00
{
LimitView( true );
2016-08-11 12:41:53 +00:00
SetMargins( 50, 80, 50, 80 );
2016-08-11 12:41:53 +00:00
wxColour grey( 96, 96, 96 );
SetColourTheme( *wxBLACK, *wxWHITE, grey );
EnableDoubleBuffer( true );
UpdateAll();
2016-08-11 12:41:01 +00:00
switch( m_type )
{
case ST_AC:
2016-08-11 12:42:01 +00:00
m_axis_x = new FREQUENCY_LOG_SCALE( wxT( "Frequency" ), mpALIGN_BOTTOM );
m_axis_y1 = new GAIN_SCALE( wxT( "Gain" ), mpALIGN_LEFT );
m_axis_y2 = new PHASE_SCALE( wxT( "Phase" ), mpALIGN_RIGHT );
2016-08-11 12:42:17 +00:00
m_axis_y2->SetMasterScale( m_axis_y1 );
break;
case ST_DC:
m_axis_x = new VOLTAGE_SCALE_X( wxT( "Voltage (sweeped)" ), mpALIGN_BOTTOM );
m_axis_y1 = new VOLTAGE_SCALE_Y( wxT( "Voltage (measured)" ), mpALIGN_LEFT );
break;
case ST_NOISE:
m_axis_x = new FREQUENCY_LOG_SCALE( wxT( "Frequency" ), mpALIGN_BOTTOM );
m_axis_y1 = new mpScaleY( wxT( "noise [(V or A)^2/Hz]" ), mpALIGN_LEFT );
break;
case ST_TRANSIENT:
m_axis_x = new TIME_SCALE( wxT( "Time" ), mpALIGN_BOTTOM );
2016-08-11 12:41:53 +00:00
m_axis_y1 = new VOLTAGE_SCALE_Y( wxT( "Voltage" ), mpALIGN_LEFT );
m_axis_y2 = new CURRENT_SCALE( wxT( "Current" ), mpALIGN_RIGHT );
2016-08-11 12:42:17 +00:00
m_axis_y2->SetMasterScale( m_axis_y1 );
break;
default:
// suppress warnings
break;
}
if( m_axis_x )
{
m_axis_x->SetTicks( false );
m_axis_x->SetNameAlign ( mpALIGN_BOTTOM );
AddLayer( m_axis_x );
}
if( m_axis_y1 )
{
m_axis_y1->SetTicks( false );
m_axis_y1->SetNameAlign ( mpALIGN_LEFT );
AddLayer( m_axis_y1 );
}
if( m_axis_y2 )
{
m_axis_y2->SetTicks( false );
m_axis_y2->SetNameAlign ( mpALIGN_RIGHT );
AddLayer( m_axis_y2 );
}
m_legend = new mpInfoLegend( wxRect( 0, 40, 200, 40 ), wxTRANSPARENT_BRUSH );
m_legend->SetVisible( false );
AddLayer( m_legend );
m_topLevel.push_back( m_legend );
2016-08-11 12:42:17 +00:00
SetColourTheme( *wxBLACK, *wxWHITE, grey );
2016-08-11 12:42:17 +00:00
EnableDoubleBuffer( true );
UpdateAll();
}
2016-08-11 12:41:14 +00:00
SIM_PLOT_PANEL::~SIM_PLOT_PANEL()
{
// ~mpWindow destroys all the added layers, so there is no need to destroy m_traces contents
2016-08-11 12:41:14 +00:00
}
bool SIM_PLOT_PANEL::IsPlottable( SIM_TYPE aSimType )
{
switch( aSimType )
{
case ST_AC:
case ST_DC:
case ST_TRANSIENT:
return true;
default:
return false;
}
}
2016-08-11 12:41:52 +00:00
bool SIM_PLOT_PANEL::AddTrace( const wxString& aName, int aPoints,
const double* aX, const double* aY, SIM_PLOT_TYPE aFlags )
2016-08-11 12:41:01 +00:00
{
TRACE* trace = NULL;
// Find previous entry, if there is one
2016-08-11 12:41:52 +00:00
auto prev = m_traces.find( aName );
2016-08-11 12:41:24 +00:00
bool addedNewEntry = ( prev == m_traces.end() );
2016-08-11 12:41:01 +00:00
2016-08-11 12:41:24 +00:00
if( addedNewEntry )
{
2016-08-11 12:42:17 +00:00
if( m_type == ST_TRANSIENT )
{
bool hasVoltageTraces = false;
for( auto tr : m_traces )
{
if( !( tr.second->GetFlags() & SPT_CURRENT ) )
{
hasVoltageTraces = true;
break;
}
}
2016-08-11 12:42:17 +00:00
if( !hasVoltageTraces )
m_axis_y2->SetMasterScale( nullptr );
else
m_axis_y2->SetMasterScale( m_axis_y1 );
}
// New entry
trace = new TRACE( aName );
trace->SetTraceColour( generateColor() );
trace->SetPen( wxPen( trace->GetTraceColour(), 2, wxPENSTYLE_SOLID ) );
m_traces[aName] = trace;
// It is a trick to keep legend & coords always on the top
for( mpLayer* l : m_topLevel )
DelLayer( l );
AddLayer( (mpLayer*) trace );
for( mpLayer* l : m_topLevel )
AddLayer( l );
}
else
{
trace = prev->second;
}
2016-08-11 12:41:14 +00:00
std::vector<double> tmp( aY, aY + aPoints );
if( m_type == ST_AC )
{
2016-08-11 12:41:52 +00:00
if( aFlags & SPT_AC_PHASE )
{
2016-08-11 12:42:17 +00:00
for( int i = 0; i < aPoints; i++ )
tmp[i] = tmp[i] * 180.0 / M_PI; // convert to degrees
}
else
{
2016-08-11 12:42:17 +00:00
for( int i = 0; i < aPoints; i++ )
tmp[i] = 20 * log( tmp[i] ) / log( 10.0 ); // convert to dB
}
}
trace->SetData( std::vector<double>( aX, aX + aPoints ), tmp );
if( aFlags & SPT_AC_PHASE || aFlags & SPT_CURRENT )
trace->SetScale( m_axis_x, m_axis_y2 );
else
trace->SetScale( m_axis_x, m_axis_y1 );
trace->SetFlags( aFlags );
UpdateAll();
2016-08-11 12:41:24 +00:00
return addedNewEntry;
}
bool SIM_PLOT_PANEL::DeleteTrace( const wxString& aName )
{
2016-08-11 12:41:52 +00:00
auto it = m_traces.find( aName );
2016-08-11 12:41:24 +00:00
if( it != m_traces.end() )
2016-08-11 12:41:24 +00:00
{
TRACE* trace = it->second;
m_traces.erase( it );
if( CURSOR* cursor = trace->GetCursor() )
DelLayer( cursor, true );
DelLayer( trace, true, true );
2016-08-11 12:42:04 +00:00
ResetScales();
2016-08-11 12:42:03 +00:00
2016-08-11 12:41:24 +00:00
return true;
}
return false;
2016-08-11 12:41:01 +00:00
}
2016-08-11 12:41:07 +00:00
2016-08-11 12:41:24 +00:00
void SIM_PLOT_PANEL::DeleteAllTraces()
2016-08-11 12:41:01 +00:00
{
2016-08-11 12:41:24 +00:00
for( auto& t : m_traces )
2016-08-11 12:41:12 +00:00
{
2016-08-11 12:41:52 +00:00
DeleteTrace( t.first );
2016-08-11 12:41:12 +00:00
}
2016-08-11 12:42:17 +00:00
m_colorIdx = 0;
m_traces.clear();
}
2016-08-11 12:41:12 +00:00
2016-08-11 12:41:14 +00:00
bool SIM_PLOT_PANEL::HasCursorEnabled( const wxString& aName ) const
{
TRACE* t = GetTrace( aName );
return t ? t->HasCursor() : false;
}
void SIM_PLOT_PANEL::EnableCursor( const wxString& aName, bool aEnable )
{
TRACE* t = GetTrace( aName );
if( t == nullptr || t->HasCursor() == aEnable )
return;
if( aEnable )
{
CURSOR* c = new CURSOR( t );
2016-08-11 12:42:14 +00:00
int plotCenter = GetMarginLeft() + ( GetXScreen() - GetMarginLeft() - GetMarginRight() ) / 2;
c->SetX( plotCenter );
t->SetCursor( c );
AddLayer( c );
}
else
{
CURSOR* c = t->GetCursor();
t->SetCursor( NULL );
DelLayer( c, true );
}
// Notify the parent window about the changes
wxQueueEvent( GetParent(), new wxCommandEvent( EVT_SIM_CURSOR_UPDATE ) );
}
2016-08-11 12:42:04 +00:00
void SIM_PLOT_PANEL::ResetScales()
{
if( m_axis_x )
m_axis_x->ResetDataRange();
if( m_axis_y1 )
m_axis_y1->ResetDataRange();
if( m_axis_y2 )
m_axis_y2->ResetDataRange();
for( auto t : m_traces )
t.second->UpdateScales();
}
wxColour SIM_PLOT_PANEL::generateColor()
{
/// @todo have a look at:
/// http://stanford.edu/~mwaskom/software/seaborn/tutorial/color_palettes.html
/// https://github.com/Gnuplotting/gnuplot-palettes
2016-08-11 12:41:12 +00:00
const unsigned long colors[] = { 0x0000ff, 0x00ff00, 0xff0000, 0x00ffff, 0xff00ff, 0xffff000, 0xffffff };
2016-08-11 12:41:14 +00:00
//const unsigned long colors[] = { 0xe3cea6, 0xb4781f, 0x8adfb2, 0x2ca033, 0x999afb, 0x1c1ae3, 0x6fbffd, 0x007fff, 0xd6b2ca, 0x9a3d6a };
2016-08-11 12:41:14 +00:00
// hls
//const unsigned long colors[] = { 0x0f1689, 0x0f7289, 0x35890f, 0x0f8945, 0x89260f, 0x890f53, 0x89820f, 0x630f89 };
2016-08-11 12:41:14 +00:00
// pastels, good for dark background
//const unsigned long colors[] = { 0x2fd8fe, 0x628dfa, 0x53d8a6, 0xa5c266, 0xb3b3b3, 0x94c3e4, 0xca9f8d, 0xac680e };
2016-08-11 12:41:14 +00:00
const unsigned int colorCount = sizeof(colors) / sizeof(unsigned long);
2016-08-11 12:41:14 +00:00
/// @todo generate shades to avoid repeating colors
return wxColour( colors[m_colorIdx++ % colorCount] );
2016-08-11 12:41:14 +00:00
}
wxDEFINE_EVENT( EVT_SIM_CURSOR_UPDATE, wxCommandEvent );