kicad/pcbnew/router/pns_walkaround.cpp

381 lines
11 KiB
C++

/*
* KiRouter - a push-and-(sometimes-)shove PCB router
*
* Copyright (C) 2013-2014 CERN
* Copyright (C) 2016-2023 KiCad Developers, see AUTHORS.txt for contributors.
* Author: Tomasz Wlostowski <tomasz.wlostowski@cern.ch>
*
* 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, see <http://www.gnu.org/licenses/>.
*/
#include <optional>
#include <geometry/shape_line_chain.h>
#include "pns_walkaround.h"
#include "pns_optimizer.h"
#include "pns_router.h"
#include "pns_debug_decorator.h"
#include "pns_solid.h"
namespace PNS {
void WALKAROUND::start( const LINE& aInitialPath )
{
m_iteration = 0;
for( int pol = 0 ; pol < MaxWalkPolicies; pol++)
{
m_currentResult.status[ pol ] = ST_IN_PROGRESS;
m_currentResult.lines[ pol ] = aInitialPath;
m_currentResult.lines[ pol ].ClearLinks();
}
}
NODE::OPT_OBSTACLE WALKAROUND::nearestObstacle( const LINE& aPath )
{
COLLISION_SEARCH_OPTIONS opts;
opts.m_kindMask = m_itemMask;
if( ! m_restrictedSet.empty() )
{
opts.m_filter = [ this ] ( const ITEM* item ) -> bool
{
if( m_restrictedSet.find( item ) != m_restrictedSet.end() )
return true;
return false;
};
}
opts.m_useClearanceEpsilon = false;
return m_world->NearestObstacle( &aPath, opts );
}
void WALKAROUND::RestrictToCluster( bool aEnabled, const TOPOLOGY::CLUSTER& aCluster )
{
m_restrictedVertices.clear();
m_restrictedSet.clear();
if( aEnabled )
{
for( ITEM* item : aCluster.m_items )
{
m_restrictedSet.insert( item );
if ( item->HasHole() )
m_restrictedSet.insert( item->Hole() );
}
}
for( ITEM* item : aCluster.m_items )
{
if( SOLID* solid = dyn_cast<SOLID*>( item ) )
m_restrictedVertices.push_back( solid->Anchor( 0 ) );
}
}
static wxString policy2string ( WALKAROUND::WALK_POLICY policy )
{
switch(policy)
{
case WALKAROUND::WP_CCW: return wxT("ccw");
case WALKAROUND::WP_CW: return wxT("cw");
case WALKAROUND::WP_SHORTEST: return wxT("shortest");
}
return wxT("?");
}
bool WALKAROUND::singleStep()
{
TOPOLOGY topo( m_world );
TOPOLOGY::CLUSTER pendingClusters[MaxWalkPolicies];
for( int i = 0; i < MaxWalkPolicies; i++ )
{
if( !m_enabledPolicies[i] )
continue;
auto& line = m_currentResult.lines[ i ];
auto& status = m_currentResult.status[ i ];
PNS_DBG( Dbg(), AddItem, &line, WHITE, 10000, wxString::Format( "current (policy %d, stat %d)", i, status ) );
if( status != ST_IN_PROGRESS )
continue;
auto obstacle = nearestObstacle( line );
if( !obstacle )
{
m_currentResult.status[ i ] = ST_DONE;
PNS_DBG( Dbg(), Message, wxString::Format( "no-more-colls pol %d st %d", i, status ) );
continue;
}
pendingClusters[ i ] = topo.AssembleCluster( obstacle->m_item, line.Layer(), 0.0, line.Net() );
PNS_DBG( Dbg(), AddItem, obstacle->m_item, BLUE, 10000, wxString::Format( "col-item owner-depth %d cl-items=%d", static_cast<const NODE*>( obstacle->m_item->Owner() )->Depth(), (int) pendingClusters[i].m_items.size() ) );
}
DIRECTION_45::CORNER_MODE cornerMode = Settings().GetCornerMode();
auto processCluster = [ & ] ( TOPOLOGY::CLUSTER& aCluster, LINE& aLine, bool aCw ) -> bool
{
PNS_DBG( Dbg(), BeginGroup, wxString::Format( "cluster-details [cw %d]", aCw?1:0 ), 1 );
for( auto& clItem : aCluster.m_items )
{
int clearance = m_world->GetClearance( clItem, &aLine, false );
SHAPE_LINE_CHAIN hull = clItem->Hull( clearance + 1000, aLine.Width(), aLine.Layer() );
if( cornerMode == DIRECTION_45::MITERED_90 || cornerMode == DIRECTION_45::ROUNDED_90 )
{
BOX2I bbox = hull.BBox();
hull.Clear();
hull.Append( bbox.GetLeft(), bbox.GetTop() );
hull.Append( bbox.GetRight(), bbox.GetTop() );
hull.Append( bbox.GetRight(), bbox.GetBottom() );
hull.Append( bbox.GetLeft(), bbox.GetBottom() );
}
LINE tmp( aLine );
bool stat = aLine.Walkaround( hull, tmp.Line(), aCw );
PNS_DBG( Dbg(), AddShape, &hull, YELLOW, 10000, wxString::Format( "hull stat %d", stat?1:0 ) );
PNS_DBG( Dbg(), AddItem, &tmp, RED, 10000, wxString::Format( "walk stat %d", stat?1:0 ) );
PNS_DBG( Dbg(), AddItem, clItem, WHITE, 10000, wxString::Format( "item stat %d", stat?1:0 ) );
if( !stat )
{
PNS_DBGN( Dbg(), EndGroup );
return false;
}
aLine.SetShape( tmp.CLine() );
}
PNS_DBGN( Dbg(), EndGroup );
return true;
};
if ( m_enabledPolicies[WP_CW] )
{
bool stat = processCluster( pendingClusters[ WP_CW ], m_currentResult.lines[ WP_CW ], true );
if( !stat )
m_currentResult.status[ WP_CW ] = ST_STUCK;
}
if ( m_enabledPolicies[WP_CCW] )
{
bool stat = processCluster( pendingClusters[ WP_CCW ], m_currentResult.lines[ WP_CCW ], false );
if( !stat )
m_currentResult.status[ WP_CCW ] = ST_STUCK;
}
if( m_enabledPolicies[WP_SHORTEST] )
{
LINE& line = m_currentResult.lines[WP_SHORTEST];
LINE path_cw( line ), path_ccw( line );
auto st_cw = processCluster( pendingClusters[WP_SHORTEST], path_cw, true );
auto st_ccw = processCluster( pendingClusters[WP_SHORTEST], path_ccw, false );
bool cw_coll = st_cw ? m_world->CheckColliding( &path_cw ).has_value() : false;
bool ccw_coll = st_ccw ? m_world->CheckColliding( &path_ccw ).has_value() : false;
double lengthFactorCw = (double) path_cw.CLine().Length() / (double) m_initialLength;
double lengthFactorCcw = (double) path_ccw.CLine().Length() / (double) m_initialLength;
PNS_DBG( Dbg(), AddItem, &path_cw, RED, 10000, wxString::Format( "shortest-cw stat %d lf %.1f", st_cw?1:0, lengthFactorCw ) );
PNS_DBG( Dbg(), AddItem, &path_ccw, BLUE, 10000, wxString::Format( "shortest-ccw stat %d lf %.1f", st_ccw?1:0, lengthFactorCcw ) );
std::optional<LINE> shortest;
std::optional<LINE> shortest_alt;
if( st_cw && st_ccw )
{
if( !cw_coll && !ccw_coll || ( cw_coll && ccw_coll) )
{
if( path_cw.CLine().Length() > path_ccw.CLine().Length() )
{
shortest = path_ccw;
shortest_alt = path_cw;
}
else
{
shortest = path_cw;
shortest_alt = path_ccw;
}
}
else if( !cw_coll )
shortest = path_cw;
else if( !ccw_coll )
shortest = path_ccw;
}
else if( st_ccw )
shortest = path_ccw;
else if( st_cw )
shortest = path_cw;
bool anyColliding = false;
if( m_lastShortestCluster && shortest.has_value() )
{
for( auto& item : m_lastShortestCluster->m_items )
{
if( shortest->Collide( item, m_world, shortest->Layer() ) )
{
anyColliding = true;
break;
}
}
PNS_DBG( Dbg(), Message, wxString::Format("check-back cc %d items %d coll %d", (int) pendingClusters[ WP_SHORTEST ].m_items.size(), (int) m_lastShortestCluster->m_items.size(), anyColliding ? 1: 0 ) );
}
if ( anyColliding )
{
shortest = shortest_alt;
}
if( !shortest )
{
m_currentResult.status[WP_SHORTEST] = ST_STUCK;
}
else
{
m_currentResult.lines[WP_SHORTEST] = *shortest;
}
m_lastShortestCluster = pendingClusters[ WP_SHORTEST ];
}
return ST_IN_PROGRESS;
}
const WALKAROUND::RESULT WALKAROUND::Route( const LINE& aInitialPath )
{
RESULT result;
m_initialLength = aInitialPath.CLine().Length();
// special case for via-in-the-middle-of-track placement
#if 0
if( aInitialPath.PointCount() <= 1 )
{
if( aInitialPath.EndsWithVia() && m_world->CheckColliding( &aInitialPath.Via(),
m_itemMask ) )
{
// fixme restult
}
//return RESULT( STUCK, STUCK );
return RESULT(); //( DONE, DONE, aInitialPath, aInitialPath );
}
#endif
start( aInitialPath );
PNS_DBG( Dbg(), AddItem, &aInitialPath, WHITE, 10000, wxT( "initial-path" ) );
while( m_iteration < m_iterationLimit )
{
singleStep();
bool stillInProgress = false;
for( int pol = 0; pol < MaxWalkPolicies; pol++ )
{
if (!m_enabledPolicies[pol])
continue;
auto& st = m_currentResult.status[pol];
auto& ln = m_currentResult.lines[pol];
double lengthFactor = (double) ln.CLine().Length() / (double) aInitialPath.CLine().Length();
// In some situations, there isn't a trivial path (or even a path at all). Hitting the
// iteration limit causes lag, so we can exit out early if the walkaround path gets very long
// compared with the initial path. If the length exceeds the initial length times this factor,
// fail out.
if( m_lengthLimitOn )
{
if( st != ST_DONE && lengthFactor > m_lengthExpansionFactor )
st = ST_ALMOST_DONE;
}
PNS_DBG( Dbg(), Message, wxString::Format( "check-wp iter %d st %d i %d lf %.1f", m_iteration, st, pol, lengthFactor ) );
if ( st == ST_IN_PROGRESS )
stillInProgress = true;
}
if( !stillInProgress )
break;
m_iteration++;
}
for( int pol = 0; pol < MaxWalkPolicies; pol++ )
{
auto& st = m_currentResult.status[pol];
const auto& ln = m_currentResult.lines[pol].CLine();
m_currentResult.lines[pol].ClearLinks();
if( st == ST_IN_PROGRESS )
st = ST_ALMOST_DONE;
if( ln.SegmentCount() < 1 || ln.CPoint( 0 ) != aInitialPath.CPoint( 0 ) )
{
st = ST_STUCK;
}
if( ln.PointCount() > 0 && ln.CPoint( -1 ) != aInitialPath.CPoint( -1 ) )
{
st = ST_ALMOST_DONE;
}
PNS_DBG( Dbg(), Message, wxString::Format( "stat=%d", st ) );
}
return m_currentResult;
}
void WALKAROUND::SetAllowedPolicies( std::vector<WALK_POLICY> aPolicies)
{
for( int i = 0; i < MaxWalkPolicies; i++ )
m_enabledPolicies[i] = false;
for ( auto p : aPolicies )
m_enabledPolicies[p] = true;
}
}