Examples of Jama.Matrix

Jama.Matrix
Jama = Java Matrix class.
The Java Matrix Class provides the fundamental operations of numerical linear algebra. Various constructors create Matrices from two dimensional arrays of double precision floating point numbers. Various "gets" and "sets" provide access to submatrices and matrix elements. Several methods implement basic matrix arithmetic, including matrix addition and multiplication, matrix norms, and element-by-element array operations. Methods for reading and printing matrices are also included. All the operations in this version of the Matrix Class involve real matrices. Complex matrices may be handled in a future version.
Five fundamental matrix decompositions, which consist of pairs or triples of matrices, permutation vectors, and the like, produce results in five decomposition classes. These decompositions are accessed by the Matrix class to compute solutions of simultaneous linear equations, determinants, inverses and other matrix functions. The five decompositions are:
- Cholesky Decomposition of symmetric, positive definite matrices.
- LU Decomposition of rectangular matrices.
- QR Decomposition of rectangular matrices.
- Singular Value Decomposition of rectangular matrices.
- Eigenvalue Decomposition of both symmetric and nonsymmetric square matrices.
Example of use:
Solve a linear system A x = b and compute the residual norm, ||b - A x||.
```
 double[][] vals = {{1.,2.,3},{4.,5.,6.},{7.,8.,10.}}; Matrix A = new Matrix(vals); Matrix b = Matrix.random(3,1); Matrix x = A.solve(b); Matrix r = A.times(x).minus(b); double rnorm = r.normInf(); 
```
@author The MathWorks, Inc. and the National Institute of Standards and Technology. @version 5 August 1998

    PowerControlledCC cancellationCarrier = new PowerControlledCC(allCarriers, cancellationCarriers);    
    
    int usedSize = allCarriers-cancellationCarriers.size();
    
    for (int maxDegrees = 0; maxDegrees < cancellationCarriers.size(); maxDegrees++) {
      Matrix displayMatrix = new Matrix(usedSize, allCarriers);
      for (int basisVector = 0; basisVector < usedSize; basisVector++) {
        double[] usedData = new double[usedSize];
        usedData[basisVector] = 1;
        double[] allData = cancellationCarrier.getAssembledValuesMaxDegrees(usedData, maxDegrees);
        for (int i = 0; i < allCarriers; i++) {
          displayMatrix.set(basisVector, i, allData[i]);  
        }
      }
      System.out.println("Display Matrix");
      displayMatrix.print(15, 5);


      PdfGenerator.generatePDF(
          "CancellationCarrier_" + allCarriers + "used_" + cancellationCarriers.size() + "cc_" + maxDegrees + "degrees.pdf", 
          displayMatrix, 
          "Cancellation Carrier",

View Full Code Here

  }
      }


      // solve H d = -g, evaluate error at new location
      //double[] d = DoubleMatrix.solve(H, g);
      double[] d = (new Matrix(H)).lu().solve(new Matrix(g, nparm)).getRowPackedCopy();
      //double[] na = DoubleVector.add(a, d);
      double[] na = (new Matrix(a, nparm)).plus(new Matrix(d, nparm)).getRowPackedCopy();
      double e1 = chiSquared(x, na, y, s, f);


      if (verbose > 0) {
  System.out.println("\n\niteration "+iter+" lambda = "+lambda);
  System.out.print("a = ");
        (new Matrix(a, nparm)).print(10, 2);
  if (verbose > 1) {
          System.out.print("H = "); 
          (new Matrix(H)).print(10, 2);
          System.out.print("g = "); 
          (new Matrix(g, nparm)).print(10, 2);
          System.out.print("d = "); 
          (new Matrix(d, nparm)).print(10, 2);
  }
  System.out.print("e0 = " + e0 + ": ");
  System.out.print("moved from ");
        (new Matrix(a, nparm)).print(10, 2);
  System.out.print("e1 = " + e1 + ": ");
  if (e1 < e0) {
    System.out.print("to ");
          (new Matrix(na, nparm)).print(10, 2);
  }
  else {
    System.out.println("move rejected");
  }
      }

