d****y 发帖数: 1650 | 1 matlab的kernel是很多数学库,比如LAPACK,BLAS等,这些是没有做多核优化的,你可
以考虑装上MKL,然后让matlab使用MKL做数学计算。 |
|
s***e 发帖数: 403 | 2 看出我是谁的给点面子,就别说破了。如果哪位师兄师姐师弟师妹有个什么软件公司的
码工内推,不妨给我试试,应该不会丢大家的脸。
干了几年,老板想让我拿硕士毕业,貌似不是很想给我博士学位。实验室里做的东西基
本上不适合找工作,所以现在求CS方面的机会。
技能:6+年C/C++经验,熟悉C++11,熟悉stl/boost,熟悉过程式/面向对象/模板元风
格编程,了解lambda编程,基本上可以短期轻松上手任何不是LISP系列的语言。了解
windows下的MFC开发,熟悉qt4库,熟悉linux系统和shell编程,熟悉基本算法和数据
结构,熟悉常用设计模式。有独立开发经验。有并行开发(pthread, openmp, intel
tbb)和性能优化(Vtune)的经验。能看懂汇编,但是自己写就不行了。有建模经验。
了解计算机基本架构(差不多是看完CSAPP那本书的水平)。了解一些操作系统,网络
方面的知识。还有一些相对垃圾的技能比如会用lex/yacc,会一点ruby,会用blas和
lapack,会打星际2……
地方不限,北美除了犯罪率高的地方以外,随便到哪里都可以。要求工资能养的起自己
... 阅读全帖 |
|
d*t 发帖数: 28 | 3 In lower version of MatLab, there is a function called "flops"
which counts the floating number operations.
Since version 6, because of the incorporation of LAPACK,
such a counting is no longer practical. |
|
|
w****a 发帖数: 186 | 5 Current TNT package is really weak. Don't believe it, though I'm
using it as a basis. |
|
s*******n 发帖数: 38 | 6 【 以下文字转载自 Statistics 讨论区 】
发信人: samaritan (Good+Samaritan), 信区: Statistics
标 题: C++统计library
发信站: BBS 未名空间站 (Fri Oct 19 15:38:45 2007)
要用C++编一些统计的程序。想找library。 那个比较好?gsl? Lapack? |
|
|
a****t 发帖数: 720 | 8 纯数background,没什么编程经验,现有的code是matlab的,准备改用c++,lapack。
这个问题以前我run了一些小的case,最近开会,类似的问题被n多牛人提及,于是打算
seriously的做做,看能不能form一个conjecture。
能推荐具体的cpu+motherboard配置吗?不知道i5 2320,大家觉得如何?motherboard
彻底没概念。
lot. |
|
p**o 发帖数: 3409 | 9 matlab曾经自己写过底层库,后来就直接用intel mkl里的blas和
lapack来对付稠密矩阵了。mkl是用fortran写的,也提供了c接口。
我没看过matlab的封装,不清楚它封的是哪个接口。
稀疏矩阵的话,我的理解是matlab会调umfpack来分块求解
(调amd来做preordering),其中每一块稠密矩阵又调mkl来算。
自己写/封装eig,性能超过matlab的可能性很小。
Fortran写的。 |
|
t****t 发帖数: 6806 | 10 matlab的矩阵计算核心应该本来就是公共的库, blas, lapack之类的, 不过matlab都给
你包装好了, 优化都做好, 比如说在intel的CPU上就会用MKL, 所以矩阵规模越大越合
算, 如果都是标量就慢了. |
|
p****s 发帖数: 32405 | 11 undefined reference to what?
是不是某个函数/变量没定义就在用? |
|
s****t 发帖数: 698 | 12
貌似是用了某个lib的函数,没link上吧。
如果文件不多,一个一个定位一下吧。 |
|
m*****r 发帖数: 130 | 13 你的知道是啥symbol找不到吧,然后看看在那个lib里,放到你的连接是参数就好了 |
|
t******1 发帖数: 334 | 14 谢谢楼上,已经找到答案了,就在intel mkl userguide里,我没有仔细读,有linking
这一章的。
不过话说,一周之前,我还只知道gcc *.c,什么link, flag什么都不懂,这也算是学
习的过程了。 |
|
E*V 发帖数: 17544 | 15 sudo emerge -av octave
Password:
These are the packages that would be merged, in order:
Calculating dependencies... done!
[ebuild N ] media-libs/qhull-2003.1-r1 708 kB
[ebuild N ] sci-visualization/gnuplot-4.2.5-r1 USE="X readline -doc -
emacs -gd -ggi -latex -lua -pdf -plotutils (-svga) -wxwidgets -xemacs" 2,832
kB
[ebuild N ] sci-mathematics/glpk-4.39 USE="mysql -doc -examples -gmp -
odbc" 2,390 kB
[ebuild N ] app-admin/eselect-lapack-0.1 0 kB
[ebuild N ] app-admin/esel |
|
j******u 发帖数: 1968 | 16 初学者
弄个p3的老破机器,装whatever,学习基本命令,装emule下载,偶尔写点latex,主要
是word没有办法运行。
要不就是摆酷,弄个windows下面的虚拟机,装个几百兆硬盘,呵呵
就算家里有好几台电脑,也把linux装在最破的上面。
初学者具有一定经验,我们叫初学者2。0算了
发现一些程序很好,能自动执行,能远程操控,能写程序,能干工作,可是发现主要瓶
颈是计算功能,cpu太弱,内存太小。于是升级系统,还摆脱不了windows,弄个一般一
点的电脑,把原来的windows让出一些空间,弄个双启动,没事倒腾一下grub等等
bootloader,呵呵
这个时期的主要问题停留在如何安装打印机,如何弄display或者双头显示,如何安装module。wireless等等,偶尔冒充高手编译一下核心,更多的情况下是失败,换会默认的,不吱声了。
中级用户,发现大部分工作linux都能作,而且干的很出色,而且有很多的subroutine
能调用,而且知道了很多lapack之类的库和python这种好东东。于是选择一个中高档的
电脑,安装linux,同时把windows虚拟机装上准 |
|
o**n 发帖数: 1249 | 17 我怀疑默认是不是从repository下来的都默认是不支持long double的?
怎么能简单看出已经安装的lib是否支持long double呢?多谢 |
|
m*****r 发帖数: 3822 | 18 数值计算的,用到了blas和lapack
我现在分别用了mkl,goto和atlas,只有atlas的时候是正常的。
这里面会是什么原因造成的? |
|
m*****r 发帖数: 3822 | 19 大概找到问题了。是lapack库multi-threaded的问题。
atlas我link了single-thread的lib所以没问题。现在link mkl的single library也
没问题了。不过目前goto和atlas里用single-thread都会报错,好在mkl还能用,
而且single和multi performance差不多。
thanks anyway! |
|
v*****r 发帖数: 1119 | 20 check stest.out file 看有没有动静呀。 |
|
o**o 发帖数: 3964 | 21 ATLAS is a auto-configured BLAS.
the OP actually needs a BLAS function called "xGEMM", x can be
S, D, Z, depending on data type. the function does
C := alpha * A * B + beta * C, for general matrices A B C.
LAPACK is a level higher than BLAS, probably doesn't have this
low-level arithmetic. |
|
mw 发帖数: 525 | 22 【 以下文字转载自 CS 讨论区 】
发信人: mw (mw), 信区: CS
标 题: 大家编程用得algebra lib是什么?
发信站: BBS 未名空间站 (Sat Aug 13 07:18:03 2005)
是都用lapack 吗?
还是有什么更通用的? |
|
o**o 发帖数: 3964 | 23 Lapack没有稀疏矩阵的么? 4kx4k的double complex矩阵求逆大概在半个小时到一个小时
. |
|
i******e 发帖数: 171 | 24 嗯, 如果本身用c,c++写得话,用GSL就方便.如果时间重要的话,还是要用LAPACK.
时 |
|
|
g*********e 发帖数: 42 | 26 clapack 和 lapack++
从网上down的包里有vc的工程文件,
用vc7.1编译,编译出的lib文件放在编译器的
lib目录里,dll是动态链接库,可以放在程序
目录里。 |
|
l*******r 发帖数: 322 | 27 请教哪位有经验的?
没有看到windows下编译的readme
尝试用cygwin也是一知半解
那位出手搭救?谢谢了 |
|
d***q 发帖数: 1119 | 28 标准库又不包括这种数值计算用的库
用lapack倒是不错 |
|
|
h****r 发帖数: 29 | 30 想请教一下做过Singular Value Decomposition (SVD) 的朋友。
调用Lapack的SVD子程序dgesvd。作测试时,4*80的方阵没问题,可是为什么对4*200
的方阵却出现
“ Segmentation fault”?比如,对最简单的情况,设方阵为T, T11=T22=T33=T44=
1,其余的
元素取为零。调用dgesvd却同样给出前面的错误。哪位有经验的朋友给讲讲?多谢! |
|
y****e 发帖数: 23939 | 31 debug 一下,看看到底是什么错
Fedora 8 是有些地方不好,不过我们还没有遇到BLAS的问题 |
|
k******n 发帖数: 35 | 32 I suggest to use ATLAS or ACML instead. |
|
a**a 发帖数: 416 | 33 In the level of application, you can use multi-processe to use multi-
processor.
In the other hand, if your CPU is of Intel processor, you might want to try
to use Intel's MKL as your BLAS/LAPACK engine, which might drive CPU to its
full use. In the application level, you might want to learn about the
parallel
programming paradigms, such as MPI. In the engine level, the BLAS engine
will try to take advantage of the CPU for you. |
|
r****t 发帖数: 10904 | 34 你到底看了 Makefile 没有,开头有很多你要改的东西。
/usr/local/lib 下面要有这些列出来的库,/usr/local/include 下面要有这些库的头
文件。没有的话或者按位置装好,或者装好再改这几行。
初用 make, 用 make -d 然后看看出错信息,理解不了贴上来问。你现在的问题是 gcc 在你给的 -I 选项里面列的目录里面找不到 GSL 的头文件,你就给个含有头文件的 -I path 就可以过这一步了。
# The path to the external libraries: metis, taucs, lapack, blas and zlib
LEXTPREFIX = /usr/local
LEXTLIBS = -L$(LEXTPREFIX)/lib
IEXTLIBS = -I$(LEXTPREFIX)/include
这还没列 GSL, 你要是 GSL 装在
/home/me/GSL
上面有两行你要改成
LEXTLIBS = -L$(LEXTPREFIX)/lib -L/home/me/GSL/lib
IEXTLIBS = -I$(LEXTP |
|
s*****k 发帖数: 604 | 35 我以前用过一个数值计算的C++的库,里面也用到lapack,
那里面就是假设 std::complex 和 fortan77里面的
complex16 的 memory layout是一样的。好像在大部分平台上
也确实如此。 |
|
|
s*****k 发帖数: 604 | 37 还是BLAS/LAPACK 好用,
还有针对CPU优化的版本 |
|
A**********e 发帖数: 3102 | 38 就那仨:blas,cblas,lapack。其实核心是同一个东西。 |
|
f**y 发帖数: 138 | 39 Is it easy to implement 'anova (analysis of variance)' or 'glm (general
linear model)' with those blas or lapack libs? |
|
z*****a 发帖数: 3809 | 40 Why don't you just solve for A'X' = B'?
dgetrs没
可是
是这 |
|
p*****e 发帖数: 310 | 41 因为我已经有了A的PLU,用于主要的计算,这个XA=B就是为了解另一个矩阵,它的出现
次数要比主要的计算频率少。
如果我单独为了它重新进行A的LU分解,那么和从A的PLU求inv(A)时间差不多了,除非A
'的PLU能方便的从A的PLU得到。
是发 |
|
z*****a 发帖数: 3809 | 42 dgetrs can be used to solve both AX=B and A'X=B when given the factorization
of A computed by dgetrf, you just need to pass the right flag.
非A |
|
j****i 发帖数: 305 | 43 I don't have access to the computer I was testing lapack. I can try what you
suggested and I think it might work, since
max has no "(" followed.
But what about the complex class in c++, which I cannot and do not want to
change. In that class the std::abs() function is used. There got to be a way
to solve this systematically. |
|
v****s 发帖数: 1112 | 44 oh man, can't believe u r using a single machine to compute this!!! we
usually qsub to our supercomputer with 512 nodes.
ru using LAPACK? try some sampling tek to reduce it. |
|
|
H***a 发帖数: 735 | 46 还是用库吧, 别的没用过, blas/lapack很好用了. 自己写的话碰到逆阵,稀疏矩阵你的
算法多半不是优化的. |
|
w***g 发帖数: 5958 | 47 这个我研究过一点点.我觉得没什么懒可偷的.矩阵运算要做得快就得老老实实学BLAS和
LAPACK.因为那个是industry standard,Intel提供的library用的也都是这些接口.C++那
些都是奇技淫巧. |
|
|
n**f 发帖数: 121 | 49 How do you build Lapack?
I know binaries exist for many operating systems. But I am still curious how
you build it from source, since it might be safer. Do you go through the
f2c -> gcc or other compiler path? |
|
n******t 发帖数: 4406 | 50 喀喀喀,问题是你这个问题,java 哪里比C++开发起来快了?
而且真的care 矩阵运算速度,都是用lapack之类的东西。 |
|