我现在写了一个fortran code，f90文件，算法的过程当中需要用到
1. 类似imsl的sort_real_int的函数，就是排序，然后取出排序后的向量下标，也类似R的order函数
2. 需要抽正态分布随机数
现在想用GCC编译然后用R调用，写成一个R函数
我自己写的时候是用了imsl的函数，但编译过程中google半天发现使用gcc和imsl的链接比较messy，fortran自带的又没有这两个函数的实现，不知道有没有其他解决方案呢？mkl可以吗？
跪求指点！