View Full Code Here

      System.err.println("Exception caught: " + ex.getMessage());
      System.exit(1); 
    }


    System.out.print("desired solution "); 
    (new Matrix(areal, areal.length)).print(10, 2);


    System.exit(0);
  } //main

View Full Code Here

    s[i] = 1.;
    i++;
  }
      }
      System.out.print("quad x= "); 
      (new Matrix(x)).print(10, 2);


      System.out.print("quad y= "); 
      (new Matrix(y,npts)).print(10, 2);




      o[0] = x;
      o[1] = a;
      o[2] = y;

View Full Code Here

  }
  
  @Override
  public void invert() throws LMAMatrix.InvertException {
    try {
      Matrix m = inverse();
      setMatrix(0, this.getRowDimension() - 1, 0, getColumnDimension() - 1, m);
    }
    catch (RuntimeException e) {
      StringWriter s = new StringWriter();
      PrintWriter p = new PrintWriter(s);

View Full Code Here

      return false;
    }


    final SlicingTransform t1t2 = t1.concatenate( t2 );


    final Matrix mt1 = new Matrix( t1.getMatrix() );
    final Matrix mt2 = new Matrix( t2.getMatrix() );
    final Matrix mt1t2 = new Matrix( t1t2.getMatrix() );


    if ( mt1.times( mt2 ).minus( mt1t2 ).normF() > 0.1 )
    {
      System.out.println( "=======================" );
      System.out.println( "t1: " + t1.numSourceDimensions() + " -> " + t1.numTargetDimensions() + " (n -> m)" );
      System.out.println( "t2: " + t2.numSourceDimensions() + " -> " + t2.numTargetDimensions() + " (n -> m)" );
      System.out.println( "t1t2: " + t1t2.numSourceDimensions() + " -> " + t1t2.numTargetDimensions() + " (n -> m)" );


      System.out.print( "t1 = " );
      mt1.print( 1, 0 );
      System.out.print( "t2 = " );
      mt2.print( 1, 0 );
      System.out.print( "t1t2 = " );
      mt1t2.print( 1, 0 );
      System.out.print( "t1 x t2 = " );
      mt1.times( mt2 ).print( 1, 0 );


      System.out.println( "wrong result" );
      System.out.println( "=======================" );

View Full Code Here

      return false;
    }


    final SlicingTransform t1t2 = t2.preConcatenate( t1 );


    final Matrix mt1 = new Matrix( t1.getMatrix() );
    final Matrix mt2 = new Matrix( t2.getMatrix() );
    final Matrix mt1t2 = new Matrix( t1t2.getMatrix() );


    if ( mt1.times( mt2 ).minus( mt1t2 ).normF() > 0.1 )
    {
      System.out.println( "=======================" );
      System.out.println( "t1: " + t1.numSourceDimensions() + " -> " + t1.numTargetDimensions() + " (n -> m)" );
      System.out.println( "t2: " + t2.numSourceDimensions() + " -> " + t2.numTargetDimensions() + " (n -> m)" );
      System.out.println( "t1t2: " + t1t2.numSourceDimensions() + " -> " + t1t2.numTargetDimensions() + " (n -> m)" );


      System.out.print( "t1 = " );
      mt1.print( 1, 0 );
      System.out.print( "t2 = " );
      mt2.print( 1, 0 );
      System.out.print( "t1t2 = " );
      mt1t2.print( 1, 0 );
      System.out.print( "t1 x t2 = " );
      mt1.times( mt2 ).print( 1, 0 );


      System.out.println( "wrong result" );
      System.out.println( "=======================" );

View Full Code Here

      return false;
    }


    final MixedTransform t1t2 = t1.concatenate( t2 );


    final Matrix mt1 = new Matrix( t1.getMatrix() );
    final Matrix mt2 = new Matrix( t2.getMatrix() );
    final Matrix mt1t2 = new Matrix( t1t2.getMatrix() );


    if ( mt1.times( mt2 ).minus( mt1t2 ).normF() > 0.1 )
    {
      System.out.println( "=======================" );
      System.out.println( "t1: " + t1.numSourceDimensions() + " -> " + t1.numTargetDimensions() + " (n -> m)" );
      System.out.println( "t2: " + t2.numSourceDimensions() + " -> " + t2.numTargetDimensions() + " (n -> m)" );
      System.out.println( "t1t2: " + t1t2.numSourceDimensions() + " -> " + t1t2.numTargetDimensions() + " (n -> m)" );


      System.out.print( "t1 = " );
      mt1.print( 1, 0 );
      System.out.print( "t2 = " );
      mt2.print( 1, 0 );
      System.out.print( "t1t2 = " );
      mt1t2.print( 1, 0 );
      System.out.print( "t1 x t2 = " );
      mt1.times( mt2 ).print( 1, 0 );


      System.out.println( "wrong result" );
      System.out.println( "=======================" );

View Full Code Here

      return false;
    }


    final MixedTransform t1t2 = t2.preConcatenate( t1 );


    final Matrix mt1 = new Matrix( t1.getMatrix() );
    final Matrix mt2 = new Matrix( t2.getMatrix() );
    final Matrix mt1t2 = new Matrix( t1t2.getMatrix() );


    if ( mt1.times( mt2 ).minus( mt1t2 ).normF() > 0.1 )
    {
      System.out.println( "=======================" );
      System.out.println( "t1: " + t1.numSourceDimensions() + " -> " + t1.numTargetDimensions() + " (n -> m)" );
      System.out.println( "t2: " + t2.numSourceDimensions() + " -> " + t2.numTargetDimensions() + " (n -> m)" );
      System.out.println( "t1t2: " + t1t2.numSourceDimensions() + " -> " + t1t2.numTargetDimensions() + " (n -> m)" );


      System.out.print( "t1 = " );
      mt1.print( 1, 0 );
      System.out.print( "t2 = " );
      mt2.print( 1, 0 );
      System.out.print( "t1t2 = " );
      mt1t2.print( 1, 0 );
      System.out.print( "t1 x t2 = " );
      mt1.times( mt2 ).print( 1, 0 );


      System.out.println( "wrong result" );
      System.out.println( "=======================" );

View Full Code Here

  public static void main( final String[] args )
  {
    final MixedTransformConcatenateTest test = new MixedTransformConcatenateTest();
    test.setUp();


    final Matrix m_tr1 = new Matrix( test.tr1.getMatrix() );
    final Matrix m_tr2 = new Matrix( test.tr2.getMatrix() );
    final Matrix m_perm1 = new Matrix( test.perm1.getMatrix() );
    final Matrix m_rot1 = new Matrix( test.rot1.getMatrix() );
    final Matrix m_proj1 = new Matrix( test.proj1.getMatrix() );
    final Matrix m_proj2 = new Matrix( test.proj2.getMatrix() );
    final Matrix m_comp1 = new Matrix( test.comp1.getMatrix() );
    final Matrix m_slicing1 = new Matrix( test.slicing1.getMatrix() );


    System.out.print( "tr1 = " );
    m_tr1.print( 1, 0 );


    System.out.print( "tr2 = " );
    m_tr2.print( 1, 0 );


    System.out.print( "perm1 = " );
    m_perm1.print( 1, 0 );


    System.out.print( "rot1 = " );
    m_rot1.print( 1, 0 );


    System.out.print( "proj1 = " );
    m_proj1.print( 1, 0 );


    System.out.print( "proj2 = " );
    m_proj2.print( 1, 0 );


    System.out.print( "comp1 = " );
    m_comp1.print( 1, 0 );


    System.out.print( "m_slicing1 = " );
    m_slicing1.print( 1, 0 );
  }

View Full Code Here

0 1 2 3 4 5 6 7 8 9

TOP

Related Classes of Jama.Matrix

br.edu.utfpr.cm.JGitMinerWeb.util.matrix.DoubleDenseMatrix

ca.nengo.dynamics.impl.CanonicalModel

ca.nengo.dynamics.impl.ImpulseIntegral

ca.nengo.math.impl.PolynomialCurveFitter

ca.nengo.math.impl.WeightedCostApproximator

ca.nengo.math.impl.WeightedCostApproximatorTest

ca.nengo.model.nef.impl.AdditiveGaussianExpressModel

ca.nengo.model.nef.impl.DecodableEnsembleImpl

caterpillarssa.SpectrumAnalysis

com.heatonresearch.aifh.regression.TrainLeastSquares

All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.