Генетический алгоритм глобальной трассировки (стр. 1 из 3)

О.Б. Лебедев

1 Введение

Основной целью задачи глобальной трассировки является равномерное и целесообразное распределение ресурсов коммутационного поля для создания благоприятных условий для последующей детальной трассировки.

Большинство алгоритмов, глобальной трассировки осуществляют последовательное построение соединений на укрупненной модели КП (волновые, лучевые, базирующиеся на построении деревьев Штейнера). [1,2,3,4,5,6,7]. Хотя на каждом шаге для каждого текущего состояния среды алгоритмы дают неплохие результаты, «камнем преткновения» является последовательность трассируемых соединений. Цепи, проложенные раньше ,могут блокировать цепи ,прокладываемые позже.

Другим недостатком является то, что большинство алгоритмов используют критерии в большей степени учитывающие параметры соединений (например: общая длина) и в меньшей степени параметры коммутационного поля, что не совсем согласуется с главной целью глобальной трассировки.

Исходя из этих соображений, в работе используется комбинаторный подход, основанный на методах генетической адаптации, при котором в один и тот же момент времени рассматриваются все соединения, а критерий учитывает распределение ресурсов КП.

При разработке генетических процедур основное влияние уделялось разработке с учетом знаний о предметной области методов кодирования решений, модификации генетических операторов и организации эволюционного процесса.

2. Проблемная формулировка, термины и обозначения

Для решения задачи глобальной трассировки используется графовая модель G = (X, U). Коммутационное поле разбивается на области. Вершины графа x_iÎ Х соответствуют областям на КП. Если области соседние, то вершины x_i и x_к, соответствующие этим областям, связываются ребром u_j. Пусть задано множество цепей Т={t_i| i = 1, 2, ...}. Для каждой цепи определяется множество областей, в которых существуют контакты, связываемые этой цепью. На графовой модели G, множеству областей, связываемых цепью x_i, соответствует подмножество вершин Х_iÌ Х. Для каждого ребра u_i,связывающего вершины, задается вес a_j, равный пропускной способности между областями, соответствующими вершинам x_i и x_k. Будем считать, что граф G метризирован, то есть каждая вершина имеет координаты. Координаты вершины принимаются равными координатам центра соответствующей области. Если области имеют один и тот же размер, то граф G представляет собой ортогональную решетку.

Далее для каждой цепи t_i на множестве вершин Х_i графа G строится минимальное связывающее дерево D_i с помощью алгоритма Прима D_i={r_k | k=1,2, … ,n_k}, где r_k – ребро минимального связывающего дерева.

Для каждого ребра r_kÎD_i формируется набор V_k вариантов v_ik маршрутов, связывающих на графе G соответствующие вершины. Формирование возможных маршрутов осуществляется следующим образом. Для ребра r_kÎD_i, связывающего x_nÎG и x_mÎG, определяется множество вершин Х_kÌX, смежных вершинам х_nи х_m ребра r_k,. Через множество вершин Х_к , а так же через вершины x_n и x_m проводятся новые вертикальные и горизонтальные линии.Отметим, что эти линии проходят по ребрам ортогонального графа G. В узлах пересечения этих линий лежат некоторые вершины х_iÎ Х. Эти вершины являются узловыми для формирования вариантов. Будем считать, что варианты маршрутов проходят по тем ребрам графа G, которые лежат на этих линиях, рис 1.

Например: сформируем набор вариантов для ребра г_к, связывающего вершины х_i и х_j. Пронумеруем узлы пересечения вертикальных и горизонтальных линий. Вершина х_i лежит в узле 3, а вершина х_j в узле 10. Для данного ребра г_k существует 10 вариантов прохождения маршрута V_k = {v_k1, v_k2, v_k3, v_k4, v_k5, v_k6, v_k7, v_k8, v_k9, v_k10}

v_k1 = {3, 2, 5, 8, 11, 10}; v_k2 = {3, 6, 5, 8, 11, 10}; v_k3 = {3, 6, 9, 8, 11, 10}; v_k4 = {3, 6, 9, 12, 11, 10}; v_k5 = {3, 2, 1, 4, 7, 10}; v_k6 = {3, 2, 5, 4, 7, 10}; v_k7 = {3, 6, 5, 4, 7, 10}; v_k8 = {3, 2, 5, 8, 7, 10}; v_k9 = {3, 6, 5, 8, 7, 10}; v_k10 = {3, 6, 9, 8, 7, 10}

Такой способ обеспечивает максимальное совпадение вариантов реализации ребра r_kс вариантами ребер, смежных ребру r_k. Формирование варианта осуществляется из следующих соображений:

вариант формируется таким образом, чтобы он был минимальной длины.

варианты формируются, чтобы обеспечивалось максимально возможное совпадение вариантов различных ребер одной цепи друг с другом.

