Commit 5832df5

committed

avoid 128-bit random numbers

1 parent a05809e commit 5832df5Copy full SHA for 5832df5

File tree

-3

lines changed

-3

lines changed

Lines changed: 6 additions & 3 deletions

Original file line number	Diff line number	Diff line change
`@@ -70,14 +70,17 @@ module test_linalg_least_squares`
`70`	`70`
`71`	`71`	`type(linalg_state_type) :: state`
`72`	`72`	`integer(ilp), parameter :: n = 12, m = 3`
	`73`	`+ real :: Arnd(n,m),xrnd(m)`
`73`	`74`	`${rt}$ :: xsol(m),x(m),y(n),A(n,m)`
`74`	`75`
`75`	`76`	`! Random coefficient matrix and solution`
`76`		`- call random_number(A)`
`77`		`- call random_number(xsol)`
	`77`	`+ call random_number(Arnd)`
	`78`	`+ call random_number(xrnd)`
`78`	`79`
`79`	`80`	`! Compute rhs`
`80`		`- y = matmul(A,xsol)`
	`81`	`+ A = real(Arnd,${rk}$)`
	`82`	`+ xsol = real(xrnd,${rk}$)`
	`83`	`+ y = matmul(A,xsol)`
`81`	`84`
`82`	`85`	`! Find polynomial`
`83`	`86`	`x = lstsq(A,y,err=state)`

Comments

(0)