15. Пакет LinearAlgebra. Элементарные операции с матрицами и векторами. Назначение и загрузка пакета LinearAlgebra
В новых реализациях систем Maple была сделана ставка на использование давно апробированных быстрых алгоритмов линейной алгебры, предложенных создателями Number Algorithm Group (NAG). Эти алгоритмы издавна применяются на больших ЭВМ и суперкомпьютерах, обеспечивая ускорение численных матричных операций от нескольких раз до нескольких десятков раз. Их применение обеспечивает эффективное использование систем символьной математики в решении задач, сводящихся к задачам линейной алгебры. В числе таких задач многочисленные задачи теоретической электротехники, механики многих объектов, моделирования электронных устройств и т.д.
В Maple 9.5/10 использование алгоритмов NAG реализуется пакетом LinearAlgebra. Для его загрузки используются следующие команды (файл NAG):
> restart; with(LinearAlgebra):
> infolevel[LinearAlgebra]:=1;
infolevelLinearAlgebra:= 1
Многие функции этого пакета (их большой список опущен) повторяет по назначению функции более старого пакета linalg, описанного выше. Поэтому мы не будем останавливаться на их повторном описании. Главное то, что эти функции задействуют возможности быстрых алгоритмов NAG и, в отличие от функций пакета linalg, ориентированы на численные расчеты в том формате обработки вещественных чисел, который характерен для применяемой компьютерной платформы. Знающий матричные методы читатель легко поймет назначение функций пакета LinearAlgebra по их составным названиям. Например, DeleteColumn означает удаление столбца матрицы, ToeplitzMatrix означает создание матрицы Теплица, ZeroMatrix — создание матрицы с нулевыми элементами и т.д. Все имена функций этого пакета начинаются с заглавной буквы.
Примеры матричных операций с применением пакета LinearAlgebra
Применение алгоритмов NAG особенно эффективно в том случае, когда используется встроенная в современные микропроцессоры арифметика чисел с плавающей запятой. С помощью специального флага такую арифметику можно отключать или включать:
> UseHardwareFloats := false; # use software floats
Use Hardware Floats := false
> UseHardwareFloats := true; # default behaviour
UseHardwareFloats := true
Матрицы в новом пакете линейной алгебры могут задаваться в угловых скобках, как показано ниже:
> М1:=<<1|2>,<4|5>>; М2:=<<1|2.>, <4|5>>;
После этого можно выполнять с ними типовые матричные операции. Например, можно инвертировать (обращать) матрицы:
> М1^(-1); М2^(-1);
MatrixInverse: "calling external function"
MatrixInverse: "NAG" hw_f07adf
MatrixInverse: "NAG" hw_f07ajf
Обратите внимание, что Maple теперь выдает информационные сообщения о новых условиях реализации операции инвертирования матриц с вещественными элементами и, в частности, об использовании алгоритмов NAG и арифметики, встроенной в сопроцессор.
Следующий пример иллюстрирует создание двух случайных матриц M1 и М2 и затем их умножение:
> M1:=RandomMatrix(2,3); М2:=RandomMatrix(3,3);
Multiply(M1,M2,'inplace'); M1;
Параметр inplace в функции умножения обеспечивает помещение результата умножения матриц на место исходной матрицы М1 — излюбленный прием создателей быстрых матричных алгоритмов NAG. Поскольку матрицы M1 и М2 заданы как случайные, то при повторении этого примера результаты, естественно, будут иными, чем приведенные.
Другой пример иллюстрирует проведение хорошо известной операции LU-разложения над матрицей М, созданной функцией Matrix:
> M:=Matrix([[14,-8,1],[-11,-4,18],[3,12,19]], datatype=float);
LUDecomposition(М,output=['NAG'],inplace);
ipiv:=%[1];
M;
LUDecomposition: "calling external function"
LUDecomposition: "NAG" hw_f07adf
