kicad/pcbnew/ratsnest_data.cpp

1316 lines
36 KiB
C++
Raw Normal View History

2013-11-25 15:50:03 +00:00
/*
* This program source code file is part of KICAD, a free EDA CAD application.
*
* Copyright (C) 2013-2015 CERN
* Copyright (C) 2016 KiCad Developers, see AUTHORS.txt for contributors.
2013-11-25 15:50:03 +00:00
* @author Maciej Suminski <maciej.suminski@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 2
* 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:
* http://www.gnu.org/licenses/old-licenses/gpl-2.0.html
* or you may search the http://www.gnu.org website for the version 2 license,
* or you may write to the Free Software Foundation, Inc.,
* 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA
*/
/**
* @file ratsnest_data.cpp
* @brief Class that computes missing connections on a PCB.
*/
#ifdef USE_OPENMP
#include <omp.h>
#endif /* USE_OPENMP */
2013-11-25 15:50:03 +00:00
#include <ratsnest_data.h>
#include <class_board.h>
#include <class_module.h>
#include <class_pad.h>
#include <class_track.h>
#include <class_zone.h>
#include <functional>
using namespace std::placeholders;
2013-11-25 15:50:03 +00:00
#include <geometry/shape_poly_set.h>
2013-11-25 15:50:03 +00:00
#include <cassert>
#include <algorithm>
#include <limits>
2015-06-04 12:54:06 +00:00
#ifdef PROFILE
#include <profile.h>
#endif
static uint64_t getDistance( const RN_NODE_PTR& aNode1, const RN_NODE_PTR& aNode2 )
2013-11-25 15:50:03 +00:00
{
// Drop the least significant bits to avoid overflow
int64_t x = ( aNode1->GetX() - aNode2->GetX() ) >> 16;
int64_t y = ( aNode1->GetY() - aNode2->GetY() ) >> 16;
// We do not need sqrt() here, as the distance is computed only for comparison
return ( x * x + y * y );
}
static bool sortDistance( const RN_NODE_PTR& aOrigin, const RN_NODE_PTR& aNode1,
2013-11-25 15:50:03 +00:00
const RN_NODE_PTR& aNode2 )
{
return getDistance( aOrigin, aNode1 ) < getDistance( aOrigin, aNode2 );
}
static bool sortWeight( const RN_EDGE_PTR& aEdge1, const RN_EDGE_PTR& aEdge2 )
2013-11-25 15:50:03 +00:00
{
return aEdge1->GetWeight() < aEdge2->GetWeight();
2013-11-25 15:50:03 +00:00
}
bool sortArea( const RN_POLY& aP1, const RN_POLY& aP2 )
{
return aP1.m_bbox.GetArea() < aP2.m_bbox.GetArea();
}
bool operator==( const RN_NODE_PTR& aFirst, const RN_NODE_PTR& aSecond )
{
return aFirst->GetX() == aSecond->GetX() && aFirst->GetY() == aSecond->GetY();
}
bool operator!=( const RN_NODE_PTR& aFirst, const RN_NODE_PTR& aSecond )
{
return aFirst->GetX() != aSecond->GetX() || aFirst->GetY() != aSecond->GetY();
}
RN_NODE_AND_FILTER operator&&( const RN_NODE_FILTER& aFilter1, const RN_NODE_FILTER& aFilter2 )
{
return RN_NODE_AND_FILTER( aFilter1, aFilter2 );
}
RN_NODE_OR_FILTER operator||( const RN_NODE_FILTER& aFilter1, const RN_NODE_FILTER& aFilter2 )
{
return RN_NODE_OR_FILTER( aFilter1, aFilter2 );
}
static bool isEdgeConnectingNode( const RN_EDGE_PTR& aEdge, const RN_NODE_PTR& aNode )
2013-11-25 15:50:03 +00:00
{
return aEdge->GetSourceNode() == aNode || aEdge->GetTargetNode() == aNode;
2013-11-25 15:50:03 +00:00
}
static std::vector<RN_EDGE_MST_PTR>* kruskalMST( RN_LINKS::RN_EDGE_LIST& aEdges,
std::vector<RN_NODE_PTR>& aNodes )
2013-11-25 15:50:03 +00:00
{
unsigned int nodeNumber = aNodes.size();
unsigned int mstExpectedSize = nodeNumber - 1;
unsigned int mstSize = 0;
bool ratsnestLines = false;
2013-11-25 15:50:03 +00:00
// The output
std::vector<RN_EDGE_MST_PTR>* mst = new std::vector<RN_EDGE_MST_PTR>;
2013-11-25 15:50:03 +00:00
mst->reserve( mstExpectedSize );
// Set tags for marking cycles
2016-09-19 15:25:31 +00:00
std::unordered_map<RN_NODE_PTR, int> tags;
2013-11-25 15:50:03 +00:00
unsigned int tag = 0;
for( RN_NODE_PTR& node : aNodes )
{
node->SetTag( tag );
2013-11-25 15:50:03 +00:00
tags[node] = tag++;
}
2013-11-25 15:50:03 +00:00
// Lists of nodes connected together (subtrees) to detect cycles in the graph
std::vector<std::list<int> > cycles( nodeNumber );
for( unsigned int i = 0; i < nodeNumber; ++i )
cycles[i].push_back( i );
// Kruskal algorithm requires edges to be sorted by their weight
aEdges.sort( sortWeight );
while( mstSize < mstExpectedSize && !aEdges.empty() )
{
RN_EDGE_PTR& dt = aEdges.front();
2013-11-25 15:50:03 +00:00
int srcTag = tags[dt->GetSourceNode()];
int trgTag = tags[dt->GetTargetNode()];
2013-11-25 15:50:03 +00:00
// Check if by adding this edge we are going to join two different forests
if( srcTag != trgTag )
{
2015-06-04 12:54:07 +00:00
// Because edges are sorted by their weight, first we always process connected
// items (weight == 0). Once we stumble upon an edge with non-zero weight,
// it means that the rest of the lines are ratsnest.
if( !ratsnestLines && dt->GetWeight() != 0 )
ratsnestLines = true;
2013-11-25 15:50:03 +00:00
// Update tags
if( ratsnestLines )
2014-03-05 13:57:14 +00:00
{
for( auto it = cycles[trgTag].begin(); it != cycles[trgTag].end(); ++it )
{
tags[aNodes[*it]] = srcTag;
}
2014-03-05 13:57:14 +00:00
// Do a copy of edge, but make it RN_EDGE_MST. In contrary to RN_EDGE,
// RN_EDGE_MST saves both source and target node and does not require any other
// edges to exist for getting source/target nodes
RN_EDGE_MST_PTR newEdge = std::make_shared<RN_EDGE_MST>( dt->GetSourceNode(),
dt->GetTargetNode(),
dt->GetWeight() );
assert( newEdge->GetSourceNode()->GetTag() != newEdge->GetTargetNode()->GetTag() );
assert( newEdge->GetWeight() > 0 );
2014-03-05 13:57:14 +00:00
mst->push_back( newEdge );
2013-11-25 15:50:03 +00:00
++mstSize;
}
else
{
//for( it = cycles[trgTag].begin(), itEnd = cycles[trgTag].end(); it != itEnd; ++it )
//for( auto it : cycles[trgTag] )
for( auto it = cycles[trgTag].begin(); it != cycles[trgTag].end(); ++it )
{
tags[aNodes[*it]] = srcTag;
aNodes[*it]->SetTag( srcTag );
}
// Processing a connection, decrease the expected size of the ratsnest MST
--mstExpectedSize;
}
// Move nodes that were marked with old tag to the list marked with the new tag
cycles[srcTag].splice( cycles[srcTag].end(), cycles[trgTag] );
2013-11-25 15:50:03 +00:00
}
// Remove the edge that was just processed
aEdges.erase( aEdges.begin() );
}
// Probably we have discarded some of edges, so reduce the size
mst->resize( mstSize );
return mst;
}
void RN_NET::validateEdge( RN_EDGE_MST_PTR& aEdge )
2013-11-25 15:50:03 +00:00
{
RN_NODE_PTR source = aEdge->GetSourceNode();
RN_NODE_PTR target = aEdge->GetTargetNode();
bool update = false, changed = false;
2013-11-25 15:50:03 +00:00
// If any of nodes belonging to the edge has the flag set,
// change it to the closest node that has flag cleared
// note: finding the right nodes can be done iteratively to get the best results,
// but it is not likely to be worth the time cost
do
2013-11-25 15:50:03 +00:00
{
if( changed || source->GetNoLine() )
2013-11-25 15:50:03 +00:00
{
changed = false;
std::list<RN_NODE_PTR> closest = GetClosestNodes( target,
LINE_TARGET_SAME_TAG( source->GetTag() ) );
if( !closest.empty() )
2013-11-25 15:50:03 +00:00
{
RN_NODE_PTR& node = closest.front();
if( node != source )
{
changed = true;
update = true;
source = node;
}
2013-11-25 15:50:03 +00:00
}
}
if( changed || target->GetNoLine() )
2013-11-25 15:50:03 +00:00
{
changed = false;
std::list<RN_NODE_PTR> closest = GetClosestNodes( source,
LINE_TARGET_SAME_TAG( target->GetTag() ) );
if( !closest.empty() )
2013-11-25 15:50:03 +00:00
{
RN_NODE_PTR& node = closest.front();
if( node != target )
{
changed = true;
update = true;
target = node;
}
2013-11-25 15:50:03 +00:00
}
}
}
while( changed );
assert( source->GetTag() >= 0 && target->GetTag() >= 0 );
assert( source->GetTag() != target->GetTag() );
assert( source != target );
2013-11-25 15:50:03 +00:00
// Replace an invalid edge with new, valid one
if( update )
2013-11-25 15:50:03 +00:00
aEdge.reset( new RN_EDGE_MST( source, target ) );
}
void RN_NET::removeNode( RN_NODE_PTR& aNode, const BOARD_CONNECTED_ITEM* aParent )
{
aNode->RemoveParent( aParent );
if( m_links.RemoveNode( aNode ) )
{
clearNode( aNode );
m_dirty = true;
}
}
void RN_NET::removeEdge( RN_EDGE_MST_PTR& aEdge, const BOARD_CONNECTED_ITEM* aParent )
{
// Save nodes, so they can be cleared later
RN_NODE_PTR start = aEdge->GetSourceNode();
RN_NODE_PTR end = aEdge->GetTargetNode();
start->RemoveParent( aParent );
end->RemoveParent( aParent );
// Connection has to be removed before running RemoveNode(),
// as RN_NODE influences the reference counter
m_links.RemoveConnection( aEdge );
// Remove nodes associated with the edge. It is done in a safe way, there is a check
// if nodes are not used by other edges.
if( m_links.RemoveNode( start ) )
clearNode( start );
if( m_links.RemoveNode( end ) )
clearNode( end );
m_dirty = true;
}
2013-11-25 15:50:03 +00:00
const RN_NODE_PTR& RN_LINKS::AddNode( int aX, int aY )
{
RN_NODE_SET::iterator node;
bool wasNewElement;
2016-09-19 15:25:31 +00:00
std::tie( node, wasNewElement ) = m_nodes.emplace( std::make_shared<RN_NODE>( aX, aY ) );
2013-11-25 15:50:03 +00:00
return *node;
}
bool RN_LINKS::RemoveNode( const RN_NODE_PTR& aNode )
2013-11-25 15:50:03 +00:00
{
if( aNode->GetRefCount() == 0 )
{
2013-11-25 15:50:03 +00:00
m_nodes.erase( aNode );
return true;
}
return false;
2013-11-25 15:50:03 +00:00
}
RN_EDGE_MST_PTR RN_LINKS::AddConnection( const RN_NODE_PTR& aNode1, const RN_NODE_PTR& aNode2,
unsigned int aDistance )
2013-11-25 15:50:03 +00:00
{
2015-06-04 12:54:07 +00:00
assert( aNode1 != aNode2 );
RN_EDGE_MST_PTR edge = std::make_shared<RN_EDGE_MST>( aNode1, aNode2, aDistance );
m_edges.push_back( edge );
2013-11-25 15:50:03 +00:00
return edge;
2013-11-25 15:50:03 +00:00
}
void RN_NET::compute()
{
const RN_LINKS::RN_NODE_SET& boardNodes = m_links.GetNodes();
const RN_LINKS::RN_EDGE_LIST& boardEdges = m_links.GetConnections();
// Special cases do not need complicated algorithms (actually, it does not work well with
// the Delaunay triangulator)
if( boardNodes.size() <= 2 )
2013-11-25 15:50:03 +00:00
{
m_rnEdges.reset( new std::vector<RN_EDGE_MST_PTR>( 0 ) );
2013-11-25 15:50:03 +00:00
// Check if the only possible connection exists
2015-03-10 14:02:27 +00:00
if( boardEdges.size() == 0 && boardNodes.size() == 2 )
2013-11-25 15:50:03 +00:00
{
2016-09-19 15:25:31 +00:00
RN_LINKS::RN_NODE_SET::const_iterator last = ++boardNodes.begin();
2013-11-25 15:50:03 +00:00
// There can be only one possible connection, but it is missing
RN_EDGE_MST_PTR edge = std::make_shared<RN_EDGE_MST>( *boardNodes.begin(), *last );
edge->GetSourceNode()->SetTag( 0 );
edge->GetTargetNode()->SetTag( 1 );
m_rnEdges->push_back( edge );
}
else
{
// Set tags to nodes as connected
for( RN_NODE_PTR node : boardNodes )
node->SetTag( 0 );
2013-11-25 15:50:03 +00:00
}
2013-11-25 15:50:03 +00:00
return;
}
// Move and sort (sorting speeds up) all nodes to a vector for the Delaunay triangulation
std::vector<RN_NODE_PTR> nodes( boardNodes.size() );
std::partial_sort_copy( boardNodes.begin(), boardNodes.end(), nodes.begin(), nodes.end() );
TRIANGULATOR triangulator;
triangulator.CreateDelaunay( nodes.begin(), nodes.end() );
2016-09-13 12:58:33 +00:00
std::unique_ptr<RN_LINKS::RN_EDGE_LIST> triangEdges( triangulator.GetEdges() );
2013-11-25 15:50:03 +00:00
// Compute weight/distance for edges resulting from triangulation
RN_LINKS::RN_EDGE_LIST::iterator eit, eitEnd;
for( eit = (*triangEdges).begin(), eitEnd = (*triangEdges).end(); eit != eitEnd; ++eit )
(*eit)->SetWeight( getDistance( (*eit)->GetSourceNode(), (*eit)->GetTargetNode() ) );
2013-11-25 15:50:03 +00:00
// Add the currently existing connections list to the results of triangulation
std::copy( boardEdges.begin(), boardEdges.end(), std::front_inserter( *triangEdges ) );
// Get the minimal spanning tree
m_rnEdges.reset( kruskalMST( *triangEdges, nodes ) );
}
void RN_NET::clearNode( const RN_NODE_PTR& aNode )
{
if( !m_rnEdges )
return;
std::vector<RN_EDGE_MST_PTR>::iterator newEnd;
2013-11-25 15:50:03 +00:00
// Remove all ratsnest edges for associated with the node
newEnd = std::remove_if( m_rnEdges->begin(), m_rnEdges->end(),
std::bind( isEdgeConnectingNode, _1, std::cref( aNode ) ) );
2013-11-25 15:50:03 +00:00
m_rnEdges->resize( std::distance( m_rnEdges->begin(), newEnd ) );
}
RN_POLY::RN_POLY( const SHAPE_POLY_SET* aParent,
int aSubpolygonIndex,
RN_LINKS& aConnections, const BOX2I& aBBox ) :
m_subpolygonIndex( aSubpolygonIndex ),
m_bbox( aBBox ),
m_parentPolyset( aParent )
2013-11-25 15:50:03 +00:00
{
const VECTOR2I& p = aParent->CVertex( 0, aSubpolygonIndex );
m_node = aConnections.AddNode( p.x, p.y );
2013-11-25 15:50:03 +00:00
2015-06-04 12:54:07 +00:00
// Mark it as not appropriate as a destination of ratsnest edges
2013-11-25 15:50:03 +00:00
// (edges coming out from a polygon vertex look weird)
m_node->SetNoLine( true );
2013-11-25 15:50:03 +00:00
}
bool RN_POLY::HitTest( const RN_NODE_PTR& aNode ) const
{
VECTOR2I p( aNode->GetX(), aNode->GetY() );
2013-11-25 15:50:03 +00:00
return m_parentPolyset->Contains( p, m_subpolygonIndex );
2013-11-25 15:50:03 +00:00
}
void RN_NET::Update()
{
// Add edges resulting from nodes being connected by zones
processZones();
processPads();
2013-11-25 15:50:03 +00:00
compute();
for( RN_EDGE_MST_PTR& edge : *m_rnEdges )
2013-11-25 15:50:03 +00:00
validateEdge( edge );
m_dirty = false;
}
bool RN_NET::AddItem( const D_PAD* aPad )
2013-11-25 15:50:03 +00:00
{
// Ratsnest is not computed for non-copper pads
if( ( aPad->GetLayerSet() & LSET::AllCuMask() ).none() )
return false;
2015-08-04 10:15:47 +00:00
RN_NODE_PTR node = m_links.AddNode( aPad->GetPosition().x, aPad->GetPosition().y );
node->AddParent( aPad );
m_pads[aPad].m_Node = node;
2013-11-25 15:50:03 +00:00
m_dirty = true;
return true;
2013-11-25 15:50:03 +00:00
}
bool RN_NET::AddItem( const VIA* aVia )
2013-11-25 15:50:03 +00:00
{
RN_NODE_PTR node = m_links.AddNode( aVia->GetPosition().x, aVia->GetPosition().y );
node->AddParent( aVia );
m_vias[aVia] = node;
2013-11-25 15:50:03 +00:00
m_dirty = true;
return true;
2013-11-25 15:50:03 +00:00
}
bool RN_NET::AddItem( const TRACK* aTrack )
2013-11-25 15:50:03 +00:00
{
if( aTrack->GetStart() == aTrack->GetEnd() )
return false;
2013-11-25 15:50:03 +00:00
RN_NODE_PTR start = m_links.AddNode( aTrack->GetStart().x, aTrack->GetStart().y );
RN_NODE_PTR end = m_links.AddNode( aTrack->GetEnd().x, aTrack->GetEnd().y );
start->AddParent( aTrack );
end->AddParent( aTrack );
m_tracks[aTrack] = m_links.AddConnection( start, end );
m_dirty = true;
return true;
2013-11-25 15:50:03 +00:00
}
bool RN_NET::AddItem( const ZONE_CONTAINER* aZone )
2013-11-25 15:50:03 +00:00
{
// Prepare a list of polygons (every zone can contain one or more polygons)
const SHAPE_POLY_SET& polySet = aZone->GetFilledPolysList();
2015-06-04 12:54:07 +00:00
// This ensures that we record aZone as added even if it contains no polygons.
(void) m_zones[aZone];
for( int i = 0; i < polySet.OutlineCount(); ++i )
2013-11-25 15:50:03 +00:00
{
const SHAPE_LINE_CHAIN& path = polySet.COutline( i );
2014-03-05 13:57:14 +00:00
RN_POLY poly = RN_POLY( &polySet, i, m_links, path.BBox() );
m_zones[aZone].m_Polygons.push_back( poly );
2013-11-25 15:50:03 +00:00
}
m_dirty = true;
return true;
2013-11-25 15:50:03 +00:00
}
bool RN_NET::RemoveItem( const D_PAD* aPad )
2013-11-25 15:50:03 +00:00
{
PAD_NODE_MAP::iterator it = m_pads.find( aPad );
2013-11-25 15:50:03 +00:00
if( it == m_pads.end() )
return false;
2013-11-25 15:50:03 +00:00
RN_PAD_DATA& pad_data = it->second;
removeNode( pad_data.m_Node, aPad );
2013-11-25 15:50:03 +00:00
for( RN_EDGE_MST_PTR& edge : pad_data.m_Edges )
removeEdge( edge, aPad );
m_pads.erase( aPad );
return true;
2013-11-25 15:50:03 +00:00
}
bool RN_NET::RemoveItem( const VIA* aVia )
2013-11-25 15:50:03 +00:00
{
VIA_NODE_MAP::iterator it = m_vias.find( aVia );
2013-11-25 15:50:03 +00:00
if( it == m_vias.end() )
return false;
2013-11-25 15:50:03 +00:00
removeNode( it->second, aVia );
m_vias.erase( it );
return true;
2013-11-25 15:50:03 +00:00
}
bool RN_NET::RemoveItem( const TRACK* aTrack )
2013-11-25 15:50:03 +00:00
{
TRACK_EDGE_MAP::iterator it = m_tracks.find( aTrack );
2013-11-25 15:50:03 +00:00
if( it == m_tracks.end() )
return false;
removeEdge( it->second, aTrack );
m_tracks.erase( it );
return true;
2013-11-25 15:50:03 +00:00
}
bool RN_NET::RemoveItem( const ZONE_CONTAINER* aZone )
2013-11-25 15:50:03 +00:00
{
ZONE_DATA_MAP::iterator it = m_zones.find( aZone );
if( it == m_zones.end() )
return false;
RN_ZONE_DATA& zoneData = it->second;
2013-11-25 15:50:03 +00:00
// Remove all subpolygons that make the zone
std::deque<RN_POLY>& polygons = zoneData.m_Polygons;
for( RN_POLY& polygon : polygons )
removeNode( polygon.GetNode(), aZone );
polygons.clear();
// Remove all connections added by the zone
std::deque<RN_EDGE_MST_PTR>& edges = zoneData.m_Edges;
for( RN_EDGE_MST_PTR edge : edges )
removeEdge( edge, aZone );
edges.clear();
m_zones.erase( it );
return true;
2013-11-25 15:50:03 +00:00
}
const RN_NODE_PTR RN_NET::GetClosestNode( const RN_NODE_PTR& aNode ) const
{
const RN_LINKS::RN_NODE_SET& nodes = m_links.GetNodes();
RN_LINKS::RN_NODE_SET::const_iterator it, itEnd;
unsigned int minDistance = std::numeric_limits<unsigned int>::max();
RN_NODE_PTR closest;
for( it = nodes.begin(), itEnd = nodes.end(); it != itEnd; ++it )
{
2014-01-31 13:41:15 +00:00
RN_NODE_PTR node = *it;
2013-11-25 15:50:03 +00:00
// Obviously the distance between node and itself is the shortest,
// that's why we have to skip it
2014-01-31 13:41:15 +00:00
if( node != aNode )
2013-11-25 15:50:03 +00:00
{
2014-01-31 13:41:15 +00:00
unsigned int distance = getDistance( node, aNode );
2013-11-25 15:50:03 +00:00
if( distance < minDistance )
{
minDistance = distance;
2014-01-31 13:41:15 +00:00
closest = node;
2013-11-25 15:50:03 +00:00
}
}
}
return closest;
}
2013-12-20 10:54:48 +00:00
const RN_NODE_PTR RN_NET::GetClosestNode( const RN_NODE_PTR& aNode,
const RN_NODE_FILTER& aFilter ) const
2013-11-25 15:50:03 +00:00
{
const RN_LINKS::RN_NODE_SET& nodes = m_links.GetNodes();
RN_LINKS::RN_NODE_SET::const_iterator it, itEnd;
unsigned int minDistance = std::numeric_limits<unsigned int>::max();
RN_NODE_PTR closest;
for( it = nodes.begin(), itEnd = nodes.end(); it != itEnd; ++it )
{
RN_NODE_PTR node = *it;
2013-11-25 15:50:03 +00:00
// Obviously the distance between node and itself is the shortest,
// that's why we have to skip it
if( node != aNode && aFilter( node ) )
2013-11-25 15:50:03 +00:00
{
unsigned int distance = getDistance( node, aNode );
2013-11-25 15:50:03 +00:00
if( distance < minDistance )
{
minDistance = distance;
closest = node;
2013-11-25 15:50:03 +00:00
}
}
}
return closest;
}
std::list<RN_NODE_PTR> RN_NET::GetClosestNodes( const RN_NODE_PTR& aNode, int aNumber ) const
{
std::list<RN_NODE_PTR> closest;
const RN_LINKS::RN_NODE_SET& nodes = m_links.GetNodes();
// Copy nodes
std::copy( nodes.begin(), nodes.end(), std::back_inserter( closest ) );
2013-11-25 15:50:03 +00:00
// Sort by the distance from aNode
closest.sort( std::bind( sortDistance, std::cref( aNode ), _1, _2 ) );
2013-11-25 15:50:03 +00:00
2015-07-30 11:49:36 +00:00
// aNode should not be returned in the results
closest.remove( aNode );
2013-11-25 15:50:03 +00:00
// Trim the result to the asked size
if( aNumber > 0 )
closest.resize( std::min( (size_t)aNumber, nodes.size() ) );
return closest;
}
std::list<RN_NODE_PTR> RN_NET::GetClosestNodes( const RN_NODE_PTR& aNode,
2013-12-20 10:54:48 +00:00
const RN_NODE_FILTER& aFilter, int aNumber ) const
2013-11-25 15:50:03 +00:00
{
std::list<RN_NODE_PTR> closest;
const RN_LINKS::RN_NODE_SET& nodes = m_links.GetNodes();
// Copy filtered nodes
std::copy_if( nodes.begin(), nodes.end(), std::back_inserter( closest ), std::cref( aFilter ) );
2013-11-25 15:50:03 +00:00
// Sort by the distance from aNode
closest.sort( std::bind( sortDistance, std::cref( aNode ), _1, _2 ) );
2013-11-25 15:50:03 +00:00
2015-07-30 11:49:36 +00:00
// aNode should not be returned in the results
closest.remove( aNode );
2013-11-25 15:50:03 +00:00
// Trim the result to the asked size
if( aNumber > 0 )
closest.resize( std::min( static_cast<size_t>( aNumber ), nodes.size() ) );
return closest;
}
void RN_NET::AddSimple( const BOARD_CONNECTED_ITEM* aItem )
{
for( RN_NODE_PTR node : GetNodes( aItem ) )
{
// Block all nodes, so they do not become targets for dynamic ratsnest lines
AddBlockedNode( node );
// Filter out junctions
if( node->GetRefCount() == 1 )
m_simpleNodes.insert( node );
}
}
2013-11-25 15:50:03 +00:00
std::list<RN_NODE_PTR> RN_NET::GetNodes( const BOARD_CONNECTED_ITEM* aItem ) const
{
std::list<RN_NODE_PTR> nodes;
switch( aItem->Type() )
2013-11-25 15:50:03 +00:00
{
case PCB_PAD_T:
{
PAD_NODE_MAP::const_iterator it = m_pads.find( static_cast<const D_PAD*>( aItem ) );
2013-11-25 15:50:03 +00:00
if( it != m_pads.end() )
nodes.push_back( it->second.m_Node );
}
break;
2013-11-25 15:50:03 +00:00
case PCB_VIA_T:
{
VIA_NODE_MAP::const_iterator it = m_vias.find( static_cast<const VIA*>( aItem ) );
2013-11-25 15:50:03 +00:00
if( it != m_vias.end() )
nodes.push_back( it->second );
}
break;
case PCB_TRACE_T:
{
TRACK_EDGE_MAP::const_iterator it = m_tracks.find( static_cast<const TRACK*>( aItem ) );
2014-01-31 13:41:15 +00:00
if( it != m_tracks.end() )
2015-06-04 12:54:07 +00:00
{
nodes.push_back( it->second->GetSourceNode() );
nodes.push_back( it->second->GetTargetNode() );
}
}
break;
case PCB_ZONE_AREA_T:
{
ZONE_DATA_MAP::const_iterator itz = m_zones.find( static_cast<const ZONE_CONTAINER*>( aItem ) );
2015-06-04 12:54:07 +00:00
if( itz != m_zones.end() )
{
const std::deque<RN_POLY>& polys = itz->second.m_Polygons;
2015-06-04 12:54:07 +00:00
for( std::deque<RN_POLY>::const_iterator it = polys.begin(); it != polys.end(); ++it )
nodes.push_back( it->GetNode() );
}
2014-01-31 13:41:15 +00:00
}
break;
default:
break;
2013-11-25 15:50:03 +00:00
}
return nodes;
}
2015-03-10 08:36:04 +00:00
void RN_NET::GetAllItems( std::list<BOARD_CONNECTED_ITEM*>& aOutput, RN_ITEM_TYPE aType ) const
{
if( aType & RN_PADS )
{
2016-09-13 12:58:33 +00:00
for( auto it : m_pads )
aOutput.push_back( const_cast<D_PAD*>( it.first ) );
2015-03-10 08:36:04 +00:00
}
if( aType & RN_VIAS )
{
2016-09-13 12:58:33 +00:00
for( auto it : m_vias )
aOutput.push_back( const_cast<VIA*>( it.first ) );
2015-03-10 08:36:04 +00:00
}
if( aType & RN_TRACKS )
{
2016-09-13 12:58:33 +00:00
for( auto it : m_tracks )
aOutput.push_back( const_cast<TRACK*>( it.first ) );
2015-03-10 08:36:04 +00:00
}
if( aType & RN_ZONES )
{
2016-09-13 12:58:33 +00:00
for( auto it : m_zones )
aOutput.push_back( const_cast<ZONE_CONTAINER*>( it.first ) );
2015-03-10 08:36:04 +00:00
}
}
void RN_NET::ClearSimple()
{
for( const RN_NODE_PTR& node : m_blockedNodes )
node->SetNoLine( false );
m_blockedNodes.clear();
m_simpleNodes.clear();
}
void RN_NET::GetConnectedItems( const BOARD_CONNECTED_ITEM* aItem,
std::list<BOARD_CONNECTED_ITEM*>& aOutput,
2015-03-10 08:36:04 +00:00
RN_ITEM_TYPE aTypes ) const
{
std::list<RN_NODE_PTR> nodes = GetNodes( aItem );
assert( !nodes.empty() );
int tag = nodes.front()->GetTag();
assert( tag >= 0 );
if( aTypes & RN_PADS )
{
for( PAD_NODE_MAP::const_iterator it = m_pads.begin(); it != m_pads.end(); ++it )
{
if( it->second.m_Node->GetTag() == tag )
aOutput.push_back( const_cast<D_PAD*>( it->first ) );
}
}
if( aTypes & RN_VIAS )
{
for( VIA_NODE_MAP::const_iterator it = m_vias.begin(); it != m_vias.end(); ++it )
{
if( it->second->GetTag() == tag )
aOutput.push_back( const_cast<VIA*>( it->first ) );
}
}
if( aTypes & RN_TRACKS )
{
for( TRACK_EDGE_MAP::const_iterator it = m_tracks.begin(); it != m_tracks.end(); ++it )
{
if( it->second->GetTag() == tag )
aOutput.push_back( const_cast<TRACK*>( it->first ) );
}
}
if( aTypes & RN_ZONES )
{
for( ZONE_DATA_MAP::const_iterator it = m_zones.begin(); it != m_zones.end(); ++it )
{
for( const RN_EDGE_MST_PTR& edge : it->second.m_Edges )
{
if( edge->GetTag() == tag )
{
aOutput.push_back( const_cast<ZONE_CONTAINER*>( it->first ) );
break;
}
}
}
}
}
void RN_DATA::AddSimple( const BOARD_ITEM* aItem )
2013-11-25 15:50:03 +00:00
{
if( aItem->IsConnected() )
{
const BOARD_CONNECTED_ITEM* item = static_cast<const BOARD_CONNECTED_ITEM*>( aItem );
int net = item->GetNetCode();
2013-11-25 15:50:03 +00:00
// Do not process orphaned & unconnected items
if( net <= NETINFO_LIST::UNCONNECTED )
return;
2013-11-25 15:50:03 +00:00
m_nets[net].AddSimple( item );
}
else if( aItem->Type() == PCB_MODULE_T )
{
const MODULE* module = static_cast<const MODULE*>( aItem );
2013-11-25 15:50:03 +00:00
for( const D_PAD* pad = module->Pads().GetFirst(); pad; pad = pad->Next() )
AddSimple( pad );
}
2013-11-25 15:50:03 +00:00
}
void RN_DATA::AddBlocked( const BOARD_ITEM* aItem )
{
if( aItem->IsConnected() )
{
const BOARD_CONNECTED_ITEM* item = static_cast<const BOARD_CONNECTED_ITEM*>( aItem );
int net = item->GetNetCode();
// Do not process orphaned & unconnected items
if( net <= NETINFO_LIST::UNCONNECTED )
return;
// Block all nodes belonging to the item
for( RN_NODE_PTR node : m_nets[net].GetNodes( item ) )
m_nets[net].AddBlockedNode( node );
}
else if( aItem->Type() == PCB_MODULE_T )
{
const MODULE* module = static_cast<const MODULE*>( aItem );
for( const D_PAD* pad = module->Pads().GetFirst(); pad; pad = pad->Next() )
AddBlocked( pad );
}
}
void RN_DATA::GetConnectedItems( const BOARD_CONNECTED_ITEM* aItem,
std::list<BOARD_CONNECTED_ITEM*>& aOutput,
2015-03-10 08:36:04 +00:00
RN_ITEM_TYPE aTypes ) const
{
int net = aItem->GetNetCode();
if( net < 1 )
return;
assert( net < (int) m_nets.size() );
m_nets[net].GetConnectedItems( aItem, aOutput, aTypes );
}
2015-03-10 08:36:04 +00:00
void RN_DATA::GetNetItems( int aNetCode, std::list<BOARD_CONNECTED_ITEM*>& aOutput,
RN_ITEM_TYPE aTypes ) const
{
if( aNetCode < 1 )
return;
assert( aNetCode < (int) m_nets.size() );
m_nets[aNetCode].GetAllItems( aOutput, aTypes );
}
bool RN_DATA::AreConnected( const BOARD_CONNECTED_ITEM* aItem, const BOARD_CONNECTED_ITEM* aOther )
{
int net1 = aItem->GetNetCode();
int net2 = aOther->GetNetCode();
if( net1 < 1 || net2 < 1 || net1 != net2 )
return false;
assert( net1 < (int) m_nets.size() && net2 < (int) m_nets.size() );
// net1 == net2
std::list<RN_NODE_PTR> items1 = m_nets[net1].GetNodes( aItem );
std::list<RN_NODE_PTR> items2 = m_nets[net1].GetNodes( aOther );
assert( !items1.empty() && !items2.empty() );
return ( items1.front()->GetTag() == items2.front()->GetTag() );
}
int RN_DATA::GetUnconnectedCount() const
{
int count = 0;
for( unsigned i = 0; i < m_nets.size(); ++i )
{
const std::vector<RN_EDGE_MST_PTR>* unconnected = m_nets[i].GetUnconnected();
if( unconnected )
count += unconnected->size();
}
return count;
}
2013-11-25 15:50:03 +00:00
void RN_NET::processZones()
{
for( ZONE_DATA_MAP::iterator it = m_zones.begin(); it != m_zones.end(); ++it )
2013-11-25 15:50:03 +00:00
{
const ZONE_CONTAINER* zone = it->first;
RN_ZONE_DATA& zoneData = it->second;
// Reset existing connections
for( RN_EDGE_MST_PTR edge : zoneData.m_Edges )
2013-11-25 15:50:03 +00:00
m_links.RemoveConnection( edge );
zoneData.m_Edges.clear();
LSET layers = zone->GetLayerSet();
2013-11-25 15:50:03 +00:00
// Compute new connections
2015-06-04 12:54:07 +00:00
RN_LINKS::RN_NODE_SET candidates = m_links.GetNodes();
2016-09-19 15:25:31 +00:00
RN_LINKS::RN_NODE_SET::const_iterator point, pointEnd;
2013-11-25 15:50:03 +00:00
// Sorting by area should speed up the processing, as smaller polygons are computed
// faster and may reduce the number of points for further checks
std::sort( zoneData.m_Polygons.begin(), zoneData.m_Polygons.end(), sortArea );
2013-11-25 15:50:03 +00:00
for( std::deque<RN_POLY>::iterator poly = zoneData.m_Polygons.begin(),
polyEnd = zoneData.m_Polygons.end(); poly != polyEnd; ++poly )
2013-11-25 15:50:03 +00:00
{
2015-06-04 12:54:07 +00:00
const RN_NODE_PTR& node = poly->GetNode();
2013-11-25 15:50:03 +00:00
point = candidates.begin();
pointEnd = candidates.end();
while( point != pointEnd )
{
if( *point != node && ( (*point)->GetLayers() & layers ).any()
&& poly->HitTest( *point ) )
2013-11-25 15:50:03 +00:00
{
//(*point)->AddParent( zone ); // do not assign parent for helper links
2015-06-04 12:54:06 +00:00
RN_EDGE_MST_PTR connection = m_links.AddConnection( node, *point );
zoneData.m_Edges.push_back( connection );
2013-11-25 15:50:03 +00:00
// This point already belongs to a polygon, we do not need to check it anymore
point = candidates.erase( point );
pointEnd = candidates.end();
}
else
{
++point;
}
}
}
}
}
void RN_NET::processPads()
{
for( PAD_NODE_MAP::iterator it = m_pads.begin(); it != m_pads.end(); ++it )
{
const D_PAD* pad = it->first;
RN_NODE_PTR node = it->second.m_Node;
std::deque<RN_EDGE_MST_PTR>& edges = it->second.m_Edges;
// Reset existing connections
for( RN_EDGE_MST_PTR edge : edges )
m_links.RemoveConnection( edge );
LSET layers = pad->GetLayerSet();
const RN_LINKS::RN_NODE_SET& candidates = m_links.GetNodes();
2016-09-19 15:25:31 +00:00
RN_LINKS::RN_NODE_SET::const_iterator point, pointEnd;
point = candidates.begin();
pointEnd = candidates.end();
while( point != pointEnd )
{
if( *point != node && ( (*point)->GetLayers() & layers ).any() &&
pad->HitTest( wxPoint( (*point)->GetX(), (*point)->GetY() ) ) )
{
//(*point)->AddParent( pad ); // do not assign parent for helper links
RN_EDGE_MST_PTR connection = m_links.AddConnection( node, *point );
edges.push_back( connection );
}
++point;
}
}
}
bool RN_DATA::Add( const BOARD_ITEM* aItem )
2013-11-25 15:50:03 +00:00
{
int net = NETINFO_LIST::ORPHANED;
if( aItem->IsConnected() )
{
net = static_cast<const BOARD_CONNECTED_ITEM*>( aItem )->GetNetCode();
}
else if( aItem->Type() == PCB_MODULE_T )
{
const MODULE* module = static_cast<const MODULE*>( aItem );
for( const D_PAD* pad = module->Pads().GetFirst(); pad; pad = pad->Next() )
{
net = pad->GetNetCode();
2014-06-03 14:08:23 +00:00
// Do not process orphaned items
if( net <= NETINFO_LIST::ORPHANED )
continue;
Add( pad );
}
return true;
}
else if( aItem->Type() == PCB_NETINFO_T )
{
int netCount = m_board->GetNetCount();
if( (unsigned) netCount > m_nets.size() )
m_nets.resize( netCount );
return true;
}
else
{
return false;
}
2013-11-25 15:50:03 +00:00
if( net < 0 )
return false;
// Autoresize is necessary e.g. for module editor
if( net >= (int) m_nets.size() )
m_nets.resize( net + 1 );
2013-11-25 15:50:03 +00:00
switch( aItem->Type() )
{
case PCB_PAD_T:
return m_nets[net].AddItem( static_cast<const D_PAD*>( aItem ) );
break;
2013-11-25 15:50:03 +00:00
case PCB_TRACE_T:
return m_nets[net].AddItem( static_cast<const TRACK*>( aItem ) );
break;
2013-11-25 15:50:03 +00:00
case PCB_VIA_T:
return m_nets[net].AddItem( static_cast<const VIA*>( aItem ) );
break;
2013-11-25 15:50:03 +00:00
case PCB_ZONE_AREA_T:
return m_nets[net].AddItem( static_cast<const ZONE_CONTAINER*>( aItem ) );
break;
2013-11-25 15:50:03 +00:00
default:
break;
}
return false;
2013-11-25 15:50:03 +00:00
}
bool RN_DATA::Remove( const BOARD_ITEM* aItem )
2013-11-25 15:50:03 +00:00
{
int net = NETINFO_LIST::ORPHANED;
2013-11-25 15:50:03 +00:00
if( aItem->IsConnected() )
{
net = static_cast<const BOARD_CONNECTED_ITEM*>( aItem )->GetNetCode();
}
else if( aItem->Type() == PCB_MODULE_T )
{
const MODULE* module = static_cast<const MODULE*>( aItem );
for( const D_PAD* pad = module->Pads().GetFirst(); pad; pad = pad->Next() )
2013-11-25 15:50:03 +00:00
{
net = pad->GetNetCode();
2014-06-03 14:08:23 +00:00
// Do not process orphaned items
if( net <= NETINFO_LIST::ORPHANED )
continue;
Remove( pad );
2013-11-25 15:50:03 +00:00
}
return true;
2013-11-25 15:50:03 +00:00
}
else
{
return false;
}
if( net < 0 )
return false;
// Autoresize is necessary e.g. for module editor
if( net >= (int) m_nets.size() )
{
m_nets.resize( net + 1 );
return false; // if it was resized, then surely the item had not been added before
}
switch( aItem->Type() )
{
case PCB_PAD_T:
return m_nets[net].RemoveItem( static_cast<const D_PAD*>( aItem ) );
break;
case PCB_TRACE_T:
return m_nets[net].RemoveItem( static_cast<const TRACK*>( aItem ) );
break;
case PCB_VIA_T:
return m_nets[net].RemoveItem( static_cast<const VIA*>( aItem ) );
break;
case PCB_ZONE_AREA_T:
return m_nets[net].RemoveItem( static_cast<const ZONE_CONTAINER*>( aItem ) );
break;
default:
break;
}
return false;
}
bool RN_DATA::Update( const BOARD_ITEM* aItem )
{
if( Remove( aItem ) )
{
bool res = Add( aItem );
assert( res );
(void) res;
return true;
}
return false;
2013-11-25 15:50:03 +00:00
}
void RN_DATA::ProcessBoard()
{
2015-06-04 12:54:07 +00:00
int netCount = m_board->GetNetCount();
2013-11-25 15:50:03 +00:00
m_nets.clear();
2015-06-04 12:54:07 +00:00
m_nets.resize( netCount );
int netCode;
2013-11-25 15:50:03 +00:00
// Iterate over all items that may need to be connected
for( MODULE* module = m_board->m_Modules; module; module = module->Next() )
{
for( D_PAD* pad = module->Pads().GetFirst(); pad; pad = pad->Next() )
{
2014-03-06 08:42:16 +00:00
netCode = pad->GetNetCode();
2015-06-04 12:54:07 +00:00
assert( netCode >= 0 && netCode < netCount );
if( netCode >= 0 && netCode < netCount )
m_nets[netCode].AddItem( pad );
}
2013-11-25 15:50:03 +00:00
}
for( TRACK* track = m_board->m_Track; track; track = track->Next() )
{
2014-03-06 08:42:16 +00:00
netCode = track->GetNetCode();
2015-06-04 12:54:07 +00:00
assert( netCode >= 0 && netCode < netCount );
if( netCode >= 0 && netCode < netCount )
{
if( track->Type() == PCB_VIA_T )
m_nets[netCode].AddItem( static_cast<VIA*>( track ) );
else if( track->Type() == PCB_TRACE_T )
m_nets[netCode].AddItem( track );
}
2013-11-25 15:50:03 +00:00
}
for( int i = 0; i < m_board->GetAreaCount(); ++i )
{
ZONE_CONTAINER* zone = m_board->GetArea( i );
2014-03-06 09:43:40 +00:00
2014-03-06 08:42:16 +00:00
netCode = zone->GetNetCode();
2015-06-04 12:54:07 +00:00
assert( netCode >= 0 && netCode < netCount );
if( netCode >= 0 && netCode < netCount )
m_nets[netCode].AddItem( zone );
2013-11-25 15:50:03 +00:00
}
Recalculate();
2013-11-25 15:50:03 +00:00
}
void RN_DATA::Recalculate( int aNet )
{
2014-06-03 14:08:23 +00:00
unsigned int netCount = m_board->GetNetCount();
if( aNet <= 0 && netCount > 1 ) // Recompute everything
2013-11-25 15:50:03 +00:00
{
2015-06-04 12:54:06 +00:00
#ifdef PROFILE
PROF_COUNTER totalRealTime;
2015-06-04 12:54:06 +00:00
#endif
2014-06-03 14:08:23 +00:00
unsigned int i;
#ifdef USE_OPENMP
2014-03-05 13:57:14 +00:00
#pragma omp parallel shared(netCount) private(i)
2013-11-25 15:50:03 +00:00
{
2014-03-05 13:57:14 +00:00
#pragma omp for schedule(guided, 1)
#else /* USE_OPENMP */
{
#endif
2014-03-05 13:57:14 +00:00
// Start with net number 1, as 0 stands for not connected
for( i = 1; i < netCount; ++i )
{
if( m_nets[i].IsDirty() )
updateNet( i );
}
2014-03-05 13:57:14 +00:00
} /* end of parallel section */
2015-06-04 12:54:06 +00:00
#ifdef PROFILE
totalRealTime.Stop();
wxLogDebug( "Recalculate all nets: %.1f ms", totalRealTime.msecs() );
2015-06-04 12:54:06 +00:00
#endif /* PROFILE */
2013-11-25 15:50:03 +00:00
}
else if( aNet > 0 ) // Recompute only specific net
2013-11-25 15:50:03 +00:00
{
updateNet( aNet );
}
}
void RN_DATA::updateNet( int aNetCode )
{
assert( aNetCode < (int) m_nets.size() );
2014-02-05 10:55:04 +00:00
if( aNetCode < 1 || aNetCode > (int) m_nets.size() )
return;
m_nets[aNetCode].ClearSimple();
m_nets[aNetCode].Update();
}