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;
}
