/* * This program source code file is part of KiCad, a free EDA CAD application. * * Copyright (C) 2016-2018 CERN * Copyright (C) 2018 KiCad Developers, see AUTHORS.txt for contributors. * * @author Tomasz Wlostowski * @author Maciej Suminski * * 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 * 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, * or you may write to the Free Software Foundation, Inc., * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA */ #include "ngspice.h" #include "spice_reporter.h" #include // LOCALE_IO #include #include #include #include using namespace std; static const wxChar* const traceNgspice = wxT( "KICAD_NGSPICE" ); NGSPICE::NGSPICE() : m_error( false ) { init_dll(); } NGSPICE::~NGSPICE() { } void NGSPICE::Init() { Command( "reset" ); } vector NGSPICE::GetPlot( const string& aName, int aMaxLen ) { LOCALE_IO c_locale; // ngspice works correctly only with C locale vector data; vector_info* vi = m_ngGet_Vec_Info( (char*) aName.c_str() ); if( vi ) { int length = aMaxLen < 0 ? vi->v_length : std::min( aMaxLen, vi->v_length ); data.reserve( length ); if( vi->v_realdata ) { for( int i = 0; i < length; i++ ) data.push_back( COMPLEX( vi->v_realdata[i], 0.0 ) ); } else if( vi->v_compdata ) { for( int i = 0; i < length; i++ ) data.push_back( COMPLEX( vi->v_compdata[i].cx_real, vi->v_compdata[i].cx_imag ) ); } } return data; } vector NGSPICE::GetRealPlot( const string& aName, int aMaxLen ) { LOCALE_IO c_locale; // ngspice works correctly only with C locale vector data; vector_info* vi = m_ngGet_Vec_Info( (char*) aName.c_str() ); if( vi ) { int length = aMaxLen < 0 ? vi->v_length : std::min( aMaxLen, vi->v_length ); data.reserve( length ); if( vi->v_realdata ) { for( int i = 0; i < length; i++ ) { data.push_back( vi->v_realdata[i] ); } } else if( vi->v_compdata ) { for( int i = 0; i < length; i++ ) { assert( vi->v_compdata[i].cx_imag == 0.0 ); data.push_back( vi->v_compdata[i].cx_real ); } } } return data; } vector NGSPICE::GetImagPlot( const string& aName, int aMaxLen ) { LOCALE_IO c_locale; // ngspice works correctly only with C locale vector data; vector_info* vi = m_ngGet_Vec_Info( (char*) aName.c_str() ); if( vi ) { int length = aMaxLen < 0 ? vi->v_length : std::min( aMaxLen, vi->v_length ); data.reserve( length ); if( vi->v_compdata ) { for( int i = 0; i < length; i++ ) { data.push_back( vi->v_compdata[i].cx_imag ); } } } return data; } vector NGSPICE::GetMagPlot( const string& aName, int aMaxLen ) { LOCALE_IO c_locale; // ngspice works correctly only with C locale vector data; vector_info* vi = m_ngGet_Vec_Info( (char*) aName.c_str() ); if( vi ) { int length = aMaxLen < 0 ? vi->v_length : std::min( aMaxLen, vi->v_length ); data.reserve( length ); if( vi->v_realdata ) { for( int i = 0; i < length; i++ ) data.push_back( vi->v_realdata[i] ); } else if( vi->v_compdata ) { for( int i = 0; i < length; i++ ) data.push_back( hypot( vi->v_compdata[i].cx_real, vi->v_compdata[i].cx_imag ) ); } } return data; } vector NGSPICE::GetPhasePlot( const string& aName, int aMaxLen ) { LOCALE_IO c_locale; // ngspice works correctly only with C locale vector data; vector_info* vi = m_ngGet_Vec_Info( (char*) aName.c_str() ); if( vi ) { int length = aMaxLen < 0 ? vi->v_length : std::min( aMaxLen, vi->v_length ); data.reserve( length ); if( vi->v_realdata ) { for( int i = 0; i < length; i++ ) data.push_back( 0.0 ); // well, that's life } else if( vi->v_compdata ) { for( int i = 0; i < length; i++ ) data.push_back( atan2( vi->v_compdata[i].cx_imag, vi->v_compdata[i].cx_real ) ); } } return data; } bool NGSPICE::LoadNetlist( const string& aNetlist ) { LOCALE_IO c_locale; // ngspice works correctly only with C locale vector lines; stringstream ss( aNetlist ); m_netlist = ""; while( !ss.eof() ) { char line[1024]; ss.getline( line, sizeof( line ) ); lines.push_back( strdup( line ) ); m_netlist += std::string( line ) + std::string( "\n" ); } lines.push_back( nullptr ); // sentinel, as requested in ngSpice_Circ description m_ngSpice_Circ( lines.data() ); for( auto line : lines ) free( line ); return true; } bool NGSPICE::Run() { LOCALE_IO c_locale; // ngspice works correctly only with C locale return Command( "bg_run" ); // bg_* commands execute in a separate thread } bool NGSPICE::Stop() { LOCALE_IO c_locale; // ngspice works correctly only with C locale return Command( "bg_halt" ); // bg_* commands execute in a separate thread } bool NGSPICE::IsRunning() { LOCALE_IO c_locale; // ngspice works correctly only with C locale return m_ngSpice_Running(); } bool NGSPICE::Command( const string& aCmd ) { LOCALE_IO c_locale; // ngspice works correctly only with C locale validate(); m_ngSpice_Command( (char*) aCmd.c_str() ); return true; } string NGSPICE::GetXAxis( SIM_TYPE aType ) const { switch( aType ) { case ST_AC: case ST_NOISE: return string( "frequency" ); break; case ST_DC: return string( "v-sweep" ); break; case ST_TRANSIENT: return string( "time" ); break; default: break; } return string( "" ); } void NGSPICE::init_dll() { if( m_initialized ) return; LOCALE_IO c_locale; // ngspice works correctly only with C locale const wxStandardPaths& stdPaths = wxStandardPaths::Get(); if( m_dll.IsLoaded() ) // enable force reload m_dll.Unload(); // Extra effort to find libngspice wxFileName dllFile( "", NGSPICE_DLL_FILE ); #if defined(__WINDOWS__) const vector dllPaths = { "", "/mingw64/bin", "/mingw32/bin" }; #elif defined(__WXMAC__) const vector dllPaths = { GetOSXKicadUserDataDir() + "/PlugIns/ngspice", GetOSXKicadMachineDataDir() + "/PlugIns/ngspice", // when running kicad.app stdPaths.GetPluginsDir() + "/sim", // when running eeschema.app wxFileName( stdPaths.GetExecutablePath() ).GetPath() + "/../../../../../Contents/PlugIns/sim" }; #else // Unix systems const vector dllPaths = { "/usr/local/lib" }; #endif #if defined(__WINDOWS__) || (__WXMAC__) for( const auto& path : dllPaths ) { dllFile.SetPath( path ); wxLogTrace( traceNgspice, "libngspice search path: %s", dllFile.GetFullPath() ); m_dll.Load( dllFile.GetFullPath(), wxDL_VERBATIM | wxDL_QUIET | wxDL_NOW ); if( m_dll.IsLoaded() ) { wxLogTrace( traceNgspice, "libngspice path found in: %s", dllFile.GetFullPath() ); break; } } if( !m_dll.IsLoaded() ) // try also the system libraries m_dll.Load( wxDynamicLibrary::CanonicalizeName( "ngspice" ) ); #else // here: __LINUX__ // First, try the system libraries m_dll.Load( NGSPICE_DLL_FILE, wxDL_VERBATIM | wxDL_QUIET | wxDL_NOW ); // If failed, try some other paths: if( !m_dll.IsLoaded() ) { for( const auto& path : dllPaths ) { dllFile.SetPath( path ); wxLogTrace( traceNgspice, "libngspice search path: %s", dllFile.GetFullPath() ); m_dll.Load( dllFile.GetFullPath(), wxDL_VERBATIM | wxDL_QUIET | wxDL_NOW ); if( m_dll.IsLoaded() ) { wxLogTrace( traceNgspice, "libngspice path found in: %s", dllFile.GetFullPath() ); break; } } } #endif if( !m_dll.IsLoaded() ) throw std::runtime_error( "Missing ngspice shared library" ); m_error = false; // Obtain function pointers m_ngSpice_Init = (ngSpice_Init) m_dll.GetSymbol( "ngSpice_Init" ); m_ngSpice_Circ = (ngSpice_Circ) m_dll.GetSymbol( "ngSpice_Circ" ); m_ngSpice_Command = (ngSpice_Command) m_dll.GetSymbol( "ngSpice_Command" ); m_ngGet_Vec_Info = (ngGet_Vec_Info) m_dll.GetSymbol( "ngGet_Vec_Info" ); m_ngSpice_AllPlots = (ngSpice_AllPlots) m_dll.GetSymbol( "ngSpice_AllPlots" ); m_ngSpice_AllVecs = (ngSpice_AllVecs) m_dll.GetSymbol( "ngSpice_AllVecs" ); m_ngSpice_Running = (ngSpice_Running) m_dll.GetSymbol( "ngSpice_running" ); // it is not a typo m_ngSpice_Init( &cbSendChar, &cbSendStat, &cbControlledExit, NULL, NULL, &cbBGThreadRunning, this ); // Load a custom spinit file, to fix the problem with loading .cm files // Switch to the executable directory, so the relative paths are correct wxString cwd( wxGetCwd() ); wxFileName exeDir( stdPaths.GetExecutablePath() ); wxSetWorkingDirectory( exeDir.GetPath() ); // Find *.cm files string cmPath = findCmPath(); // __CMPATH is used in custom spinit file to point to the codemodels directory if( !cmPath.empty() ) Command( "set __CMPATH=\"" + cmPath + "\"" ); // Possible relative locations for spinit file const vector spiceinitPaths = { ".", #ifdef __WXMAC__ stdPaths.GetPluginsDir() + "/sim/ngspice/scripts", wxFileName( stdPaths.GetExecutablePath() ).GetPath() + "/../../../../../Contents/PlugIns/sim/ngspice/scripts" #endif /* __WXMAC__ */ "../share/kicad", "../share", "../../share/kicad", "../../share" }; bool foundSpiceinit = false; for( const auto& path : spiceinitPaths ) { wxLogTrace( traceNgspice, "ngspice init script search path: %s", path ); if( loadSpinit( path + "/spiceinit" ) ) { wxLogTrace( traceNgspice, "ngspice path found in: %s", path ); foundSpiceinit = true; break; } } // Last chance to load codemodel files, we have not found // spiceinit file, but we know the path to *.cm files if( !foundSpiceinit && !cmPath.empty() ) loadCodemodels( cmPath ); // Restore the working directory wxSetWorkingDirectory( cwd ); // Workarounds to avoid hang ups on certain errors // These commands have to be called, no matter what is in the spinit file Command( "unset interactive" ); Command( "set noaskquit" ); Command( "set nomoremode" ); m_initialized = true; } bool NGSPICE::loadSpinit( const string& aFileName ) { if( !wxFileName::FileExists( aFileName ) ) return false; wxTextFile file; if( !file.Open( aFileName ) ) return false; for( auto cmd = file.GetFirstLine(); !file.Eof(); cmd = file.GetNextLine() ) Command( cmd.ToStdString() ); return true; } string NGSPICE::findCmPath() const { const vector cmPaths = { #ifdef __WXMAC__ "/Applications/ngspice/lib/ngspice", "Contents/Frameworks", wxStandardPaths::Get().GetPluginsDir() + "/sim/ngspice", wxFileName( wxStandardPaths::Get().GetExecutablePath() ).GetPath() + "/../../../../../Contents/PlugIns/sim/ngspice" #endif /* __WXMAC__ */ "../lib/ngspice", "../../lib/ngspice" "lib/ngspice", "ngspice" }; for( const auto& path : cmPaths ) { wxLogTrace( traceNgspice, "ngspice code models search path: %s", path ); if( wxFileName::FileExists( path + "/spice2poly.cm" ) ) { wxLogTrace( traceNgspice, "ngspice code models found in: %s", path ); return path; } } return string(); } bool NGSPICE::loadCodemodels( const string& aPath ) { wxArrayString cmFiles; size_t count = wxDir::GetAllFiles( aPath, &cmFiles ); for( const auto& cm : cmFiles ) Command( "codemodel " + cm.ToStdString() ); return count != 0; } int NGSPICE::cbSendChar( char* what, int id, void* user ) { NGSPICE* sim = reinterpret_cast( user ); if( sim->m_reporter ) { // strip stdout/stderr from the line if( ( strncasecmp( what, "stdout ", 7 ) == 0 ) || ( strncasecmp( what, "stderr ", 7 ) == 0 ) ) what += 7; sim->m_reporter->Report( what ); } return 0; } int NGSPICE::cbSendStat( char* what, int id, void* user ) { /* NGSPICE* sim = reinterpret_cast( user ); if( sim->m_consoleReporter ) sim->m_consoleReporter->Report( what );*/ return 0; } int NGSPICE::cbBGThreadRunning( bool is_running, int id, void* user ) { NGSPICE* sim = reinterpret_cast( user ); if( sim->m_reporter ) // I know the test below seems like an error, but well, it works somehow.. sim->m_reporter->OnSimStateChange( sim, is_running ? SIM_IDLE : SIM_RUNNING ); return 0; } int NGSPICE::cbControlledExit( int status, bool immediate, bool exit_upon_quit, int id, void* user ) { // Something went wrong, reload the dll NGSPICE* sim = reinterpret_cast( user ); sim->m_error = true; return 0; } void NGSPICE::validate() { if( m_error ) { m_initialized = false; init_dll(); } } const std::string NGSPICE::GetNetlist() const { return m_netlist; } bool NGSPICE::m_initialized = false;