C++矩阵库Armadillo出现warning solve() system is singular错误的解决

  本文介绍使用C++ 语言的矩阵库Armadillo时,出现报错system is singular; attempting approx solution的解决方法。

  在之前的文章中,我们介绍过Armadillo矩阵库在Visual Studio软件C++环境中的配置方法(https://blog.csdn.net/zhebushibiaoshifu/article/details/127123511),并且也介绍过Armadillo与OpenCV矩阵数据mat、vec与Mat的相互转换(https://blog.csdn.net/zhebushibiaoshifu/article/details/128274327)。而在使用Armadillo模块加以矩阵计算时,出现了1个问题——我这里会经常出现warning: solve(): system is singular; attempting approx solution字样的报错。而且因为我代码中调用Armadillo模块加以矩阵计算的次数非常多,导致这个报错会出现非常多次,如下图所示。

  虽然这个错误倒不会引发程序中断,程序还是会继续运行,相当于这个warning只是1个输出提示信息而已;但是,如果大量出现这种报错提示信息,会使得原本程序运行的速度被拖慢。如下图所示,我这里需要循环4710万次,如果一直出现这个报错的话就会使得程序运行很慢。

  因此,需要对这个报错加以解决。

  之所以出现这个错误,是因为Armadillo模块在执行矩阵计算时,遇到了秩缺陷系统,导致某些数值求解方法可能无法找到精确的解,进而导致其需要尝试找到近似解——此时,就会触发这个报错提示。其中,当线性方程组的系数矩阵不是满秩矩阵时,其被称为秩缺陷系统

  解决的思路也很简单,一个是直接隐藏这个级别的错误提示——相当于还是有这个问题,但是不让他再输出这个错误提示了;还有一个就是针对这个出错的原因加以处理。

  首先先看第一个思路。在Armadillo模块的官方网站中,提到了取消这个错误提示的方法,如下图所示。

  我们只需要在代码最开头,在#include <armadillo>这句代码前,添加如下的代码,从而将上述错误提示信息隐藏掉,在之后每次运行程序时屏蔽这个报错。

#define ARMA_WARN_LEVEL 1

  添加上述代码,如下图所示。

  此时,再运行我这里的程序,就会发现原本的错误提示信息消失了;如下图所示,不再出现原本的错误提示信息,而是直接输出了代码的进度提示。

  第二种思路,就是从报错的源头代码入手,避免出现这个错误。

  在本文的代码中,是在调用Armadillo模块的solve()函数时出现这个错误的;如下图所示。

  对于此,我们可以在solve()函数中添加solve_opts::no_approx这个参数,使得solve()函数不要寻找近似解;如下图所示。

  当然,这种方法就会带来1个新的问题——既然我们通过这种方式,设定了solve()函数不能找近似解;那么与此同时,如果精确解也无法找到的话,就会直接报错。如下图所示,本文这个场景下,solve()函数就无法找到精确解;那么如果我用上面的方式,不让他找近似解,就会直接报错。

  从上图的报错中点进来,可以发现报错的代码就是solve()函数这里;如下图所示。

  所以,对于本文这种场景,我们就只能用第1种思路,通过隐藏报错信息来解决这一问题了。

  至此,大功告成。

欢迎关注:疯狂学习GIS

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/614758.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【Java-TesseractOCR】通过Java实现OCR

通过Java实现OCR 一、TesseractOCR二、引入pom训练集下载地址三、引入训练集三、使用 一、TesseractOCR 本文使用的是TesseractOCR进行识别 二、引入pom <dependency><groupId>net.sourceforge.tess4j</groupId><artifactId>tess4j</artifactId&…

Android 加密之 打包为arr 项目依赖或者为jar

Android 加密之 打包为arr 项目依赖或者为jar 1. 修改build.gradle plugins {//id com.android.application// 1. 修改为libraryid com.android.library }android {namespace com.dzq.iccid2compileSdk 33defaultConfig {//applicationId "com.dzq.iccid2"// 2. 注…

如何让MacOS「终端」走代理

在 MacOS 操作系统中&#xff0c;默认情况下&#xff0c;终端命令行不会通过代理进行网络连接。这导致在应用软件研发过程中&#xff0c;许多需要通过命令行下载安装的软件或依赖包无法成功安装。经常出现Failed to connect to xxx port 443 after 75329 ms: Couldnt connect t…

mysql搭建主从

mysql搭建主从: 1:拉取mysql镜像 docker pull mysql2:创建主从对应目录 3:建立一个简易的mysql docker run -it --name mytest -e MYSQL_ROOT_PASSWORD123 -d mysql4:进入这个简易的mysql;从中获取my.cnf文件 docker exec -it mytest bash5:从容器中将my.cnf拷贝到 /3306/c…

商城系统个性化功能——可视化编辑

商城系统的普及化&#xff0c;让很多中小企业和商家也开始接触商城系统管理&#xff0c;之前在使用第三方平台时&#xff0c;大多数商城系统的样式都是固定的&#xff0c;商城页面也不能按照自己的想法去调整。 现在&#xff0c;随着商城系统越来越普遍&#xff0c;商城系统功…

unity在linux环境下videoplayer 无法播放问题解决路径

1、问题 一个项目需要在linux下播放视频&#xff0c;并且视频在机器上&#xff0c;也就是要使用应用外的视频文件进行播放。 视频的格式当前提供的事avi格式&#xff0c;并且使用videoplayer 在windows下播放正常。 但是发出包之后再Ubuntu环境怎么都无法播放。 2、测试环境…

这家动画公司,女神表情灵动秒杀90%的国漫女角色!

当3D国漫市场逐渐加入“内卷”的行列&#xff0c;从大的底层创作引擎UE的运用迭代&#xff0c;到细节的人物动捕、面捕技术的实际结合&#xff0c;在这场内卷的百舸争流中&#xff0c;涌现出一家家风格各异的头部国漫制作公司&#xff1a;有整体偏写实风格的原力动画&#xff0…

SpringBoot 微服务token 传递实现

1、前言 随着微服务的流行&#xff0c;微服务之间的安全以及业务需要&#xff0c;都需要穿递token &#xff0c;而token的传递一般通过header 头来传递。从架构的角度来讲 &#xff0c;一般的企业应用都由nginx、业务网关和各个微服务组成。这个nginx 传递header 我就不讲述。下…

Unity笔记之Android打包、减小包体之类的问题

打包问题 问题1&#xff1a; 一般大部分问题就是JDK、SDK、NDK之类的问题。现在是其他的问题&#xff0c;之前遇到过&#xff0c;好久没玩android了都忘了。 这试了半天&#xff0c;结果是需要有密钥库。那就给他创建一个填一下就行了 &#xff08;在网上看了半天&#xff…

React + three.js 实现人脸动捕与3D模型表情同步

系列文章目录 React 使用 three.js 加载 gltf 3D模型 | three.js 入门React three.js 3D模型骨骼绑定React three.js 3D模型面部表情控制React three.js 实现人脸动捕与3D模型表情同步 示例项目(github)&#xff1a;https://github.com/couchette/simple-react-three-facia…

什么是One-Class SVM

1. 简介 单类支持向量机&#xff0c;简称One-Class SVM(One-Class Support Vector Machine)&#xff0c;是一种用于异常检测的监督学习算法。其主要目标是找出数据集中的异常或罕见样本&#xff0c;而不需要大量的正常样本用于训练。这使其在处理高维数据和非常稀疏的异常检测问…

好文阅读-数据库-CREATE TABLE AS

添加链接描述 收获如下&#xff1a; 1 DROP DELETE TRUNCATE对比 2 CREATE TABLE AS 这种方式创建表不会复制表的索引&#xff0c;主键&#xff0c;外键约束&#xff0c;包括自增ID。因此应该使用 CREATE TABLE LIKE。 3 数据的处理可以通过建立临时表的方式&#xff0c;之后还…