Устройство для решения нелинейных краевых задач

СОЮЗ СОВЕТСКИХСОЦИАЛИСТИЧЕСКИХРЕСПУБЛИК А 06 Е 15/20 ОПИСАН ЗОБРЕТЕНИЯ АВ ОМУ СВИДЕТЕЛЬСТ(56) Авторское свидетельство СССР М 383069, кл. б 06 0 7/48, 1970.Авторское свидетельство СССР М 918951, кл. 0 06 Р 15/328, 1982,(54) УСТРОЙСТВО ДЛЯ РЕШЕНИЯ НЕЛИНЕЙНЫХ КРАЕВЫХ ЗАДАЧ(57) Изобретение относится к аналоговой вычислительной технике и может быть применено для решения краевыхзадач, описываемых дифференциальными уравнениями в частных производных с функциональными и нелинейными зависимостями коэффициентов, методами дискретного моделирования. Целью изобретения является Голенкова, Б.Пинигин и ОО ГОСУДАРСТВЕННЫЙ КОМИТЕТПО ИЗОБРЕТЕНИЯМ И ОТКРЫТИЯМПРИ ГКНТ СССР повышение быстродействия устройства, Устройство содержит регистр 1 конфигурации, коммутатор 2, матрицу блоков 3 решения уравнений системы, внешнюю злектронновычислительную машину (ЭВМ) 4, Каждый блок 3 содержит узел 5 вычисления системы конечно-разностных уравнений, блок 6 коммутации, регистр 7 данных, дешифратор 8 адреса, шифратор 9, аналого-цифровой преобразователь 10, дешифратор 11 команд, шифратор команд, распределитель 12 импульсов, первый контроллер 13 шины, управляющие входы 14, информационные выходы-входы 15, третий контроллер 16 шины (с магистралью Х, т), второй контроллер 17 шины, коммутатор 18 памяти, блок 19 памяти, переключатель 20 локальных общих шин, первую общую шину 21, вторую общую шину 22, первую, вторую и третью локальные магистрали 231, 232 и 23 з, входы-выходы 24 первой группы для связи граничныхЯ1683028 20 30 35 4550 55 где а", ау, а - соответствующие коэффициенты связи между узловыми точками по координатам Х и У, сеточного шаблона соответственно1;,=, - в случае граничных условийрода;= 1- в случае граничных условийрода,Б =о,Ч 1, - в случае граничных условий И рода;о, в .с,Ц - в случае моделирования временной производной по методу Либмана;о,;,- диагональный элемент, преда в случае граничных условийрода или временной производной и диагонального "реобладания а.Блок-схема программы организации вычислительного процесса в устройстве содержит блоки занесения исходных данных в реша;ощие блоки, установки конфигурации области, масштабирования, занесения информации в матрицу узловых процессоров, съема решения и демасштабирования,уточнения решения, учета нелинейности,выдачи решения из решающих блоков вЭВМ,Иэ ЭВМ 4 через общие шины 21 и 22 иблоки связи с ЭВМ в блоки 19 памяти матрицы блоков 3 решения уравнений системы записывается программа работы устройства в исходные данные системы конечно-раэностных уравнений АХ = В. Для каждого коэффициента а", аУ, а, Ь, О в блоке 19 памяти отведена определенная страница памяти, длина которой равна 4 й байтов, где Й -число узлов процессоров в матрице реша,ащего блока 3, Процесс передачи данных и программ из ЭВМ 4 в блок 19 памяти осуществляется через контроллер 17 шины представляющий собой скоростной канал связи, состоящий из двух однотипных интерфейсных узлов(не показаны). Последние имеют независимое управление со стороны входа 14 управления и Э ВМ 4 соответственно в процессе пословного или блочного обмена информацией между ними,Процесс обмена информацией между ЭВМ 4 и входами 14 управления строится по принципу запрос-ответ, т.е. в любом обмене информацией существует процессор- инициатор обмена и процессор-абонент.При передаче данных ЭВМ 4 производит загрузку адресуемых регистров контроллера 17 шины со стороны общей шины 21, причем в регистры начального и конечного адресов контроллера засылаются коды начального и конечного адресов Ак передаваемого блока данных, В регистр команд и состояния контроллера заносится код команды, разрешения прерывания и разряд пуска, Выборка регистров осущесте-,яется кодом адреса на общей шине 21.Программа работы устройства состоит из сменных подпрограмм, реализующих определенные блоки алгоритма работы устройства, загрузка которых в блоки 19 памяти блоков 3 производится аналогично.Организация взаимодействия между ЭВМ 4 и блоками 3 осуществляется путем передачи адреса подпрограммы, вызываемой в блоках 3.ЭВМ 4 последовательно передает в контроллер 17 шины начальный адрес подпрограммы и код команды записи одиночного слова, Аналогично ЭВМ 4 осуществляет запуск остальных блоков 3.Для установки конфигурации области осуществляется программная настройка блоков 3 устройства для решения конкретной краевой задачи, Из ЭВМ 4 в регистр 1 конфигурации выдается код тороидальной базовой области. Из регистра 1 этот код выдается на коммутатор 2, который осуществляет коммутацию границ матрицы узлов 5 блоков 3 по координатам Х, У, образуя необходимую конфигурацию моделирующей базовой области в плоскости Х, У. По координате Е координатчые проводимости 9 узлов 5 блоков 3 через входы-выходы 25 соединены последовательно в кольцо. Кроме того, осуществляетСя программная настройка контроллера 16 шины. Контроллер 16 шины представляет собой программно настраиваемый блок, предназначенный для сопряжения данного блока 3 с общей шиной 22, которая содержит линию связл, объединяющую между собой блоки контроллера 16 шины всех блоков 3, и включает в себя набор линий управления, Контроллер 16 шины обеспечивает выдачу заявок на общую шину 22 для чтения информации из соседних блоков 3, прием и передачу данных с общей шины 22 на шины данных локальных общих шин 23, чтение и выдачу данных из блока 19 памяти своего узла 5 на шины данных по заявкам блоков 3, Для обеспечения связи блока 3 с другими в составе контроллера 16 шины имеется блок памяти реконфигурации, предназначенный для хранения значений номеров соседних блоков 3 в области моделирования краевой задачи, расположенных слева и справа по координате Х, а также снизу и сверху по координате У соответственно. Кроме того. блок памяти реконфигурации обеспечивает хранение номеров соседних блоков 3 по отношению к резервируемому решающему блоку 3 по координате Е, 168302850 при р+ц нечетном, и ЭВМ 4 последовательно передает в блоки 3 через контроллеры 17 шины начальный адрес подпрогоаммы масштабирования. Процесс масштабирования заключается в пересчете исходных данных (параметров) исходного уравнения к параметрам узлов, для чего определяются максимальные значения коэффициентов (амакс) и правой части (Ь) )акс), определяются масштабные коэффициенты по проводимости Мя по току М 1 и по напряжению МО осуществляется расчет кодов параметров Ох, Оу, Оъ Ос и токов 1) узлов 5. Определение максимальных значений и расчет кодов параметров узлов 5 осуществляется по информационным выходам-входам 15 под управлением по входам 14 в соответствии с подпрограммой масштабирования.Команды по выходам-входам 15 различаются не только количеством адресуемых операндов, одновременно участвующих в операции, но и видом самих операций. Операнды, участвующие в операциях, могут быть как скалярными, так и векторными величинами, Скалярные величины могут быть представлены в виде одиночных 32-х или 64-х разрядных слов, размещаемых в ячейках блока 19 памяти. Векторные величины представляют в виде совокупности 32-х или 64-х разрядных слов, размещаемых в последовательных ячейках блоков 19 на целочисленных границах слов. Размерность векторной величины может быть до Й слов,В системе команд по выходам-входам , 15 предусмотрено четыре формата команд, первый и второй форматы из которых используются для выполнения операций над двумя скалярными величинами и операции одной скалярной величины с элементами векторной величины. Значение скалярной величины в регистрах общего назначения, значения элементов векторных величин, размещенных в последовательных ячейках памяти, определяются совокупностью значений адресов страниц и адресов точек. Остальные форматы используются для выполнения операций над элементами двух векторных величин, размещенных в блоке 19 памяти, начиная с адресов, заданных значениями для первого и второго операндов. Результат операции в виде третьей векторной величины размещается в блоке 19 памяти. С учетом особенностей реализации цифровых итерационных алгоритмов в систему команд по выходам-входам 15 введены команды для выполнения групповых операций "Поиск максимума", "Умножение на константу", "Умножение с накоплением", "Деление", "Вычисление четырехточечного 5 10 15 20 25 30 35 40 шаблона", "Вычисление шеститочечного шаблона".В процессе работы устройства осуществляется поиск максимума коэффициентов а", а", ас соответствующих страниц памяти, определение максимального значения коэффициента амакс и вычисление. масштабного коэффициента проводимости Мя, По управляющим входам 14 и выходам-входам 15 осуществляется симметрирование коэффициентов матрицы А и нахождение кодов проводимостей Оу, Ог, От.Далее выполняется подпрограмма определения максимального значения правой части Ьасс, определения масштабных коэффициентов потоку М 1 и по напряжению М).). Блок 19 памяти по заданному сигналу выставляет информацию на локальные шины 23. Информация сигналом записи по выходу 32 с распределителя 12 импульсов переписывается в регистр 7 данных и через шифратор 9 поступает на информационные входы 26 узлов 5. Младшие разряды с регистра адреса контроллера 13 шины поступают на дешифратор 8, сигналы с выхода которого разрешают запись принятой информации. Дешифратор 11 введен для формирования "1" в младшем разряде кода проводимостей при переходе от 32-разрядной информации к усеченной разрядности кодоуправляемых проводимостей для исключения разрывов в матрице проводимостей,После получения сообщений от всех блоков 3 ЭВМ 4 последовательно передает в эти блоки начальный адрес подпрограммы уточнения аналогового решения одним из итерационных методов, например модифицированным методом последовательной верхней релаксации.Для точечной реализации метода для узловых точек с четной суммой (р+ц), где р - номер строки, о - номер столбца матрицы узловых процессоров, можно записать хк+1 = хк + (Ь - (бхк + И 1где х+1 - новое приближение решения; контроллера шины, переключателя локальных общих шин и коммутатора памятиЯ 1, Я - итерационные параметры ре- (1+1)-го узла вычисления системы конечно-разностных уравнений, третья группаФ о р м у л а и з о б р е т е н и я 5 информационных выходов-входов переклюУстройство для решения нелинейных чателя локальных общих шин каждого некраевых задач, содержащее регистр конфи- четного и четного 1-го узлов вычисленияраций, коммутатор и группу узлов вычис- системы конечно-разностных уравненийл ения системы конечно-разностных подключена к третьеи группе информациоуравнений, каждый иэ которых выполнен в 10 ных входов-выходов коммутаторами памятивиде дешифратора команд, блока коммута- соответственно (1-1)-го и (+1)-го узлов выции, шифратора, регистра данных и матри- числения системы конечно-разностныхцы блоков решения уравнений системы, уравнений, вторые группы информационпервые информационные входы и выходы ных входов вторых контроллеров шиныкоторых соединены с соответствующими 15 узлов вычисления системы конечно-раэновыходами и информационными входами стныхуравнений соединены между собой икоммутатора, подключенногоуправляющим с соответствующей группой информацивходом к выходу регистра конфигураций, онных входов регистра конфигураций ивторые информационные входы и выходы являются информационными входами устблоков матрицы 1-го ( = 2 - и - 1) узла вычис ройства, вторые группы входов-выходовления системы конечно-разностных уравне- третьих контроллеров шины узлов вычислений соединены соответственно с вторыми ния системы конечно-разностных уравневыходами и информационными входами со- ний соединены между собой, третья группаответствующих блоков матриц (1 - 1)-го и информационных выходов-входов третье(+1)-го узлов вычисления системы конечно го контроллера шины в каждом узле выразностных уравнений группы, третьи ин- числения системы конечно-разностныхформационные входы блоков матрицы уравнений подключена к четвертой групперешения уравнений системы подключены в информационных входов-выходов коммукаждом узле вычисления системы конечно- татора памяти, подключенного пятой групразностных уравнений через шифратор к 30 пой выходов и информационных входоввыходу регистра данных, а их выходы соеди- соответственно к адресным входам и выхонены с соответствующими информационны- дам блока памяти, выход аналого-цифровоми,входами блока коммутации, о т л и ч а ю- го преобразователя в каждом 1-ом узлещ е е с я тем, что, с целью повышения решения системы конечно-раэностныхбыстродействия, в каждый узел вычисления 35 уравнений соединен с соответствующим инсистемы конечно-разностных уравнений формационным входом первого контроллевве дены блок памяти, коммутатор памяти, ра шины, первыи, второи, третий ипервый, второй и третий контроллеры ши- четвертый информационные выход оны, дешифратор адреса, распределитель рого подключены к информационным вхоимпульсов, переключатель локальных об дам соответственно регистра данных,щих шин и аналого-цифровой преобраэова- дешифратора адреса, дешифратора командтель, причем в каждом -ом узле вычисления и к входу синхронизации распределителясистемы конечно-разностных уравнений импульсов, установочный вход которого сопервые группы информационных входов- единен с выходом дешифратора команд, авыходов и управляющие входы первого, вто первый, второй, третий и четвертый выходырого и третьего контроллеров шины подключены соответственно к управляюсоединены через первую локальную общую щим входам матрицы блоков решения уравшину, подключенную к управляющему вхо- нений системы, куправляющему входу узладу узла и к информациц онным выходам-вхо- коммутации, к входу запуска аналого-цифт обадам, с управляющим и ия им и информационными 50 рового преобразователя и к входу стра авходами-выходами коммутатора памяти и записи регистра данных, выходдешифратои реключателя общих шин -го узла, с вто- ра адреса подключен к адресным входамерыми группами ин, ормацф ационных входов- матрицы блоков решения уравнений сис евыходов и управляющими входами первого мы и блока коммутации,