Java CommonOps.transpose方法代码示例

import org.ejml.ops.CommonOps; //导入方法依赖的package包/类
private MatrixFormulation() {
    int numRows = response.length;
    int numCols = predictors.length + ((hasIntercept) ? 1 : 0);
    this.X = createMatrixA(numRows, numCols);
    this.Xt = new DenseMatrix64F(numCols, numRows);
    CommonOps.transpose(X, Xt);
    this.XtXInv = new DenseMatrix64F(numCols, numCols);
    this.b = new DenseMatrix64F(numCols, 1);
    this.y = new DenseMatrix64F(numRows, 1);
    solveSystem(numRows, numCols);
    this.fitted = computeFittedValues();
    this.residuals = computeResiduals();
    this.sigma2 = estimateSigma2(numCols);
    this.covarianceMatrix = new DenseMatrix64F(numCols, numCols);
    CommonOps.scale(sigma2, XtXInv, covarianceMatrix);
}
 

import org.ejml.ops.CommonOps; //导入方法依赖的package包/类
@Test
public void canComputeChiSquared()
{
	// We have to use the transpose of the table
	DenseMatrix64F m = new DenseMatrix64F(chi2);
	CommonOps.transpose(m);
	int max = m.numCols;
	double[] et = m.data;
	for (int i = 0, j = 0; i < p.length; i++)
	{
		ChiSquaredDistributionTable upperTable = ChiSquaredDistributionTable.createUpperTailed(p[i], max);
		// Use 1-p as the significance level to get the same critical values
		ChiSquaredDistributionTable lowerTable = ChiSquaredDistributionTable.createLowerTailed(1 - p[i], max);
		for (int df = 1; df <= max; df++)
		{
			double o = upperTable.getCrititalValue(df);
			double e = et[j++];
			//System.out.printf("p=%.3f,df=%d = %f\n", p[i], df, o);
			Assert.assertEquals(e, o, 1e-2);

			// The test only stores 2 decimal places so use the computed value to set upper/lower
			double upper = o * 1.01;
			double lower = o * 0.99;

			Assert.assertTrue("Upper did not reject higher", upperTable.reject(upper, df));
			Assert.assertFalse("Upper did not reject actual value", upperTable.reject(o, df));
			Assert.assertFalse("Upper did not accept lower", upperTable.reject(lower, df));

			Assert.assertTrue("Lower did not reject lower", lowerTable.reject(lower, df));
			Assert.assertFalse("Lower did not accept actual value", lowerTable.reject(o, df));
			Assert.assertFalse("Lower did not accept higher", lowerTable.reject(upper, df));
		}
	}
}
 

import org.ejml.ops.CommonOps; //导入方法依赖的package包/类
public static void main(String[] args) {
    DenseMatrix64F y = new DenseMatrix64F(1, 2, true, new double[]{1.0, 2.0});
    DenseMatrix64F yT = CommonOps.transpose(y, null);
    DenseMatrix64F m = new DenseMatrix64F(2,2);
    CommonOps.mult(yT, y, m);

    SimpleMatrix m2 = new SimpleMatrix(2,2, true, new double[]{1,2,3,4});
    org.ejml.factory.SingularValueDecomposition<DenseMatrix64F> svd = DecompositionFactory.svd(2, 2, false, true, false);
    svd.decompose(m);
    m = svd.getV(null, true);

}
 

import org.ejml.ops.CommonOps; //导入方法依赖的package包/类
@Override
public void process(double[] input, double[] output) {
	RodriguesRotationJacobian g = new RodriguesRotationJacobian();
	
	g.process(input[0],input[1],input[2]);

	DenseMatrix64F J = DenseMatrix64F.wrap(3,9,output);
	
	System.arraycopy(g.Rx.data,0,output,0,9);
	System.arraycopy(g.Ry.data,0,output,9,9);
	System.arraycopy(g.Rz.data,0,output,18,9);

	CommonOps.transpose(J);
}
 

import org.ejml.ops.CommonOps; //导入方法依赖的package包/类
public void transpose() {
    CommonOps.transpose(this.data);
}
 

import org.ejml.ops.CommonOps; //导入方法依赖的package包/类
public Matrix transpose() {
	DenseMatrix64F ret = new DenseMatrix64F(matrix.numCols, matrix.numRows);
	CommonOps.transpose(matrix, ret);
	return new EJMLDenseDoubleMatrix2D(ret);
}
 

import org.ejml.ops.CommonOps; //导入方法依赖的package包/类
/**
 * Since the sign of the homography is ambiguous a point is required to make sure the correct
 * one was selected.
 *
 * @param p test point, used to determine the sign of the matrix.
 */
protected void adjustHomographSign( PairLineNorm p , DenseMatrix64F H ) {

	CommonOps.transpose(H,H_t);

	double val = GeometryMath_F64.innerProd(p.l1, H_t, p.l2);

	if( val < 0 )
		CommonOps.scale(-1, H);
}