Commit 8e812664 authored by Alois SCHLOEGL's avatar Alois SCHLOEGL

fix 64bit index, that was represented only by int (32bit)

parent 0fe5e345
...@@ -58,7 +58,7 @@ int calc_N_motifs(gsl_matrix_float *W, int N_nodes, ...@@ -58,7 +58,7 @@ int calc_N_motifs(gsl_matrix_float *W, int N_nodes,
for(int i=0; i<N_nodes;i++) for(int i=0; i<N_nodes;i++)
row_sums[i]=column_sums[i]=0.0; row_sums[i]=column_sums[i]=0.0;
for(int i=0; i<N_nodes;i++) { for(int i=0; i<N_nodes;i++) {
int i_tda = i*tda; size_t i_tda = i*tda;
for(int j=0; j<N_nodes;j++) { for(int j=0; j<N_nodes;j++) {
double temp=thedata[i_tda+j]; double temp=thedata[i_tda+j];
row_sums[i]+=temp; row_sums[i]+=temp;
......
...@@ -474,7 +474,7 @@ int gen_corr_gaussian(int N_pops, int *N_nodes, ...@@ -474,7 +474,7 @@ int gen_corr_gaussian(int N_pops, int *N_nodes,
for(int pop1=0; pop1<N_pops; pop1++) for(int pop1=0; pop1<N_pops; pop1++)
for(int i=N_shift[pop1]; i<N_shift[pop1+1]; i++) { for(int i=N_shift[pop1]; i<N_shift[pop1+1]; i++) {
int i_tda=i*tda; size_t i_tda=i*tda;
for(int pop2=0; pop2<N_pops; pop2++) for(int pop2=0; pop2<N_pops; pop2++)
for(int j=N_shift[pop2]; j<N_shift[pop2+1]; j++) { for(int j=N_shift[pop2]; j<N_shift[pop2+1]; j++) {
......
...@@ -246,7 +246,7 @@ int gen_corr_gaussian(const int N_nodes, double sqrt_diag, double sqrt_recip, ...@@ -246,7 +246,7 @@ int gen_corr_gaussian(const int N_nodes, double sqrt_diag, double sqrt_recip,
double matrix_sum = 0.0; double matrix_sum = 0.0;
for (int i=0; i<N_nodes;i++){ for (int i=0; i<N_nodes;i++){
int i_tda=i*tda; size_t i_tda=i*tda;
for (int j=0; j<N_nodes;j++){ for (int j=0; j<N_nodes;j++){
// no connection from node onto itself // no connection from node onto itself
if(j==i) if(j==i)
...@@ -272,7 +272,7 @@ int gen_corr_gaussian(const int N_nodes, double sqrt_diag, double sqrt_recip, ...@@ -272,7 +272,7 @@ int gen_corr_gaussian(const int N_nodes, double sqrt_diag, double sqrt_recip,
matrix_sum *= 0.5; matrix_sum *= 0.5;
for (int i=0; i<N_nodes; i++){ for (int i=0; i<N_nodes; i++){
int i_tda=i*tda; size_t i_tda=i*tda;
for (int j=0; j<N_nodes; j++){ for (int j=0; j<N_nodes; j++){
// no connection from node onto itself // no connection from node onto itself
if(i==j) if(i==j)
......
...@@ -467,7 +467,7 @@ int gen_corr_gaussian(const int N_nodes[2], ...@@ -467,7 +467,7 @@ int gen_corr_gaussian(const int N_nodes[2],
for(int ntype1=0; ntype1<2; ntype1++) for(int ntype1=0; ntype1<2; ntype1++)
for(int i=N_shift[ntype1]; i<N_nodes[ntype1]+N_shift[ntype1]; i++) { for(int i=N_shift[ntype1]; i<N_nodes[ntype1]+N_shift[ntype1]; i++) {
int i_tda=i*tda; size_t i_tda=i*tda;
for(int ntype2=0; ntype2<2; ntype2++) for(int ntype2=0; ntype2<2; ntype2++)
for(int j=N_shift[ntype2]; j<N_nodes[ntype2]+N_shift[ntype2]; j++) { for(int j=N_shift[ntype2]; j<N_nodes[ntype2]+N_shift[ntype2]; j++) {
...@@ -497,7 +497,7 @@ int gen_corr_gaussian(const int N_nodes[2], ...@@ -497,7 +497,7 @@ int gen_corr_gaussian(const int N_nodes[2],
for(int ntype1=0; ntype1<2; ntype1++) for(int ntype1=0; ntype1<2; ntype1++)
for(int i=N_shift[ntype1]; i<N_nodes[ntype1]+N_shift[ntype1]; i++) { for(int i=N_shift[ntype1]; i<N_nodes[ntype1]+N_shift[ntype1]; i++) {
int i_tda=i*tda; size_t i_tda=i*tda;
for(int ntype2=0; ntype2<2; ntype2++) for(int ntype2=0; ntype2<2; ntype2++)
for(int j=N_shift[ntype2]; j<N_nodes[ntype2]+N_shift[ntype2]; j++) { for(int j=N_shift[ntype2]; j<N_nodes[ntype2]+N_shift[ntype2]; j++) {
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment