A******g
发帖数: 612 | 1 请教各位，
想解这么个简单的问题
Ax=y
A是4x4矩阵，x，y是4x1矢量
已知x，y想求解A
matlab里很简单，A=y\x
请问有没有方便使用的，免费的C＋＋库可以做这个。懒得自己implement了，太费劲。
最好有例子。
以前看过linpack的doc，好像没看懂，放弃了，觉得讲了半天不知所云，也没什么例子。
谢谢了！ |
y***d
发帖数: 2330 | 2 gsl
子。
【在 A******g 的大作中提到】
: 请教各位，
: 想解这么个简单的问题
: Ax=y
: A是4x4矩阵，x，y是4x1矢量
: 已知x，y想求解A
: matlab里很简单，A=y\x
: 请问有没有方便使用的，免费的C＋＋库可以做这个。懒得自己implement了，太费劲。
: 最好有例子。
: 以前看过linpack的doc，好像没看懂，放弃了，觉得讲了半天不知所云，也没什么例子。
: 谢谢了！
|
A******g
发帖数: 612 | 3 谢谢，请问gsl里面哪个function可以解出A，我去查相应的doc
【在 y***d 的大作中提到】
: gsl
:
: 子。
|
y***d
发帖数: 2330 | |
a****l
发帖数: 8211 | 5 I think matlab has free api interface for C? At least I used it some 10
years ago.
子。
【在 A******g 的大作中提到】
: 请教各位，
: 想解这么个简单的问题
: Ax=y
: A是4x4矩阵，x，y是4x1矢量
: 已知x，y想求解A
: matlab里很简单，A=y\x
: 请问有没有方便使用的，免费的C＋＋库可以做这个。懒得自己implement了，太费劲。
: 最好有例子。
: 以前看过linpack的doc，好像没看懂，放弃了，觉得讲了半天不知所云，也没什么例子。
: 谢谢了！
|
t****t
发帖数: 6806 | 6 the support is discontinued since they want to promote their own compiler (
from .m ).
【在 a****l 的大作中提到】
: I think matlab has free api interface for C? At least I used it some 10
: years ago.
:
: 子。
|
r*******n
发帖数: 3020 | 7 scilab?
子。
【在 A******g 的大作中提到】
: 请教各位，
: 想解这么个简单的问题
: Ax=y
: A是4x4矩阵，x，y是4x1矢量
: 已知x，y想求解A
: matlab里很简单，A=y\x
: 请问有没有方便使用的，免费的C＋＋库可以做这个。懒得自己implement了，太费劲。
: 最好有例子。
: 以前看过linpack的doc，好像没看懂，放弃了，觉得讲了半天不知所云，也没什么例子。
: 谢谢了！
|
p***o
发帖数: 1252 | 8 太dt了，4*4的方程，不用loop估计100行也能解出来。
不过再看一遍lz的问题，似乎还不是解方程 ...
【在 a****l 的大作中提到】
: I think matlab has free api interface for C? At least I used it some 10
: years ago.
:
: 子。
|
s****a
发帖数: 238 | 9 I'd rather use blitz++ as the container, which claims to be as fast as
fortran, and pull out the pointer to data with data() so that gsl or fftw
can use them.
You can also write a wrapper for gsl or fftw, which makes the usage look
more like matlab
【在 y***d 的大作中提到】
: gsl
:
: 子。
|
p**n
发帖数: 1437 | |
|
|
a****l
发帖数: 8211 | 11 this is mean....
【在 t****t 的大作中提到】
: the support is discontinued since they want to promote their own compiler (
: from .m ).
|
r*********r
发帖数: 3195 | |
l**0
发帖数: 23 | 13 确认是 知道 x, y 求 A吗？ 这样是个不定方程阿，4个方程，16个未知量，一般会有
无穷个解。 |
A******g
发帖数: 612 | 14 谢谢，再次看还是觉得麻烦
用了一下这个例子，
总是有很多 #include"gsl/gsl_***"
这个前面的 gsl/ 总是让找不到.h file,改了几个以后放弃了
一怒之下，自己写了个class做矩阵计算，还不错
用C++就是没得偷懒啊...
要不花很多时间学别人的，然后不停改设置
要不自己写
如果这个数学的东西很成为标准C++的内容就好了
Examples.html
【在 y***d 的大作中提到】
: http://www.gnu.org/software/gsl/manual/html_node/Linear-Algebra-Examples.html
|
A******g
发帖数: 612 | 15 你说得对，昨天太困没想清楚问题
睡醒想到别的办法了
【在 l**0 的大作中提到】
: 确认是 知道 x, y 求 A吗？ 这样是个不定方程阿，4个方程，16个未知量，一般会有
: 无穷个解。
|
A******g
发帖数: 612 | 16 Yes,I know how to compile matlab code into a .dll and the function can
thus be called from C/C++
However,it takes some time to wrap matlab matrix using the class they
give. This takes extra effort. And it is harder to debug.
Anyway, I wrote my own C++ class to handle matrix now.
Thanks a lot!
【在 a****l 的大作中提到】
: I think matlab has free api interface for C? At least I used it some 10
: years ago.
:
: 子。
|
A******g
发帖数: 612 | 17 I was just lazy...
This is just a small part of the bigger problem... wanna take a shortcut
if there is one...
【在 p***o 的大作中提到】
: 太dt了，4*4的方程，不用loop估计100行也能解出来。
: 不过再看一遍lz的问题，似乎还不是解方程 ...
|
p*********g
发帖数: 226 | 18 这个链接讲如何在C/C++ 中调用matlab库
http://webdocs.cs.ualberta.ca/~xinhua2/matlab.htm
里面的 PETSc 和 TAO 可以略过。
【在 a****l 的大作中提到】
: I think matlab has free api interface for C? At least I used it some 10
: years ago.
:
: 子。
|
l*****z
发帖数: 2305 | 19 intel mkl (under linux, free) |
H***a
发帖数: 735 | 20 还是用库吧, 别的没用过, blas/lapack很好用了. 自己写的话碰到逆阵,稀疏矩阵你的
算法多半不是优化的. |
|
|
G*****7
发帖数: 1759 | 21 4x4 -> use intel ipp smal matrix library
4k x 4k ->use intel mkl blas/linpack/sparse solvers
either way, intel ipp/mkl is free for academic use under linux and comes
with tons of examples.
子。
【在 A******g 的大作中提到】
: 请教各位，
: 想解这么个简单的问题
: Ax=y
: A是4x4矩阵，x，y是4x1矢量
: 已知x，y想求解A
: matlab里很简单，A=y\x
: 请问有没有方便使用的，免费的C＋＋库可以做这个。懒得自己implement了，太费劲。
: 最好有例子。
: 以前看过linpack的doc，好像没看懂，放弃了，觉得讲了半天不知所云，也没什么例子。
: 谢谢了！
|
G*****7
发帖数: 1759 | 22 they still have that interface. engine.h
you can evaluate matlab commands from c/cpp.
【在 a****l 的大作中提到】
: I think matlab has free api interface for C? At least I used it some 10
: years ago.
:
: 子。
|
G*****7
发帖数: 1759 | 23 this is called least square problem
AX=Y -> A'AX=A'Y -> X= inverse(A'A)A'Y -> in matlab X = pinv(A)*Y.
子。
【在 A******g 的大作中提到】
: 请教各位，
: 想解这么个简单的问题
: Ax=y
: A是4x4矩阵，x，y是4x1矢量
: 已知x，y想求解A
: matlab里很简单，A=y\x
: 请问有没有方便使用的，免费的C＋＋库可以做这个。懒得自己implement了，太费劲。
: 最好有例子。
: 以前看过linpack的doc，好像没看懂，放弃了，觉得讲了半天不知所云，也没什么例子。
: 谢谢了！
|
w***g
发帖数: 5958 | 24 这个我研究过一点点.我觉得没什么懒可偷的.矩阵运算要做得快就得老老实实学BLAS和
LAPACK.因为那个是industry standard,Intel提供的library用的也都是这些接口.C++那
些都是奇技淫巧.
【在 s****a 的大作中提到】
: I'd rather use blitz++ as the container, which claims to be as fast as
: fortran, and pull out the pointer to data with data() so that gsl or fftw
: can use them.
: You can also write a wrapper for gsl or fftw, which makes the usage look
: more like matlab
|
f*******a
发帖数: 663 | 25 Matlab 6.5 下带的C/C++ 库其实非常好用，虽然版本旧了点，但做矩阵运算绰绰有余
。但正是因为这个库的强大，Mathwork从7开始就完全改变了接口。
有兴趣的话可以去查一下，上手要点时间，但很方便。缺点是以后的程序都得带一堆
dll。 |
t****t
发帖数: 6806 | 26 那个C++版本的库其实一般好用, 语法不是很intuitive
C版本的就更惨了
我曾经重写了一个比较好的. 不过后来我做矩阵很少, 就不用了.
【在 f*******a 的大作中提到】
: Matlab 6.5 下带的C/C++ 库其实非常好用，虽然版本旧了点，但做矩阵运算绰绰有余
: 。但正是因为这个库的强大，Mathwork从7开始就完全改变了接口。
: 有兴趣的话可以去查一下，上手要点时间，但很方便。缺点是以后的程序都得带一堆
: dll。
|
f*******a
发帖数: 663 | 27 我倒没觉得怎么不好用，熟悉了设置使用也还蛮简单的。复合运算时也不需要什么中间
变量，编写效率挺高的。当然，C库使用是非常的不便，但个人觉得C++库已经封装的挺
好了。不知道这位兄台做了什么样的改进？
【在 t****t 的大作中提到】
: 那个C++版本的库其实一般好用, 语法不是很intuitive
: C版本的就更惨了
: 我曾经重写了一个比较好的. 不过后来我做矩阵很少, 就不用了.
|
s****a
发帖数: 238 | 28 老板不喜欢不开源的...blas我看过一阵，dsgmm什么的太吓人了，自己用总还要写
wrapper吧，不过还是谢谢指点，不能太依赖这些了。
+那
【在 w***g 的大作中提到】
: 这个我研究过一点点.我觉得没什么懒可偷的.矩阵运算要做得快就得老老实实学BLAS和
: LAPACK.因为那个是industry standard,Intel提供的library用的也都是这些接口.C++那
: 些都是奇技淫巧.
|
