tests/rand/statlib.c

490   unsigned int ui_i, ui_j, ui_k, ui_l;
503   for (ui_i = 0; ui_i < GMP_SPECT_MAXT; ui_i++)
507 	  mpz_init_set_ui (U[ui_i][ui_j], 0);
508 	  mpz_init_set_ui (V[ui_i][ui_j], 0);
510       mpz_init_set_ui (X[ui_i], 0);
511       mpz_init_set_ui (Y[ui_i], 0);
512       mpz_init (Z[ui_i]);
654       for (ui_i = 0; ui_i < ui_t; ui_i++)
656 	  mpz_mul (tmp1, V[ui_i][0], r); /* tmp1=vi1*r */
659 	  mpz_sub (V[ui_i][ui_t], tmp1, tmp2);
663 	      mpz_mul (tmp1, q, U[ui_i][ui_j]);	/* tmp=q*uij */
684 	  for (ui_i = 0; ui_i <= ui_t; ui_i++)
686 	      if (ui_i != ui_j)
688 		  vz_dot (tmp1, V[ui_i], V[ui_j], ui_t + 1); /* tmp1=dot(Vi,Vj). */
705 			  mpz_sub (V[ui_i][ui_l], V[ui_i][ui_l], tmp1); /* Vi=Vi-q*Vj */
706 			  mpz_mul (tmp1, q, U[ui_i][ui_l]);
753       for (ui_i = 0; ui_i <= ui_t; ui_i++)
755 	  vz_dot (tmp1, V[ui_i], V[ui_i], ui_t + 1);
760 	  mpz_set_f (Z[ui_i], f_tmp2);
779 	      for (ui_i = 0; ui_i <= ui_t; ui_i++)
780 		mpz_add (Y[ui_i], Y[ui_i], U[ui_k][ui_i]);
787 		  for (ui_i = 0; ui_i <= ui_t; ui_i++)
789 		      mpz_mul (tmp2, tmp1, U[ui_k][ui_i]);
790 		      mpz_sub (Y[ui_i], Y[ui_i], tmp2);
812   for (ui_i = 0; ui_i < GMP_SPECT_MAXT; ui_i++)
816 	  mpz_clear (U[ui_i][ui_j]);
817 	  mpz_clear (V[ui_i][ui_j]);
819       mpz_clear (X[ui_i]);
820       mpz_clear (Y[ui_i]);
821       mpz_clear (Z[ui_i]);