/// form (q transpose)*fvec and store the first n components in qtf.
form (q transpose)*fvec and store the first n components in qtf.
form ( q transpose ) * fvec and store the first n components in qtf .
fn transpose(&mut self) {
        self.wa4.copy_from_slice(&self.fvec);
        let mut jj = 0;
        for j in 0..self.nfree {
            let temp = self.fjac[jj];
            if temp != 0. {
                let mut sum = 0.0;
                let mut ij = jj;
                for i in j..self.m {
                    sum += self.fjac[ij] * self.wa4[i];
                    ij += 1;
                }
                let temp = -sum / temp;
                ij = jj;
                for i in j..self.m {
                    self.wa4[i] += self.fjac[ij] * temp;
                    ij += 1;
                }
            }
            self.fjac[jj] = self.wa1[j];
            jj += self.m + 1;
            self.qtf[j] = self.wa4[j];
        }
    }
fn transpose ( & mut self ) { self . wa4 . copy_from_slice ( & self . fvec ) ; let mut jj = 0 ; for j in 0 .. self . nfree { let temp = self . fjac [ jj ] ; if temp != 0. { let mut sum = 0.0 ; let mut ij = jj ; for i in j .. self . m { sum += self . fjac [ ij ] * self . wa4 [ i ] ; ij += 1 ; } let temp = - sum / temp ; ij = jj ; for i in j .. self . m { self . wa4 [ i ] += self . fjac [ ij ] * temp ; ij += 1 ; } } self . fjac [ jj ] = self . wa1 [ j ] ; jj += self . m + 1 ; self . qtf [ j ] = self . wa4 [ j ] ; } }
form (q transpose)*fvec and store the first n components in qtf.
form ( q transpose ) * fvec and store the first n components in qtf .
&mut self