On 06/11/2024 12:39, Marcin Błażejowski wrote:
Hi,
as long as I understand docs, qform() is not an option if A is identity.
But what if it is diagonal:
We introduced a faster variant in 2023c exactly for this case: the
second argument can be a vector, which is understood as the diagonal of
the middle matrix. Try appending this to your script:
<hansl>
dA = diag(A)
loop N
r3 = qform(x', dA)
endloop
printf "Total time: %f\n", $stopwatch
</hansl>
-------------------------------------------------------
Riccardo (Jack) Lucchetti
Dipartimento di Scienze Economiche e Sociali (DiSES)
Università Politecnica delle Marche
(formerly known as Università di Ancona)
r.lucchetti(a)univpm.it
http://www2.econ.univpm.it/servizi/hpp/lucchetti
-------------------------------------------------------