<< qr QD: 32 bytes decimal arithmetics qdzeros >>

DD_QD (MuPAT) >> DD_QD (MuPAT) > QD: 32 bytes decimal arithmetics > qdrand

qdrand

(Quasi) Pseudorandom QD number generator

Syntax

a = qdrand(m,n)
a = rand(QDmat)

Parameters

m,n

integers

QDmat

a matrix of QD numbers

a

matrix of size [m,n] or size(QDmat) of random QD numbers, with uniform distribution on [0,1[.

Description

qdrand() uses the rand() generator. To set the generator, please refer to the rand() page.

Examples

qdrand(1,1)    // scalar
qdrand(1,4)    // vector
qdrand(2,3)    // matrix

q = d2qd([1 2 3 ; 4 5 6]);  // template
rand(q)
--> qdrand(1,1)    // scalar
 ans  =
   2.20740856602787985822517890483961187853820925130362879048353288E-1

--> qdrand(1,4)    // vector
 ans  =
[d1]
   0.6561381   0.2445539   0.5283124   0.8468926

[d2] 10^-17 *
   1.0506091  -1.2569616   1.1229981   1.223734

[d3] 10^-34 *
  -5.6297094   3.1317322   0.5624801   2.1559544

[d4] 10^-51 *
  -31.699615  -7.739013   4.5088728   4.9954088


--> qdrand(2,3)    // matrix
 ans  =
[d1]
   0.3161073   0.5715175   0.824862
   0.5305191   0.0478015   0.5798843

[d2] 10^-17 *
   0.1625074   1.8937044   0.6539004
  -4.4689629   0.2948024   3.7193957

[d3] 10^-34 *
  -0.9364927   4.6887237  -2.8213599
   11.553239   0.6494335   8.5577923

[d4] 10^-51 *
  -1.3157813   25.767596  -17.926795
  -37.604663   2.6530177  -85.1845


--> q = d2qd([1 2 3 ; 4 5 6]);  // template
--> rand(q)
 ans  =
[d1]
   0.2039064   0.0181815   0.0105835
   0.158999    0.4098371   0.196531

[d2] 10^-17 *
  -0.4996219  -0.1327166   0.0438558
   1.0698596   2.3638705  -1.0858258

[d3] 10^-34 *
   1.9138771   0.6324264  -0.1820345
   0.1360263   0.7458167  -5.6961792

[d4] 10^-51 *
  -6.67334    -3.2283822   1.2898989
   0.1737281  -2.4198306   9.7320647
 

See Also

Authors

Copyright (C) 2011 - Tsubasa Saito

Report an issue
<< qr QD: 32 bytes decimal arithmetics qdzeros >>