大规模开源线性代数求解器(Eigen,LAPACK,Ceres)+JSim数值解算器+PlotDigitizer

云深之无迹 2021-12-07 23:59



看见一个招聘的要求,说有这个的经验最好。就搜索了一下。


https://eigen.tuxfamily.org/index.php?title=Main_Page

一个C++的计算矩阵的库


#include #include 
using Eigen::MatrixXd;
int main(){ MatrixXd m(2,2); m(0,0) = 3; m(1,0) = 2.5; m(0,1) = -1; m(1,1) = m(1,0) + m(0,1); std::cout << m << std::endl;}


没有什么依赖的库,就是标准库


https://eigen.tuxfamily.org/dox/group__TutorialMatrixClass.html


这是他的长篇教程



http://www.netlib.org/lapack/

LAPACK 是用 Fortran 90 编写的,提供用于求解联立线性方程组、线性方程组的最小二乘解、特征值问题和奇异值问题的例程。还提供了相关的矩阵分解(LU、Cholesky、QR、SVD、Schur、广义 Schur),以及相关计算,例如 Schur 分解的重新排序和估计条件数。处理密集矩阵和带状矩阵,但不处理一般稀疏矩阵。在所有领域,都为单精度和双精度实数和复数矩阵提供了类似的功能。

LAPACK 项目的最初目标是使广泛使用的 EISPACK 和 LINPACK 库在共享内存向量和并行处理器上高效运行。在这些机器上,LINPACK 和 EISPACK 效率低下,因为它们的内存访问模式不考虑机器的多层内存层次结构,从而花费太多时间移动数据而不是进行有用的浮点运算。LAPACK 通过重新组织算法以在最内层循环中使用块矩阵运算(例如矩阵乘法)来解决此问题。这些块操作可以针对每个架构进行优化,以考虑内存层次结构,从而提供一种可移植的方式来在不同的现代机器上实现高效率。我们使用术语“便携式”而不是“便携式”,因为,编写 LAPACK 例程,以便通过调用基本线性代数子程序 (BLAS) 来执行尽可能多的计算。LAPACK 一开始就被设计为利用 Level 3 BLAS——一组 Fortran 子程序的规范,用于执行各种类型的矩阵乘法和具有多个右手边的三角系统的解决方案。由于 3 级 BLAS 操作的粗粒度,它们的使用提高了许多高性能计算机的效率,特别是如果制造商提供了特殊编码的实现。

BLAS 的高效机器特定实现可用于许多现代高性能计算机。有关已知供应商或 ISV 提供的 BLAS 的详细信息,请参阅 BLAS 常见问题解答。或者,用户可以下载 ATLAS 以自动为架构生成优化的 BLAS 库。BLAS 的 Fortran 77 参考实现可从 netlib 获得;但是,不鼓励使用它,因为它的性能不如经过专门调整的实现。


但是吧,太复杂了。


http://www.netlib.org/lapack/lug/node3.html


给了一本书的地址,可以看看。


http://ceres-solver.org/

Google的计算库,看起来不错,名字叫做谷神星








反正挺多的

https://github.com/ceres-solver/ceres-solver


开源地址在此。


其实对于我这种数值算法爱好者来说,这些库真是再好不过的学习资源了。


对了,看论文的时候看到一个有趣的东西。是一个Java的数值计算包


https://www.physiome.org/


文章就是这样写的


https://www.physiome.org/jsim/docs/MML_Matlab.html
// Example for MML to Matlab translationimport nsrunit; unit conversion on;
math mml2matlab { realDomain t sec; t.min=0; t.max=1; t.delta=0.01; real C(t) mM; real k =1 sec^(-1); when(t=t.min) C=1; C:t = -k*C; }


我们写的仿真程序

options = odeset('RelTol',1e-12,'AbsTol',1e-9);[t x rInfo] = mytest( linspace(0,1,101),@ode23,options);clf;plot(t,x,'k')title('Exponential Decay')xlabel('Time, sec')ylabel('Concentration, mM')



