- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
从 Eigen 3.2.7 更新到 3.3.4 后,我遇到了 JacobiSVD.solve 的问题,它返回了一个非常错误的结果。BDCSVD 产生相同的结果。
可以使用以下代码重现该问题:
#include <Eigen/Eigen>
#include <Eigen/SVD>
int main(int argc, const char * argv[]) {
// M is calculated beforehand and does not change with different Eigen versions
Eigen::Matrix<double, 12, 12> M = Eigen::Matrix<double, 12, 12>::Zero();
M(0,0) = 27; M(0,2) = 4.3625039999999995; M(0,3) = -9; M(0,5) = -2.1812519999999997; M(0,6) = -9;
M(1,1) = 27; M(1,2) = 3.2718720000000001; M(1,4) = -9; M(1,7) = -9; M(1,8) = -1.6359360000000001;
M(2,0) = 4.3625039999999995; M(2,1) = 3.2718720000000001; M(2,2) = 2.4780489612000003;
M(2,3) = -2.1812519999999997; M(2,5) = -0.82601632039999995; M(2,7) = -1.6359360000000001;
M(2,8) = -0.82601632039999995; M(3,0) = -9; M(3,2) = -2.1812519999999997; M(3,3) = 9;
M(4,1) = -9; M(4,4) = 9; M(5,0) = -2.1812519999999997; M(5,2) = -0.82601632039999995;
M(5,5) = 0.82601632039999995; M(6,0) = -9; M(6,6) = 9; M(7,1) = -9; M(7,2) = -1.6359360000000001;
M(7,7) = 9; M(8,1) = -1.6359360000000001; M(8,2) = -0.82601632039999995; M(8,8) = 0.82601632039999995;
Eigen::JacobiSVD<Eigen::MatrixXd> svd(M, Eigen::ComputeFullU);
Eigen::Matrix<double, 12, 12> ut = svd.matrixU();
Eigen::Matrix<double, 6, 10> l_6x10;
Eigen::Matrix<double, 4, 6> dv0;
Eigen::Matrix<double, 4, 6> dv1;
Eigen::Matrix<double, 4, 6> dv2;
for(int i = 0; i < 4; i++) {
int a = 0, b = 1;
for(int j = 0; j < 6; j++) {
dv0(i, j) = ut(3 * a + 0, 11 - i) - ut(3 * b + 0, 11 - i);
dv1(i, j) = ut(3 * a + 1, 11 - i) - ut(3 * b + 1, 11 - i);
dv2(i, j) = ut(3 * a + 2, 11 - i) - ut(3 * b + 2, 11 - i);
b++;
if (b > 3) {
a++;
b = a + 1;
}
}
}
for(int i = 0; i < 6; i++) {
l_6x10(i,0) = (dv0(0, i) * dv0(0, i) + dv1(0, i) * dv1(0, i) + dv2(0, i) * dv2(0, i));
l_6x10(i,1) = 2.0f * (dv0(0, i) * dv0(1, i) + dv1(0, i) * dv1(1, i) + dv2(0, i) * dv2(1, i));
l_6x10(i,2) = (dv0(1, i) * dv0(1, i) + dv1(1, i) * dv1(1, i) + dv2(1, i) * dv2(1, i));
l_6x10(i,3) = 2.0f * (dv0(0, i) * dv0(2, i) + dv1(0, i) * dv1(2, i) + dv2(0, i) * dv2(2, i));
l_6x10(i,4) = 2.0f * (dv0(1, i) * dv0(2, i) + dv1(1, i) * dv1(2, i) + dv2(1, i) * dv2(2, i));
l_6x10(i,5) = (dv0(2, i) * dv0(2, i) + dv1(2, i) * dv1(2, i) + dv2(2, i) * dv2(2, i));
l_6x10(i,6) = 2.0f * (dv0(0, i) * dv0(3, i) + dv1(0, i) * dv1(3, i) + dv2(0, i) * dv2(3, i));
l_6x10(i,7) = 2.0f * (dv0(1, i) * dv0(3, i) + dv1(1, i) * dv1(3, i) + dv2(1, i) * dv2(3, i));
l_6x10(i,8) = 2.0f * (dv0(2, i) * dv0(3, i) + dv1(2, i) * dv1(3, i) + dv2(2, i) * dv2(3, i));
l_6x10(i,9) = (dv0(3, i) * dv0(3, i) + dv1(3, i) * dv1(3, i) + dv2(3, i) * dv2(3, i));
}
Eigen::Matrix<double, 6, 4> L_6x4;
for(int i = 0; i < 6; i++) {
L_6x4(i, 0) = l_6x10(i, 0);
L_6x4(i, 1) = l_6x10(i, 1);
L_6x4(i, 2) = l_6x10(i, 3);
L_6x4(i, 3) = l_6x10(i, 6);
}
// L_6x4 has the following values on 3.2.7 (everything else is 0):
//
// L_6x4(0,2) = 1;
// L_6x4(1,0) = 1;
// L_6x4(1,1) = 1;
// L_6x4(5,0) = -1.137432760287006;
// L_6x4(5,2) = -1.1374327602870071;
// L_6x4(5,3) = -1.1374327602870049;
//
//
// on 3.3.4 it has the following slightly different values:
//
// L_6x4(0,2) = 1;
// L_6x4(1,0) = 1;
// L_6x4(1,1) = 1;
// L_6x4(5,0) = -1.1374327602869998;
// L_6x4(5,2) = -1.1374327602870271;
// L_6x4(5,3) = -1.1374327602869889;
// Rho is calculated beforehand and does not change with different Eigen versions
Eigen::Matrix<double, 6, 1> Rho;
Rho << 0.25, 0.140625, 0, 0.390625, 0.25, 0.140625;
Eigen::JacobiSVD<Eigen::MatrixXd> l_6x4(L_6x4, Eigen::ComputeFullU | Eigen::ComputeFullV);
Eigen::Vector4d B4 = l_6x4.solve(Rho);
// The slight difference in L_6x4 apparently causes an insane difference in the result of solve.
//
// Eigen 3.2.7: B4={0.056766494562302421, 0, 0, -0.064568070601816963}
// Eigen 3.3.4: B4={-4773392957911.6992, 0, 0, -4196637484493.9165}
return 0;
}
使用 Eigen 3.3.4 计算的 B4 值让我认为这可能是一些奇怪的内存对齐问题、SDK 中的错误或希望是该程序中的错误。
我尝试使用 -DEIGEN_DONT_ALIGN_STATICALLY 或 -DEIGEN_DONT_VECTORIZE 和 -DEIGEN_DISABLE_UNALIGNED_ARRAY_ASSERT 对其进行编译。我还在这些矩阵上尝试了 Eigen::DontAlign。这些都对结果没有任何影响。
使用 ndk-r14b clang c++_static 在 Android(4.4 和 5.1)上测试和 Mac(macOs Sierra 10.12.6):“clang++ -std=c++14 -g -O0 -I/EIGEN_PATH main.cpp -o main”
最佳答案
这两个答案实际上非常相似并且产生相似的残差。这是因为您的矩阵有两个零奇异值,第三个奇异值非常接近机器精度。因此,根据舍入误差,它被认为是 0,在这种情况下,您将在解决方案的剩余 3D 子空间内获得最小范数解决方案:
{0.056766494562302421, 0, 0, -0.064568070601816963}
和 0.91 的残差。否则认为是有意义的,然后得到更小的二维子空间内的最小范数解(对应两个零奇异值):
{-4773392957911.6992, 0, 0, -4196637484493.9165}
残差较小,为 0.89。所以第二个解决方案在某种意义上更准确,但如果您更喜欢范数较小但残差较高的解决方案,则可以调整阈值,使第三个奇异值被视为零。这是通过 setThreshold 完成的:
Eigen::JacobiSVD<Eigen::MatrixXd> l_6x4;
l_6x4.setThreshold(1e-14);
l_6x4.compute(L_6x4, Eigen::ComputeFullU | Eigen::ComputeFullV);
Eigen::Vector4d B4 = l_6x4.solve(Rho);
关于c++ - Eigen 3.3 SVD.solve 返回错误值,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/47460151/
我已经使用 vue-cli 两个星期了,直到今天一切正常。我在本地建立这个项目。 https://drive.google.com/open?id=0BwGw1zyyKjW7S3RYWXRaX24tQ
您好,我正在尝试使用 python 库 pytesseract 从图像中提取文本。请找到代码: from PIL import Image from pytesseract import image_
我的错误 /usr/bin/ld: errno: TLS definition in /lib/libc.so.6 section .tbss mismatches non-TLS reference
我已经训练了一个模型,我正在尝试使用 predict函数但它返回以下错误。 Error in contrasts<-(*tmp*, value = contr.funs[1 + isOF[nn]])
根据Microsoft DataConnectors的信息我想通过 this ODBC driver 创建一个从 PowerBi 到 PostgreSQL 的连接器使用直接查询。我重用了 Micros
我已经为 SoundManagement 创建了一个包,其中有一个扩展 MediaPlayer 的类。我希望全局控制这个变量。这是我的代码: package soundmanagement; impo
我在Heroku上部署了一个应用程序。我正在使用免费服务。 我经常收到以下错误消息。 PG::Error: ERROR: out of memory 如果刷新浏览器,就可以了。但是随后,它又随机发生
我正在运行 LAMP 服务器,这个 .htaccess 给我一个 500 错误。其作用是过滤关键字并重定向到相应的域名。 Options +FollowSymLinks RewriteEngine
我有两个驱动器 A 和 B。使用 python 脚本,我在“A”驱动器中创建一些文件,并运行 powerscript,该脚本以 1 秒的间隔将驱动器 A 中的所有文件复制到驱动器 B。 我在 powe
下面的函数一直返回这个错误信息。我认为可能是 double_precision 字段类型导致了这种情况,我尝试使用 CAST,但要么不是这样,要么我没有做对...帮助? 这是错误: ERROR: i
这个问题已经有答案了: Syntax error due to using a reserved word as a table or column name in MySQL (1 个回答) 已关闭
我的数据库有这个小问题。 我创建了一个表“articoli”,其中包含商品的品牌、型号和价格。 每篇文章都由一个 id (ID_ARTICOLO)` 定义,它是一个自动递增字段。 好吧,现在当我尝试插
我是新来的。我目前正在 DeVry 在线学习中级 C++ 编程。我们正在使用 C++ Primer Plus 这本书,到目前为止我一直做得很好。我的老师最近向我们扔了一个曲线球。我目前的任务是这样的:
这个问题在这里已经有了答案: What is an undefined reference/unresolved external symbol error and how do I fix it?
我的网站中有一段代码有问题;此错误仅发生在 Internet Explorer 7 中。 我没有在这里发布我所有的 HTML/CSS 标记,而是发布了网站的一个版本 here . 如您所见,我在列中有
如果尝试在 USB 设备上构建 node.js 应用程序时在我的树莓派上使用 npm 时遇到一些问题。 package.json 看起来像这样: { "name" : "node-todo",
在 Python 中,您有 None单例,在某些情况下表现得很奇怪: >>> a = None >>> type(a) >>> isinstance(a,None) Traceback (most
这是我的 build.gradle (Module:app) 文件: apply plugin: 'com.android.application' android { compileSdkV
我是 android 的新手,我的项目刚才编译和运行正常,但在我尝试实现抽屉导航后,它给了我这个错误 FAILURE: Build failed with an exception. What wen
谁能解释一下?我想我正在做一些非常愚蠢的事情,并且急切地等待着启蒙。 我得到这个输出: phpversion() == 7.2.25-1+0~20191128.32+debian8~1.gbp108
我是一名优秀的程序员,十分优秀!