Fortran Coder

标题: mkl能求解非线性方程组吗? [打印本页]

作者: lwq1834515    时间: 2015-4-10 17:20
标题: mkl能求解非线性方程组吗?
就是类似于imsl的NEQNF()函数。刚开始学习FORTRAN,mkl资料好少啊。

作者: lwq1834515    时间: 2015-4-11 08:41
大神回答一下嘛……
作者: fcode    时间: 2015-4-11 12:23
既然对IMSL有足够的了解,那就使用 IMSL 就好了,为何非要用MKL?
作者: lwq1834515    时间: 2015-4-11 13:16
fcode 发表于 2015-4-11 12:23
既然对IMSL有足够的了解,那就使用 IMSL 就好了,为何非要用MKL?

我看了一些fortran的教程,了解了一点IMSL,可是IMSL一直配置的有问题。才打算用MKL的。
作者: fcode    时间: 2015-4-11 14:25
配置是个容易解决的事情。相比而言,我觉得 MKL 的配置更难。

我没太多用过 MKL,解非线性方程组,应该是用雅克比相关的函数的。你可能需要若干个函数配合才能实现(当然,这样更灵活)但不如 IMSL 快捷。
作者: lwq1834515    时间: 2015-4-11 18:00
本帖最后由 lwq1834515 于 2015-4-11 18:28 编辑
fcode 发表于 2015-4-11 14:25
配置是个容易解决的事情。相比而言,我觉得 MKL 的配置更难。

我没太多用过 MKL,解非线性方程组,应该是 ...

我在配置的时候出现了这样的错误:
Error        1         error #6866: Dotted string neither a defined operator nor a structure component   [IX]        
Error        2         Compilation Aborted (code 1)        
这是代码:
                                                                                                                                                                                                                              
program main

  include 'link_fnl_shared.h'
  use numerical_libraries
  implicit none

  real::a(2,2)=(/1,2,12,3/)
  real::b(2,2)=(/2,1,3,4/)
  real::c(2,2)
  
  c=a .ix.b
  write(*,*) c

  stop
end program main

如果把 use numerical_libraries改成use imsl_libraries,就出现这样的错误:
Error        1        Compilation Aborted (code 1)       

只有这一个。

是不是我配置的不对?


作者: fcode    时间: 2015-4-11 19:24
加一句 use linear_operators
作者: lwq1834515    时间: 2015-4-11 20:32
fcode 发表于 2015-4-11 19:24
加一句 use linear_operators

加上了还是说 Compilation Aborted,只加use linear_operators也是这样。

1.png (21.41 KB, 下载次数: 398)

1.png

作者: lwq1834515    时间: 2015-4-11 20:39
fcode 发表于 2015-4-11 19:24
加一句 use linear_operators

这个是output的内容:

1.png (13.84 KB, 下载次数: 382)

1.png

作者: fcode    时间: 2015-4-11 21:28
出现这种可能一般是,你的系统是32位,而下载的 IMSL 是64位。或者其他参数不匹配。
作者: lwq1834515    时间: 2015-4-12 08:31
这样啊,那只能重新安装了
谢谢耐心指教




欢迎光临 Fortran Coder (http://bbs.fcode.cn/) Powered by Discuz! X3.2