Название вопроса. Этот метод включает использование «матрицы кофакторов» или «матрицы сопряжения» и дает явные формулы для компонентов обратной квадратной матрицы. Нелегко сделать вручную для матрицы больше, чем, скажем, 3 × 33×33\times 3 . Для матрицы n × nn×Nn\times n требуется вычисление...
В программном проекте, над которым я работаю, некоторые вычисления намного проще для плотных матриц низкого ранга. В некоторых проблемных случаях используются плотные матрицы низкого ранга, но они даны мне полностью, а не как факторы, поэтому мне придется проверять ранг и фактор матрицы, если я...
Вопрос: Какие методы доступны для точного и эффективного расчета разреженной структуры матрицы конечных элементов? Информация: я работаю над решателем уравнения давления Пуассона, использую метод Галеркина с квадратичной базой Лагранжа, написанный на C, и использую PETSc для хранения разреженных...
Мне интересно, является ли алгоритм Томаса самым быстрым (доказуемо?) Решением симметричной диагонально доминирующей разреженной трехдиагональной системы с точки зрения алгоритмической сложности (не ища пакетов реализации, таких как LAPACK и т. Д.). Я знаю, что и алгоритм Томаса, и многосетка имеют...
Числа с плавающей запятой одинарной точности занимают половину памяти, и на современных машинах (кажется, даже на графических процессорах) операции могут выполняться с ними почти вдвое быстрее, чем с двойной точностью. Многие коды FDTD, которые я обнаружил, используют исключительно арифметику...
Я работаю над библиотекой матриц только для заголовков, чтобы обеспечить некоторую разумную степень возможностей линейной алгебры в максимально простом пакете, и я пытаюсь рассмотреть, каково текущее состояние техники: вычисление SVD сложная матрица. Я делаю двухфазную декомпозицию,...
Я ищу библиотеку, которая выполняет матричные операции над большими разреженными матрицами без ущерба для числовой стабильности. Матрицы будут 1000+ на 1000+, а значения матрицы будут между 0 и 1000. Я буду выполнять алгоритм исчисления индекса, поэтому буду генерировать (разреженные) векторы строк...
Учитывая две матрицы и , я хотел бы найти векторы и , такие, что В матричной форме я пытаюсь минимизировать норму Фробениуса для A - \ mbox {diag} (x) \ cdot B \ cdot \ mbox {diag} (y) = A - B \ circ (xy ^ \ top) .AAABBBxxxyyymin∑ij(Aij−xiyjBij)2.min∑ij(Aij−xiyjBij)2. \min \sum_{ij} (A_{ij} - x_i...
Мне нужно решить ту же самую разреженную линейную систему (от 300x300 до 1000x1000) со многими правыми сторонами (от 300 до 1000). В дополнение к этой первой проблеме, я также хотел бы решить различные системы, но с одинаковыми ненулевыми элементами (только с разными значениями), то есть с...
Я использую MATLAB для решения проблемы, которая включает в себя решение на каждом временном шаге, где b изменяется со временем. Прямо сейчас я делаю это, используя MATLAB :A x = bAx=b\mathbf{A} \mathbf{x}=\mathbf{b}бb\mathbf{b}mldivide x = A\b У меня есть возможность делать столько предварительных...
Учитывая систему где A ∈ R n × n , я прочитал, что, если итерация Якоби используется в качестве решателя, метод не будет сходиться, если b имеет ненулевую компоненту в нуль-пространстве A , Итак, как можно формально утверждать, что при условии, что b имеет ненулевой компонент, охватывающий нулевое...
Мне дают матрицу Q которая является симметричной, обратимой, положительно определенной и плотной. Мне нужно проверить, если det, где J является матрицей всех единиц.12 × 1212×1212 \times 12QQQdet ( Q ) = det ( 12 I- Q - J)( 1 )йе(Q)знак равнойе(12я-Q-J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ В...
Как MATLAB, например, вычисляет SVD данной матрицы? Я предполагаю, что ответ, вероятно, включает в себя вычисление собственных векторов и собственных значений A*A'. Если это так, я также хотел бы знать, как он вычисляет...
Пожалуйста, перечислите пакет Python (petsc4py и т. Д.) И разреженные прямые решатели, которые он поддерживает. Один (Community-Wiki) ответ на пакет,
Я решаю дифференциальные уравнения, которые требуют инвертировать плотные квадратные матрицы. Эта инверсия матрицы занимает большую часть моего времени вычислений, поэтому мне было интересно, использую ли я самый быстрый из доступных алгоритмов. Мой текущий выбор - numpy.linalg.inv . Из моих чисел...
Предположим следующую матрицу дается [ 0,500 - 0,333 - 0,167 - 0,500 0,667 - 0,167 - 0,500 - 0,333 0,833 ] с транспонированной T . Продукт A T A = G дает [ 0,750 - 0,334 - 0,417 - 0,334 0,667 - 0,333 - 0,417 - 0,333 0,750 ]...
Следующее матричное уравнение в Σ - для заданных матриц B и C - фигурирует в моей работе как характеристика ковариационной матрицы. Я узнал, что это уравнение известно, в частности, в теории непрерывного управления временем, как уравнение Ляпунова , и что существуют различные хорошо известные...
Есть ли надежда на эффективное решение следующей линейной системы итерационным методом? A ∈ Rn × n, x ∈ RN, b ∈ RNс n > 106A∈Rn×n,x∈Rn,b∈Rn, with n>106A \in \mathbb{R}^{n \times n}, x \in \mathbb{R}^n, b \in \mathbb{R}^n \text{, with } n > 10^6 A x = bAx=bAx=b с участием A = ( Δ - K)A=(Δ−K)...
Я новичок в пакете PETSc. У меня матрица A размером ~ 4000x4000 в формате матричного рынка, и я хочу, чтобы PETSc решил эту проблему с помощью нескольких процессоров. Я знаю, как решить систему на одном процессоре, но я не знаю, как распределить матрицу и векторы между различными процессорами. Есть...
Я пытаюсь диагонализировать некоторые плотные, плохо обусловленные матрицы. В машинной точности результаты являются неточными (возвращая отрицательные собственные значения, собственные векторы не имеют ожидаемой симметрии). Я переключился на функцию Eigensystem [] Mathematica, чтобы использовать...