/* * This program source code file is part of KiCad, a free EDA CAD application. * * Copyright (C) 2016 CERN * @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 "dialog_spice_model.h" #include #include #include #include #include // Helper function to shorten conditions static bool empty( const wxTextCtrl* aCtrl ) { return aCtrl->GetValue().IsEmpty(); } // Function to sort PWL values list static int wxCALLBACK comparePwlValues( wxIntPtr aItem1, wxIntPtr aItem2, wxIntPtr WXUNUSED( aSortData ) ) { float* t1 = reinterpret_cast( &aItem1 ); float* t2 = reinterpret_cast( &aItem2 ); if( *t1 > *t2 ) return 1; if( *t1 < *t2 ) return -1; return 0; } DIALOG_SPICE_MODEL::DIALOG_SPICE_MODEL( wxWindow* aParent, SCH_COMPONENT& aComponent, SCH_FIELDS& aFields ) : DIALOG_SPICE_MODEL_BASE( aParent ), m_component( aComponent ), m_fields( aFields ), m_spiceEmptyValidator( true ), m_notEmptyValidator( wxFILTER_EMPTY ) { m_pasValue->SetValidator( m_spiceValidator ); m_semiType->SetValidator( m_notEmptyValidator ); m_semiModel->SetValidator( m_notEmptyValidator ); m_icModel->SetValidator( m_notEmptyValidator ); m_genDc->SetValidator( m_spiceEmptyValidator ); m_genAcMag->SetValidator( m_spiceEmptyValidator ); m_genAcPhase->SetValidator( m_spiceEmptyValidator ); m_pulseInit->SetValidator( m_spiceEmptyValidator ); m_pulseNominal->SetValidator( m_spiceEmptyValidator ); m_pulseDelay->SetValidator( m_spiceEmptyValidator ); m_pulseRise->SetValidator( m_spiceEmptyValidator ); m_pulseFall->SetValidator( m_spiceEmptyValidator ); m_pulseWidth->SetValidator( m_spiceEmptyValidator ); m_pulsePeriod->SetValidator( m_spiceEmptyValidator ); m_sinOffset->SetValidator( m_spiceEmptyValidator ); m_sinAmplitude->SetValidator( m_spiceEmptyValidator ); m_sinFreq->SetValidator( m_spiceEmptyValidator ); m_sinDelay->SetValidator( m_spiceEmptyValidator ); m_sinDampFactor->SetValidator( m_spiceEmptyValidator ); m_expInit->SetValidator( m_spiceEmptyValidator ); m_expPulsed->SetValidator( m_spiceEmptyValidator ); m_expRiseDelay->SetValidator( m_spiceEmptyValidator ); m_expRiseConst->SetValidator( m_spiceEmptyValidator ); m_expFallDelay->SetValidator( m_spiceEmptyValidator ); m_expFallConst->SetValidator( m_spiceEmptyValidator ); m_pwlTimeCol = m_pwlValList->AppendColumn( "Time [s]", wxLIST_FORMAT_LEFT, 100 ); m_pwlValueCol = m_pwlValList->AppendColumn( "Value [V/A]", wxLIST_FORMAT_LEFT, 100 ); m_sdbSizerOK->SetDefault(); } bool DIALOG_SPICE_MODEL::TransferDataFromWindow() { if( !DIALOG_SPICE_MODEL_BASE::TransferDataFromWindow() ) return false; wxWindow* page = m_notebook->GetCurrentPage(); // Passive if( page == m_passive ) { if( !m_passive->Validate() ) return false; switch( m_pasType->GetSelection() ) { case 0: m_fieldsTmp[SF_PRIMITIVE] = (char) SP_RESISTOR; break; case 1: m_fieldsTmp[SF_PRIMITIVE] = (char) SP_CAPACITOR; break; case 2: m_fieldsTmp[SF_PRIMITIVE] = (char) SP_INDUCTOR; break; default: wxASSERT_MSG( false, "Unhandled passive type" ); return false; break; } m_fieldsTmp[SF_MODEL] = m_pasValue->GetValue(); } // Semiconductor else if( page == m_semiconductor ) { if( !m_semiconductor->Validate() ) return false; switch( m_semiType->GetSelection() ) { case 0: m_fieldsTmp[SF_PRIMITIVE] = (char) SP_DIODE; break; case 1: m_fieldsTmp[SF_PRIMITIVE] = (char) SP_BJT; break; case 2: m_fieldsTmp[SF_PRIMITIVE] = (char) SP_MOSFET; break; default: wxASSERT_MSG( false, "Unhandled semiconductor type" ); return false; break; } m_fieldsTmp[SF_MODEL] = m_semiModel->GetValue(); if( !empty( m_semiLib ) ) m_fieldsTmp[SF_LIB_FILE] = m_semiLib->GetValue(); } // Integrated circuit else if( page == m_ic ) { if( !m_ic->Validate() ) return false; m_fieldsTmp[SF_PRIMITIVE] = (char) SP_SUBCKT; m_fieldsTmp[SF_MODEL] = m_icModel->GetValue(); if( !empty( m_icLib ) ) m_fieldsTmp[SF_LIB_FILE] = m_icLib->GetValue(); } // Power source else if( page == m_power ) { wxString model; if( !generatePowerSource( model ) ) return false; m_fieldsTmp[SF_PRIMITIVE] = (char)( m_pwrType->GetSelection() ? SP_ISOURCE : SP_VSOURCE ); m_fieldsTmp[SF_MODEL] = model; } else { wxASSERT_MSG( false, "Unhandled model type" ); return false; } m_fieldsTmp[SF_ENABLED] = !m_disabled->GetValue() ? "Y" : "N"; // note bool inversion m_fieldsTmp[SF_NODE_SEQUENCE] = m_nodeSeqCheck->IsChecked() ? m_nodeSeqVal->GetValue() : ""; // Apply the settings for( int i = 0; i < SF_END; ++i ) { if( m_fieldsTmp.count( (SPICE_FIELD) i ) > 0 && !m_fieldsTmp.at( i ).IsEmpty() ) { getField( i ).SetText( m_fieldsTmp[i] ); } else { // Erase empty fields (having empty fields causes a warning in the properties dialog) const wxString& spiceField = NETLIST_EXPORTER_PSPICE::GetSpiceFieldName( (SPICE_FIELD) i ); auto fieldIt = std::find_if( m_fields.begin(), m_fields.end(), [&]( const SCH_FIELD& f ) { return f.GetName() == spiceField; } ); if( fieldIt != m_fields.end() ) m_fields.erase( fieldIt ); } } return true; } bool DIALOG_SPICE_MODEL::TransferDataToWindow() { const auto& spiceFields = NETLIST_EXPORTER_PSPICE::GetSpiceFields(); // Fill out the working buffer for( unsigned int idx = 0; idx < spiceFields.size(); ++idx ) { const wxString& spiceField = spiceFields[idx]; auto fieldIt = std::find_if( m_fields.begin(), m_fields.end(), [&]( const SCH_FIELD& f ) { return f.GetName() == spiceField; } ); // Do not modify the existing value, just add missing fields with default values if( fieldIt != m_fields.end() && !fieldIt->GetText().IsEmpty() ) m_fieldsTmp[idx] = fieldIt->GetText(); else m_fieldsTmp[idx] = NETLIST_EXPORTER_PSPICE::GetSpiceFieldDefVal( (SPICE_FIELD) idx, &m_component, NET_ADJUST_INCLUDE_PATHS | NET_ADJUST_PASSIVE_VALS ); } // Analyze the component fields to fill out the dialog char primitive = toupper( m_fieldsTmp[SF_PRIMITIVE][0] ); switch( primitive ) { case SP_RESISTOR: case SP_CAPACITOR: case SP_INDUCTOR: m_notebook->SetSelection( m_notebook->FindPage( m_passive ) ); m_pasType->SetSelection( primitive == SP_RESISTOR ? 0 : primitive == SP_CAPACITOR ? 1 : primitive == SP_INDUCTOR ? 2 : -1 ); m_pasValue->SetValue( m_fieldsTmp[SF_MODEL] ); break; case SP_DIODE: case SP_BJT: case SP_MOSFET: m_notebook->SetSelection( m_notebook->FindPage( m_semiconductor ) ); m_semiType->SetSelection( primitive == SP_DIODE ? 0 : primitive == SP_BJT ? 1 : primitive == SP_MOSFET ? 2 : -1 ); m_semiModel->SetValue( m_fieldsTmp[SF_MODEL] ); m_semiLib->SetValue( m_fieldsTmp[SF_LIB_FILE] ); if( !empty( m_semiLib ) ) { const wxString& libFile = m_semiLib->GetValue(); m_fieldsTmp[SF_LIB_FILE] = libFile; updateFromFile( m_semiModel, libFile, ".model" ); } break; case SP_SUBCKT: m_notebook->SetSelection( m_notebook->FindPage( m_ic ) ); m_icModel->SetValue( m_fieldsTmp[SF_MODEL] ); m_icLib->SetValue( m_fieldsTmp[SF_LIB_FILE] ); if( !empty( m_icLib ) ) { const wxString& libFile = m_icLib->GetValue(); m_fieldsTmp[SF_LIB_FILE] = libFile; updateFromFile( m_icModel, libFile, ".subckt" ); } break; case SP_VSOURCE: case SP_ISOURCE: if( !parsePowerSource( m_fieldsTmp[SF_MODEL] ) ) return false; m_notebook->SetSelection( m_notebook->FindPage( m_power ) ); m_pwrType->SetSelection( primitive == SP_ISOURCE ? 1 : 0 ); break; default: //wxASSERT_MSG( false, "Unhandled Spice primitive type" ); break; } m_disabled->SetValue( !NETLIST_EXPORTER_PSPICE::StringToBool( m_fieldsTmp[SF_ENABLED] ) ); // Check if node sequence is different than the default one if( m_fieldsTmp[SF_NODE_SEQUENCE] != NETLIST_EXPORTER_PSPICE::GetSpiceFieldDefVal( SF_NODE_SEQUENCE, &m_component, 0 ) ) { m_nodeSeqCheck->SetValue( true ); m_nodeSeqVal->SetValue( m_fieldsTmp[SF_NODE_SEQUENCE] ); } return DIALOG_SPICE_MODEL_BASE::TransferDataToWindow(); } bool DIALOG_SPICE_MODEL::parsePowerSource( const wxString& aModel ) { if( aModel.IsEmpty() ) return false; wxStringTokenizer tokenizer( aModel, " ()" ); wxString tkn = tokenizer.GetNextToken().Lower(); while( tokenizer.HasMoreTokens() ) { // Variables used for generic values processing (filling out wxTextCtrls in sequence) bool genericProcessing = false; unsigned int genericReqParamsCount = 0; std::vector genericControls; if( tkn == "dc" ) { // There might be an optional "dc" or "trans" directive, skip it if( tkn == "dc" || tkn == "trans" ) tkn = tokenizer.GetNextToken().Lower(); // DC value try { m_genDc->SetValue( SPICE_VALUE( tkn ).ToSpiceString() ); } catch( ... ) { return false; } } else if( tkn == "ac" ) { // AC magnitude try { tkn = tokenizer.GetNextToken().Lower(); m_genAcMag->SetValue( SPICE_VALUE( tkn ).ToSpiceString() ); } catch( ... ) { return false; } // AC phase (optional) try { tkn = tokenizer.GetNextToken().Lower(); m_genAcPhase->SetValue( SPICE_VALUE( tkn ).ToSpiceString() ); } catch( ... ) { continue; // perhaps another directive } } else if( tkn == "pulse" ) { m_powerNotebook->SetSelection( m_powerNotebook->FindPage( m_pwrPulse ) ); genericProcessing = true; genericReqParamsCount = 2; genericControls = { m_pulseInit, m_pulseNominal, m_pulseDelay, m_pulseRise, m_pulseFall, m_pulseWidth, m_pulsePeriod }; } else if( tkn == "sin" ) { m_powerNotebook->SetSelection( m_powerNotebook->FindPage( m_pwrSin ) ); genericProcessing = true; genericReqParamsCount = 2; genericControls = { m_sinOffset, m_sinAmplitude, m_sinFreq, m_sinDelay, m_sinDampFactor }; } else if( tkn == "exp" ) { m_powerNotebook->SetSelection( m_powerNotebook->FindPage( m_pwrExp ) ); genericProcessing = true; genericReqParamsCount = 2; genericControls = { m_expInit, m_expPulsed, m_expRiseDelay, m_expRiseConst, m_expFallDelay, m_expFallConst }; } else if( tkn == "pwl" ) { m_powerNotebook->SetSelection( m_powerNotebook->FindPage( m_pwrPwl ) ); try { while( tokenizer.HasMoreTokens() ) { tkn = tokenizer.GetNextToken(); SPICE_VALUE time( tkn ); tkn = tokenizer.GetNextToken(); SPICE_VALUE value( tkn ); addPwlValue( time.ToSpiceString(), value.ToSpiceString() ); } } catch( ... ) { return false; } } else { // Unhandled power source type wxASSERT_MSG( false, "Unhandled power source type" ); return false; } if( genericProcessing ) { try { for( unsigned int i = 0; i < genericControls.size(); ++i ) { // If there are no more tokens, let's check if we got at least required fields if( !tokenizer.HasMoreTokens() ) return ( i >= genericReqParamsCount ); tkn = tokenizer.GetNextToken().Lower(); genericControls[i]->SetValue( SPICE_VALUE( tkn ).ToSpiceString() ); } } catch( ... ) { return false; } } // Get the next token now, so if any of the branches catches an expection, try to // process it in another branch tkn = tokenizer.GetNextToken().Lower(); } return true; } bool DIALOG_SPICE_MODEL::generatePowerSource( wxString& aTarget ) const { wxString acdc, trans; wxWindow* page = m_powerNotebook->GetCurrentPage(); bool useTrans = true; // shall we use the transient command part? // Variables for generic processing bool genericProcessing = false; unsigned int genericReqParamsCount = 0; std::vector genericControls; /// DC / AC section // If SPICE_VALUE can be properly constructed, then it is a valid value try { if( !empty( m_genDc ) ) acdc += wxString::Format( "dc %s ", SPICE_VALUE( m_genDc->GetValue() ).ToSpiceString() ); } catch( ... ) { DisplayError( NULL, wxT( "Invalid DC value" ) ); return false; } try { if( !empty( m_genAcMag ) ) { acdc += wxString::Format( "ac %s ", SPICE_VALUE( m_genAcMag->GetValue() ).ToSpiceString() ); if( !empty( m_genAcPhase ) ) acdc += wxString::Format( "%s ", SPICE_VALUE( m_genAcPhase->GetValue() ).ToSpiceString() ); } } catch( ... ) { DisplayError( NULL, wxT( "Invalid AC magnitude or phase" ) ); return false; } /// Transient section if( page == m_pwrPulse ) { if( !m_pwrPulse->Validate() ) return false; genericProcessing = true; trans += "pulse"; genericReqParamsCount = 2; genericControls = { m_pulseInit, m_pulseNominal, m_pulseDelay, m_pulseRise, m_pulseFall, m_pulseWidth, m_pulsePeriod }; } else if( page == m_pwrSin ) { if( !m_pwrSin->Validate() ) return false; genericProcessing = true; trans += "sin"; genericReqParamsCount = 2; genericControls = { m_sinOffset, m_sinAmplitude, m_sinFreq, m_sinDelay, m_sinDampFactor }; } else if( page == m_pwrExp ) { if( !m_pwrExp->Validate() ) return false; genericProcessing = true; trans += "exp"; genericReqParamsCount = 2; genericControls = { m_expInit, m_expPulsed, m_expRiseDelay, m_expRiseConst, m_expFallDelay, m_expFallConst }; } else if( page == m_pwrPwl ) { if( m_pwlValList->GetItemCount() > 0 ) { trans += "pwl("; for( int i = 0; i < m_pwlValList->GetItemCount(); ++i ) { trans += wxString::Format( "%s %s ", m_pwlValList->GetItemText( i, m_pwlTimeCol ), m_pwlValList->GetItemText( i, m_pwlValueCol ) ); } trans.Trim(); trans += ")"; } } if( genericProcessing ) { bool finished = false; unsigned int paramCounter = 0; trans += "("; for( auto textCtrl : genericControls ) { if( empty( textCtrl ) ) { finished = true; if( paramCounter < genericReqParamsCount ) { if( paramCounter == 0 ) { // It is fine, no parameters were entered useTrans = false; break; } DisplayError( NULL, wxString::Format( wxT( "You need to specify at least the " "first %d parameters for the transient source" ), genericReqParamsCount ) ); return false; } } else if( finished ) { DisplayError( NULL, wxT( "You cannot leave interleaved blank " "spaces for the transient source" ) ); return false; } trans += wxString::Format( "%s ", textCtrl->GetValue() ); ++paramCounter; } trans.Trim(); trans += ")"; } aTarget = acdc; if( useTrans ) aTarget += trans; aTarget.Trim( false ); aTarget.Trim( true ); return true; } void DIALOG_SPICE_MODEL::updateFromFile( wxComboBox* aComboBox, const wxString& aFilePath, const wxString& aKeyword ) { wxString curValue = aComboBox->GetValue(); const wxString keyword( aKeyword.Lower() ); wxFileName filePath( aFilePath ); if( !filePath.Exists() ) { // Look for the file in the project path filePath.SetPath( Prj().GetProjectPath() + filePath.GetPath() ); } wxTextFile file; if( !file.Open( filePath.GetFullPath() ) ) return; aComboBox->Clear(); // Process the file, looking for components for( wxString line = file.GetFirstLine().Lower(); !file.Eof(); line = file.GetNextLine() ) { wxStringTokenizer tokenizer( line, " " ); while( tokenizer.HasMoreTokens() ) { wxString token = tokenizer.GetNextToken().Lower(); if( token == keyword ) { token = tokenizer.GetNextToken(); if( !token.IsEmpty() ) aComboBox->Append( token ); } } } // Restore the previous value or if there is none - pick the first one from the loaded library if( !curValue.IsEmpty() ) aComboBox->SetValue( curValue ); else if( aComboBox->GetCount() > 0 ) aComboBox->SetSelection( 0 ); } SCH_FIELD& DIALOG_SPICE_MODEL::getField( int aFieldType ) { const wxString& spiceField = NETLIST_EXPORTER_PSPICE::GetSpiceFieldName( (SPICE_FIELD) aFieldType ); auto fieldIt = std::find_if( m_fields.begin(), m_fields.end(), [&]( const SCH_FIELD& f ) { return f.GetName() == spiceField; } ); // Found one, so return it if( fieldIt != m_fields.end() ) return *fieldIt; // Create a new field with requested name m_fields.emplace_back( wxPoint(), m_fields.size(), &m_component, spiceField ); return m_fields.back(); } bool DIALOG_SPICE_MODEL::addPwlValue( const wxString& aTime, const wxString& aValue ) { // TODO execute validators if( aTime.IsEmpty() || aValue.IsEmpty() ) return false; long idx = m_pwlValList->InsertItem( m_pwlTimeCol, aTime ); m_pwlValList->SetItem( idx, m_pwlValueCol, aValue ); // There is no wxString::ToFloat, but we need to guarantee it fits in 4 bytes double timeD; float timeF; m_pwlTime->GetValue().ToDouble( &timeD ); timeF = timeD; // Store the time value, so the entries can be sorted m_pwlValList->SetItemData( idx, *reinterpret_cast( &timeF ) ); // Sort items by timestamp m_pwlValList->SortItems( comparePwlValues, -1 ); return true; } void DIALOG_SPICE_MODEL::onSemiSelectLib( wxCommandEvent& event ) { wxString searchPath = wxFileName( m_semiLib->GetValue() ).GetPath(); if( searchPath.IsEmpty() ) searchPath = Prj().GetProjectPath(); wxFileDialog openDlg( this, wxT( "Select library" ), searchPath, "", "Spice library file (*.lib)|*.lib;*.LIB|Any file|*", wxFD_OPEN | wxFD_FILE_MUST_EXIST ); if( openDlg.ShowModal() == wxID_CANCEL ) return; wxFileName libPath( openDlg.GetPath() ); // Try to convert the path to relative to project if( libPath.MakeRelativeTo( Prj().GetProjectPath() ) && !libPath.GetFullPath().StartsWith( ".." ) ) m_semiLib->SetValue( libPath.GetFullPath() ); else m_semiLib->SetValue( openDlg.GetPath() ); updateFromFile( m_semiModel, openDlg.GetPath(), ".model" ); m_semiModel->Popup(); } void DIALOG_SPICE_MODEL::onSelectIcLib( wxCommandEvent& event ) { wxString searchPath = wxFileName( m_icLib->GetValue() ).GetPath(); if( searchPath.IsEmpty() ) searchPath = Prj().GetProjectPath(); wxFileDialog openDlg( this, wxT( "Select library" ), searchPath, "", "Spice library file (*.lib)|*.lib;*.LIB|Any file|*", wxFD_OPEN | wxFD_FILE_MUST_EXIST ); if( openDlg.ShowModal() == wxID_CANCEL ) return; wxFileName libPath( openDlg.GetPath() ); // Try to convert the path to relative to project if( libPath.MakeRelativeTo( Prj().GetProjectPath() ) && !libPath.GetFullPath().StartsWith( ".." ) ) m_icLib->SetValue( libPath.GetFullPath() ); else m_icLib->SetValue( openDlg.GetPath() ); updateFromFile( m_icModel, openDlg.GetPath(), ".subckt" ); m_icModel->Popup(); } void DIALOG_SPICE_MODEL::onPwlAdd( wxCommandEvent& event ) { addPwlValue( m_pwlTime->GetValue(), m_pwlValue->GetValue() ); } void DIALOG_SPICE_MODEL::onPwlRemove( wxCommandEvent& event ) { long idx = m_pwlValList->GetNextItem( -1, wxLIST_NEXT_ALL, wxLIST_STATE_SELECTED ); m_pwlValList->DeleteItem( idx ); }