将C++代码转换为 R 代码以生成数据

Question

我希望根据以下代码为我的硕士论文创建一些数据，但我只熟悉 R，有关以下代码的作用或在 R 中实现它的一些指导将不胜感激。

case 9: // generates a 20-dim affine space  in 100- dimensions
dim=20;
// allocate memory for the output
plhs[0] = mxCreateDoubleMatrix(dim, num, mxREAL);
data = mxGetPr(plhs[0]);
para=new double*[20];
for(i=0;i<20;i++)
para[i]=new double[long(num)];
for(i=0;i<20;i++)
for(j=0;j<num;j++)
para[i][j]=5.0*uni()-2.5;
vectors=new double*[20];
for(i=0;i<20;i++)
vectors[i]=new double[dim];
for(i=0;i<20;i++)
for(j=0;j<dim;j++)
{
if(i==j)
vectors[i][j]=1;
else
vectors[i][j]=0;//0.25*(uni()-0.125);
}
//vectors[i][j]=10*(uni()-0.5);
//generate the figure
for(i=0;i<num;i++)
{
for(j=0;j<dim;j++)
{
data[i*dim+j]=0;
for(k=0;k<20;k++)
data[i*dim+j]+=para[k][i]*vectors[k][j];
}
}
mexPrintf("Generated %i data points of a 20-dimensional affine subspace in R^20n",num);
mexPrintf("The correct dimension of this submanifold is 20n");
for(i=0;i<20;i++)
delete para[i];
delete para;
break;

num是我想要生成的点数，dim维度，因此它们本质上是我的输入。

我编辑了我的帖子，因为我看到需要整个代码才能正确回答，很抱歉给您带来不便。

Answer 1

这里有一个翻译，也许你可以尝试一下

vectors <- matrix(nrow = 20,ncol = dim)
for (i in 1:20) {
for (j in 1:dim) {
vectors[i,j] <- ifelse(i==j,1,0)
}
} 
for (i in 1:20) {
for (j in 1:dim) {
data[i*(dim-1+j)] <- 0
for (k in 1:20) {
data[i*(dim-1)+j] <- data[i*(dim-1)+j] + para[k,i]*vectors[k,j]
}
}
}