GSLライブラリを使用してCコードを記述しています。私はGSLライブラリを使って自分のコードに変数を割り振り、解放し、割り振る際に、ある種の矛盾（少なくともCの限られた知識から）に気付いています。最初の割り当てがループ内にあるときは、うまく動作します。別の関数で int i;
for(i=1; i < 101; i++){
gsl_matrix * W = gsl_matrix_all
`で正常に動作しません：いくつかの方法で # include <gsl/gsl_statistics_double.h>
# include <iostream>
using namespace std;
int main()
{
double y = 50.2944, yc = 63.2128;
double pearson_corr = gsl_stats_co
this postを見て、Rcppで行列をサブセット化しようとしました。 RcppArmadilloで ： // [[Rcpp::depends(RcppArmadillo)]]
#include "RcppArmadillo.h"
// [[Rcpp::export]]
arma::mat submatrix(const arma::mat& m1in, int fromin, int t
私は明らかによく知られている問題に苦しんでいます。クラスメンバー関数で定義されたODEシステムがあり、GSLソルバの1つで解決/統合したいと考えています。私の統合ルーチンに続いて class my_model{
...
public:
int NEQ = 4;
double y[4], dydt[4];
double params[25];
int od