Пусть имеется некоторое решение задачи глобальной трассировки, заключающееся в том, что для всех ребер г_k всех МСД-цепей, выбраны варианты их реализации. Введем некоторые обозначения. Обозначим через b_i число цепей, проходящих по ребру u_i графа G. Для каждого ребра u_i графа G введем параметр с_i, равный с_i = a_i - b_i.

Найдем в графе G ребро, у которого с_i имеет минимальное значение, и обозначим его через с_min, то есть c_min®"i[c_min£c_i].

Для нашей задачи цель оптимизации – максимизация параметра с_min. При перезагрузках параметр с_min может принимать отрицательное значение. Для оценки качества удобнее использовать критерий, который имеет только положительные значения. Обозначим через a_m – максимально возможное значение показателя a_i среди всех ребер u_i графа G, т.е. a_m® ("_i)[ a_m ³a_i]. В силу этого величина a_m– c_min никогда не может быть отрицательной, т.е. a_m– c_m³ 0. Отметим, что a_m является константой.

В качестве фитнесса (оценки качества) индивидуальности будет использоваться критерий

Цель оптимизации – максимизация F₁.По своей сути максимизация F₁ совпадает с максимизацией c_min.

В качестве другого критерия оптимизации можно использовать другую характеристику. Обозначим через g число ребер графа G для которых c_i имеет отрицательное значение. Очевидно, что целью генетического поиска должен быть поиск решения (индивидуальности) у которого g имеет минимальное значение.

Пусть m – число всех ребер графа G. Введем фитнесс F₂=m - g.

Цель оптимизации в таком случае – максимизация F₂

Если при оптимизации по критерию F₁ не удается найти такого решения при котором у всех ребер u_iÎG, соответствующие им с_i имеют положительные значения, то оптимизация по критерию F₂ минимизирует число таких ребер. Это делается из следующих соображений: чтобы обеспечить 100% реализацию (разводки) соединений, потребуется меньшее число деформаций областей или границ.

Обозначим через j число ребер всех цепей, для которых выбранные и реализованные варианты соединений включают ребра u_i графа G с отрицательными значениями с_i.

Пусть N_d – общее число ребер всех максимальных связывающих деревьев.

Введем третий критерий оптимизации: F₃=N_d - j

Цель оптимизации заключается в максимизации F₃. Это обеспечивает минимальное число соединений, требующих перетрассировки, т.е. оптимизация по F₃ уменьшает число перетрассируемых соединений. Если при оптимизации по F₁, имеются ребра с отрицательными с_i, то в зависимости от метода дотрассировки, деформации областей или границ, используется критерий F₂ или F₃.

3. Разработка генетических процедур для задачи глобальной трассировки

3.1. Кодирование и декодирование хромосом

Для решения задачи глобальной трассировки используются генетические методы оптимизации. Представим решение в виде хромосомы. Кодирование осуществляется следующим образом. Хромосома состоит из генов. Количество генов в хромосоме H_i равно количеству ребер минимальных связывающих деревьев для всех цепей, расположенных на КП. Значением гена является номер варианта из заданного набора вариантов маршрутов, связывающих на графе G соответствующие вершины.

Например, дано КП, на котором расположено множество цепей Т = {t₁, t₂, t₃} , рис. 2.

Так как мы используем графовую модель, то КП можно представить соответственно рис. 3.

Для цепи t₁ множество связуемых вершин – Х₁ = {x₆, x₂, х₁₁}. Для цепи t₂ множество связуемых вершин – Х₂ = {x₇, x₉, x₁₃ }. Для цепи t₃ множество связуемых вершин – Х₃ = {x₅, x₁₄}. С помощью алгоритма Прима для каждой цепи строится минимальное связывающее дерево МСД. Для каждой из цепей это выглядит так:рис. 4.

После этого для каждого ребра r_ij МСД формируется набор вариантов маршрутов, связывающих на графе G соответствующие вершины.

Ребро г₁₁ (то есть первое ребро МСД для цепи 1) имеет два варианта прохождения маршрута r₁₁ = {v₁₁₁, v₁₁₂}: v₁₁₁={x₆,x₁,x₂}, v₁₁₂={x₆,x₇,x₂}.

Ребро r₁₂ (второе ребро цепи 1) имеет один вариант V₁₂₁={x₆,x₁₁}

Ребро г₂₁ (то есть первое ребро МСД для цепи 2) имеет два варианта прохождения маршрута r₂₁={v₂₁₁,v₂₁₂} : v₂₁₁={x₇,x₈,x₁₃} v₂₁₂={x₇,x₁₂,x₁₃}.

Ребро r₂₂ имеет два варианта v₂₂₁={x₁₃,x₈,x₉}, v₂₂₂={x₁₃,x₁₄,x₉}.