个人中心
文章发布
退出
作者热门文章
- iOS/Objective-C 元类和类别
- objective-c - -1001 错误，当 NSURLSession 通过 httpproxy 和/etc/hosts
- java - 使用网络类获取 url 地址
- ios - 推送通知中不播放声音
26
4
最近我一直致力于使用 Fortran 解决一些数学问题。我想到一个线性矩阵方程:
AX-XA=B
其中 A 和 B 是已知的 n*n 矩阵，X 是需要求解的矩阵。我知道这看起来像一个典型的李雅普诺夫方程。但是，为了空间友好地求解这个方程，可以使用迭代的方式来得到当克罗内克积时X的数值结果:
A*I-I*A(T)
是非奇异的(T 表示转置)。
不幸的是，这不是我的情况。我的问题中的克罗内克积是奇异的。因此，我不能使用空间友好的迭代方法来解决问题，而是使用非常耗费空间的方法来显式生成和存储 Kronecker 产品!!!然后我使用矩阵算法的 Moore-Penrose pseudoinverse 来生成这个巨大矩阵的伪逆，并使用库矩阵向量乘法例程来求解方程。
对于小矩阵，这没问题。但是当矩阵的大小增加时(例如当 n 增加到几百时)我的计算机内存已经完全用完了。任何人都可以帮助解决这个问题，以便我可以使用空间友好的算法来解决这样的问题吗？非常感谢。
小号
最佳答案
这不是严格意义上的答案，但它可能会帮助您取得一些进步。如果我正确理解你的问题，你想解决一个西尔维斯特方程，即这里 http://en.wikipedia.org/wiki/Sylvester_equation
你有没有试过LAPACK，也许:http://www.netlib.org/lapack/lawn41/node112.html -> 搜索 TGSY2
或者可能是这样的:http://manpages.ubuntu.com/manpages/precise/man3/stgsyl.3lapack.html
甚至可能:http://www8.cs.umu.se/~isak/recsy/
祝你好运:)
关于algorithm - 如何求解线性矩阵方程 : AX-XA=B efficiently?，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/22987516/
26
4
0
我有两种结构，Header 和Session，它们都符合协议(protocol)TimelineItem。 我有一个 Array 由 TimelineItem 组成，如下所示: [Header1, S
这个问题在这里已经有了答案: Multiple assignment and evaluation order in Python (11 个答案) 关闭 6 年前。 我刚接触python所以想问你
我试图找到一种方法来在 R 中获取 A、A、A、A、B、B、B、B、B 的所有可能的唯一排列的列表。 组合最初被认为是获得解决方案的方法，因此组合的答案。 最佳答案 我认为这就是你所追求的。 @bil
我怎样才能将两个给定的向量混合成一个新的向量，它以交替的顺序保存它们的值。 (f [a a] [b b]) ; > [a b a b] 这是我想到的: (flatten (map vector [:a
这是我的第一个问题，我开始学习Python。之间有区别吗: a, b = b, a + b 和 a = b b = a + b 当您在下面的示例中编写它时，它会显示不同的结果。 def fib(n):
这个问题在这里已经有了答案: Why is there an injected class name? (1 个回答) 12 个月前关闭。 我不知道如何解释: namespace A { struct
我尝试了一些代码来交换 Java 中的两个整数，而不使用第三个变量，使用 XOR。 这是我尝试过的两个交换函数: package lang.numeric; public class SwapVars
假设类 B 扩展类 A，并且我想为 B 声明一个变量。什么更有效？为什么？ B b或 A b . 最佳答案 您混淆了两个不同的概念。 class B extends A { } 意味着B 是 A .
我不确定这个问题的标题是什么，这也可能是一个重复的问题。所以请相应地指导。 我是 python 编程的新手。我有这个简单的代码来生成斐波那契数列。 1: def fibo(n): 2: a =
我在谷歌上搜索了有关 dynamic_cast 的内容，我发现显式地将基类对象转换为派生类指针可能是不安全的。但是当我运行一些示例代码来检查它时，我没有收到任何错误。请在下面找到我的代码: class
这个问题在这里已经有了答案: What is this weird colon-member (" : ") syntax in the constructor? (14 个答案) 关闭 8 年前。
在不重现产生非整数值的表达式的情况下实现以下目标的惯用方法是什么(在我的真实情况下，该值是在我不想重现的冗长查询之后计算为百分比的): SELECT * FROM SomeTable WHERE 1/
在析构中，这两个代码的结果确实不同。我不确定为什么。 提示说 const [b,a] = [a,b] 将导致 a,b 的值为 undefined (从左到右的简单分配规则)。我不明白为什么会这样。 l
C++ Templates - The Complete Guide, 2nd Edition介绍max模板: template T max (T a, T b) { // if b < a th
我最近开始学习代码(Java)，并根据第 15.17.3 节在 Oracle 网站上查找了模运算符。以下链接: http://docs.oracle.com/javase/specs/jls/se8/
无法理解以下行为。 d1 := &data{1}; 的区别d1 和 d2 := 数据{1}； &d1。两者都是指针，对吧？但他们的行为不同。这里发生了什么 package main import "f
这个问题在这里已经有了答案: How to make loop infinite with "x = y && x != y"? (4 个回答) How can i define variables
在我的程序中，当我调试我的代码时，它似乎在我生成的代码中的某处 X1=['[a,a,a]','[b,b,b]'] 还有我生成的其他地方 X2=[[a,a,a],[b,b,b]] 当我想添加这两个列表然
我试图使用递归将两个整数相乘，并意外编写了这段代码: //the original version int multiply(int a, int b) { if ( !b ) retu
我有一个列表中数字之间所有可能的操作组合: list = ['2','7','8'] 7+8*2 8+7*2 2*8+7 2+8*7 2-8*7 8-2/7 etc 我想知道是否可以说像 ('7*2+
我是一名优秀的程序员,十分优秀！