ZVEC *zv_rand(ZVEC *x) #endif { if ( ! x ) error(E_NULL,"zv_rand"); mrandlist((Real *)(x->ve),2*x->dim); return x; }
/* v_rand -- initialises x to be a random vector, components independently & uniformly ditributed between 0 and 1 */ VEC *v_rand(VEC *x) { /* int i; */ if ( ! x ) error(E_NULL,"v_rand"); /* for ( i = 0; i < x->dim; i++ ) */ /* x->ve[i] = rand()/((Real)MAX_RAND); */ /* x->ve[i] = mrand(); */ mrandlist(x->ve,x->dim); return x; }
/* m_rand -- initialises A to be a random vector, components independently & uniformly distributed between 0 and 1 */ MAT *m_rand(MAT *A) { int i /* , j */; if ( ! A ) error(E_NULL,"m_rand"); for ( i = 0; i < A->m; i++ ) /* for ( j = 0; j < A->n; j++ ) */ /* A->me[i][j] = rand()/((Real)MAX_RAND); */ /* A->me[i][j] = mrand(); */ mrandlist(A->me[i],A->n); return A; }