[英]rcpp - R session Aborted
library(Rcpp)
library(geoR)
elevationd=as.matrix(data.frame(xcoords=elevation$coords[,1],
ycoords=elevation$coords[,2], elev=elevation$data))
elevationd
cppFunction('void a(NumericMatrix data){
int nr = data.nrow();
int nc = data.ncol();
NumericVector tmp;
for (int i; i<nr; i++){
tmp[i] = data(i,2);
}
NumericMatrix mat(nr, nr);
for (int i; i<nr; i++){
for (int j; j<nr; j++){
mat(i,j) = (tmp[i] - tmp[j])*(tmp[i] - tmp[j]);
}
}
}')
a(elevationd)
我有这样的AR代码。
但是当我实现此操作时,我收到“ R会话已中止。R遇到致命错误。会话已终止。” 信息。
但我找不到我的代码有任何问题。
我想知道这是否是硬件问题.....不是我的代码TT
ㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡㅡ
到德克
这是我编辑的代码。
cppFunction('NumericMatrix f1(NumericMatrix data){
int nr = data.nrow();
int nc = data.ncol();
NumericMatrix mat(nr, nr);
for (int i=0; i<nr; i++){
for(int j=0; j<nr; j++){
mat(i,j) = (data(i,2)-data(j,2))*(data(i,2)-data(j,2));
}
}
return mat;
}')
罗兰(Roland)已经给了您一个很好的提示,您也有两个错误的类型,并错过了矢量初始化的大小。
以下是修复的版本,作为C ++文件,您可以提供源代码并自动执行相关的R代码:
R> Rcpp::sourceCpp("/tmp/jiwon.cpp")
R> library(geoR)
--------------------------------------------------------------
Analysis of Geostatistical Data
For an Introduction to geoR go to http://www.leg.ufpr.br/geoR
geoR version 1.7-5.2 (built on 2016-05-02) is now loaded
--------------------------------------------------------------
R> elevationd <- as.matrix(data.frame(xcoords=elevation$coords[,1],
+ ycoords=elevation$coords[,2],
+ .... [TRUNCATED]
R> str(elevationd)
num [1:52, 1:3] 0.3 1.4 2.4 3.6 5.7 1.6 2.9 3.4 3.4 4.8 ...
- attr(*, "dimnames")=List of 2
..$ : chr [1:52] "1" "2" "3" "4" ...
..$ : chr [1:3] "xcoords" "ycoords" "elev"
R> str(a(elevationd))
num [1:52, 1:52] 0 5929 13225 32400 4900 ...
R>
这是修复的代码,解决了您的一些错误:
#include <Rcpp.h>
using namespace Rcpp;
// [[Rcpp::export]]
NumericMatrix a(NumericMatrix data){
int nr = data.nrow();
int nc = data.ncol();
NumericVector tmp(nr);
for (int i=0; i<nr; i++){
tmp[i] = data(i,2);
}
NumericMatrix mat(nr, nr);
for (int i=0; i<nr; i++){
for (int j=0; j<nr; j++){
mat(i,j) = (tmp[i] - tmp[j])*(tmp[i] - tmp[j]);
}
}
return mat;
}
/*** R
library(geoR)
elevationd <- as.matrix(data.frame(xcoords=elevation$coords[,1],
ycoords=elevation$coords[,2],
elev=elevation$data))
str(elevationd)
str(a(elevationd))
*/
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.