http://plotdigitizer.sourceforge.net

Plot Digitizer 是一个 Java 程序,用于数字化功能数据的扫描图。在报告和参考文献中,数据通常以功能 XY 型散点图或折线图的形式呈现。为了使用这些数据,必须以某种方式将其数字化。该程序将允许您获取绘图的扫描图像(GIF、JPEG 或 PNG 格式),只需在每个数据点上单击鼠标即可快速将绘图中的值数字化。然后可以将这些数字保存到文本文件中,并在您需要的任何地方使用。Plot Digitizer 适用于线性和对数轴刻度。除了数字化数据图中的点外,该程序还可用于数字化其他类型的扫描数据(例如缩放绘图或正字照片)。



评论
  • 遇到部分串口工具不支持1500000波特率,这时候就需要进行修改,本文以触觉智能RK3562开发板修改系统波特率为115200为例,介绍瑞芯微方案主板Linux修改系统串口波特率教程。温馨提示:瑞芯微方案主板/开发板串口波特率只支持115200或1500000。修改Loader打印波特率查看对应芯片的MINIALL.ini确定要修改的bin文件#查看对应芯片的MINIALL.ini cat rkbin/RKBOOT/RK3562MINIALL.ini修改uart baudrate参数修改以下目
    Industio_触觉智能 2024-12-03 11:28 84浏览
  • TOF多区传感器: ND06   ND06是一款微型多区高集成度ToF测距传感器,其支持24个区域(6 x 4)同步测距,测距范围远达5m,具有测距范围广、精度高、测距稳定等特点。适用于投影仪的无感自动对焦和梯形校正、AIoT、手势识别、智能面板和智能灯具等多种场景。                 如果用ND06进行手势识别,只需要经过三个步骤: 第一步&
    esad0 2024-12-04 11:20 50浏览
  •         温度传感器的精度受哪些因素影响,要先看所用的温度传感器输出哪种信号,不同信号输出的温度传感器影响精度的因素也不同。        现在常用的温度传感器输出信号有以下几种:电阻信号、电流信号、电压信号、数字信号等。以输出电阻信号的温度传感器为例,还细分为正温度系数温度传感器和负温度系数温度传感器,常用的铂电阻PT100/1000温度传感器就是正温度系数,就是说随着温度的升高,输出的电阻值会增大。对于输出
    锦正茂科技 2024-12-03 11:50 106浏览
  • RDDI-DAP错误通常与调试接口相关,特别是在使用CMSIS-DAP协议进行嵌入式系统开发时。以下是一些可能的原因和解决方法: 1. 硬件连接问题:     检查调试器(如ST-Link)与目标板之间的连接是否牢固。     确保所有必要的引脚都已正确连接,没有松动或短路。 2. 电源问题:     确保目标板和调试器都有足够的电源供应。     检查电源电压是否符合目标板的规格要求。 3. 固件问题: &n
    丙丁先生 2024-12-01 17:37 100浏览
  • 最近几年,新能源汽车愈发受到消费者的青睐,其销量也是一路走高。据中汽协公布的数据显示,2024年10月,新能源汽车产销分别完成146.3万辆和143万辆,同比分别增长48%和49.6%。而结合各家新能源车企所公布的销量数据来看,比亚迪再度夺得了销冠宝座,其10月新能源汽车销量达到了502657辆,同比增长66.53%。众所周知,比亚迪是新能源汽车领域的重要参与者,其一举一动向来为外界所关注。日前,比亚迪汽车旗下品牌方程豹汽车推出了新车方程豹豹8,该款车型一上市就迅速吸引了消费者的目光,成为SUV
    刘旷 2024-12-02 09:32 119浏览
  • 概述 说明(三)探讨的是比较器一般带有滞回(Hysteresis)功能,为了解决输入信号转换速率不够的问题。前文还提到,即便使能滞回(Hysteresis)功能,还是无法解决SiPM读出测试系统需要解决的问题。本文在说明(三)的基础上,继续探讨为SiPM读出测试系统寻求合适的模拟脉冲检出方案。前四代SiPM使用的高速比较器指标缺陷 由于前端模拟信号属于典型的指数脉冲,所以下降沿转换速率(Slew Rate)过慢,导致比较器检出出现不必要的问题。尽管比较器可以使能滞回(Hysteresis)模块功
    coyoo 2024-12-03 12:20 108浏览
  • 戴上XR眼镜去“追龙”是种什么体验?2024年11月30日,由上海自然博物馆(上海科技馆分馆)与三湘印象联合出品、三湘印象旗下观印象艺术发展有限公司(下简称“观印象”)承制的《又见恐龙》XR嘉年华在上海自然博物馆重磅开幕。该体验项目将于12月1日正式对公众开放,持续至2025年3月30日。双向奔赴,恐龙IP撞上元宇宙不久前,上海市经济和信息化委员会等部门联合印发了《上海市超高清视听产业发展行动方案》,特别提到“支持博物馆、主题乐园等场所推动超高清视听技术应用,丰富线下文旅消费体验”。作为上海自然
    电子与消费 2024-11-30 22:03 98浏览
  • 作为优秀工程师的你,已身经百战、阅板无数!请先醒醒,新的项目来了,这是一个既要、又要、还要的产品需求,ARM核心板中一个处理器怎么能实现这么丰富的外围接口?踌躇之际,你偶阅此文。于是,“潘多拉”的魔盒打开了!没错,USB资源就是你打开新世界得钥匙,它能做哪些扩展呢?1.1  USB扩网口通用ARM处理器大多带两路网口,如果项目中有多路网路接口的需求,一般会选择在主板外部加交换机/路由器。当然,出于成本考虑,也可以将Switch芯片集成到ARM核心板或底板上,如KSZ9897、
    万象奥科 2024-12-03 10:24 68浏览
  • 光伏逆变器是一种高效的能量转换设备,它能够将光伏太阳能板(PV)产生的不稳定的直流电压转换成与市电频率同步的交流电。这种转换后的电能不仅可以回馈至商用输电网络,还能供独立电网系统使用。光伏逆变器在商业光伏储能电站和家庭独立储能系统等应用领域中得到了广泛的应用。光耦合器,以其高速信号传输、出色的共模抑制比以及单向信号传输和光电隔离的特性,在光伏逆变器中扮演着至关重要的角色。它确保了系统的安全隔离、干扰的有效隔离以及通信信号的精准传输。光耦合器的使用不仅提高了系统的稳定性和安全性,而且由于其低功耗的
    晶台光耦 2024-12-02 10:40 120浏览
  • 11-29学习笔记11-29学习笔记习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习笔记&记录学习习笔记&记学习学习笔记&记录学习学习笔记&记录学习习笔记&记录学习学习笔记&记录学习学习笔记记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&学习学习笔记&记录学习学习笔记&记录学习学习笔记&记
    youyeye 2024-12-02 23:58 71浏览
  • 当前,智能汽车产业迎来重大变局,随着人工智能、5G、大数据等新一代信息技术的迅猛发展,智能网联汽车正呈现强劲发展势头。11月26日,在2024紫光展锐全球合作伙伴大会汽车电子生态论坛上,紫光展锐与上汽海外出行联合发布搭载紫光展锐A7870的上汽海外MG量产车型,并发布A7710系列UWB数字钥匙解决方案平台,可应用于数字钥匙、活体检测、脚踢雷达、自动泊车等多种智能汽车场景。 联合发布量产车型,推动汽车智能化出海紫光展锐与上汽海外出行达成战略合作,联合发布搭载紫光展锐A7870的量产车型
    紫光展锐 2024-12-03 11:38 101浏览
我要评论
0
点击右上角,分享到朋友圈 我知道啦
请使用浏览器分享功能 我知道啦