srand(0);

//srand(time(0));

num_users = nu; num_movies = nm; num_features = nf;

movie_features = new float[(nm+nu)*nf];

user_features = movie_features + (nm*nf);

for (int i=0; i<(nm+nu)*nf; i++)

movie_features[i] = rndn();

movie_avg = new double[nu+nm];

load_averages(movie_avg);

user_avg = movie_avg + nm;

average = 3.6033;

float sum = average + movie_avg[movie] + user_avg[user];

sum += cblas_sdot(num_features, &movie_features[movie*num_features],

1, &user_features[user*num_features], 1);

if (sum > 5) sum = 5;

if (sum < 1) sum = 1;

return sum;

return 3*((double)rand()/RAND_MAX/5 - .1);

// the below uses the boost library to get normal random numbers

// not really necessary