?????????????? ?????? ? ????????? ?????????? ??????? ?????

алгоритма имеем матрицу сопротивлений ветвей zv, массив источников
тока ветвей jv, массив источников напряжений ветвей ev.
В пятом блоке формируются топологические ортогональные матрицы:
контурная и узловая. Для этой цели используется подпрограмма
derevo. Вычисление контурной матрицы выполняется инвертированием
узловой матрицы. Для инвертирования матрицы используются подпрограммы
библиотеки SLATEC : sgefa и sgedi .
В шестом блоке на основе метода контурных токов вычисляются массивы
напряжений ветвей uv и токов ветвей iv. Вычисления производятся
по тем же формулам, что и в алгоритме на SciLAB раздела 2.3.7. Для инвертирования
матрицы контурных сопротивлений использована подпрограмма
inverse4, [16],текст которой приводится за основным алгоритмом.
Для перемножения матриц используется встроенная Фортрановская
функция matmul.
В седьмом блоке выполняется вывод результатов расчета — токов и
напряжений ветвей в конец входного файла.
Для оценки времени счета используется обращение к встроенной
Фортран-подпрограмме CPU_TIME При первом обращении в переменную
rt1 записывается процессорное время прошедшее с начала загрузки
программы. Следующее обращение к этой подпрограмме происходит
перед формированием топологических матриц, процессорное время записывается
в rt2 и последнее обращение происходит после расчетов,
перед выводом результатов в файл, процессорное время записывается в
rt3. Таким образом, разница rt2 и rt1 соответствует подготовке информации
для расчета, разница rt3 и rt2 соответствует процессорному времени,
затрачиваемому непосредственно на расчет схемы. В DOS-окно
время на подготовку и время на расчет выводятся соответственно под
именами time1 и time2.
87