Hello Every1,
I have problem with pointers in my code.
I want to simulate the array according to the probabilities of the elements in the array.
I am using GSL library for that.
The syntax it uses is
I dont know how i should pass the probabilities from the arrays storing them.Code:gsl_ran_poisson (r, prob); // r is the envi setup ,prob is the probability for the Bernoulli trail
I have 3D arrays with values.
What data types are r and prob? Are they the pointers you're having problems with?
Do you have a function prototype or definition for this? Any return type on this function?
We need all the details!
r has nothing to do here, its just to set up environment for gsl lib.. (Random Number Distribution Examples - GNU Scientific Library -- Reference Manual)
prob is floating value....
it should return me array of integers
which it doesnt do
Read the docs - it doesn't return an array of integers. You need to fill the array by calling the function in a loop...Originally Posted by satty
r has nothing to do here, its just to set up environment for gsl lib.. (Random Number Distribution Examples - GNU Scientific Library -- Reference Manual)
prob is floating value....
it should return me array of integers
which it doesnt do
Code:#include <cmath> #include <complex> bool euler_flip(bool value) { return std::pow ( std::complex<float>(std::exp(1.0)), std::complex<float>(0, 1) * std::complex<float>(std::atan(1.0) *(1 << (value + 2))) ).real() < 0; }
I meant r doesnt have anything to do with the issue i have....
also i know i shd use for loop....
dat anyway I am using....
my question was how can i pass 3D array values to prob for gsl function
It's not entirely clear what you are trying to do. It sounds like you need to cycle through the array, passing each element as the second parameter of gsl_ran_poisson, storing the result somewhere, obviously. Anyway, it might help if you gave a more complete description of the problem...
Code:#include <cmath> #include <complex> bool euler_flip(bool value) { return std::pow ( std::complex<float>(std::exp(1.0)), std::complex<float>(0, 1) * std::complex<float>(std::atan(1.0) *(1 << (value + 2))) ).real() < 0; }
array1 has floating values, that i want to use as probabilites for my program.Code:for (i=0; i<iter; i++){ for (j=0; j<nodes; j++){ for (k=0; k<nodes; k++){ array2[i][j][k]=0; prob = array1[j][k]; ????? array2[i][j][k] = gsl_ran_bernoulli(r, prob);
but it is a 2D array, i dont know how shd i pass it to prob....
Without knowing the exact implementation of your program, it's impossible to say. What do the arrays represent, and what is the significance of the dimensions being used?
Code:#include <cmath> #include <complex> bool euler_flip(bool value) { return std::pow ( std::complex<float>(std::exp(1.0)), std::complex<float>(0, 1) * std::complex<float>(std::atan(1.0) *(1 << (value + 2))) ).real() < 0; }
array1 just has values from a file which looks like this
## node1 node1 value
n array 2 has to be one dimension higher than the array1 one because i want to iterate it for number of times...for bernouli trial
In that case, it seems like the code you posted may be the correct scheme. Not 100% sure, though. Have you run it to see if it gives satisfactory results?array1 just has values from a file which looks like this
## node1 node1 value
n array 2 has to be one dimension higher than the array1 one because i want to iterate it for number of times...for bernouli trial
Code:#include <cmath> #include <complex> bool euler_flip(bool value) { return std::pow ( std::complex<float>(std::exp(1.0)), std::complex<float>(0, 1) * std::complex<float>(std::atan(1.0) *(1 << (value + 2))) ).real() < 0; }
Post the prototype of gsl_ran_poisson() so it's clear what are its arguments and return values.
