Package cern.colt.matrix

Examples of cern.colt.matrix.DoubleMatrix1D


  }
 
  protected DoubleMatrix1D getExpectedMatrix1DSlowly(LearnerGraph gr)
  {
    int size=gr.getStateNumber()*(gr.getStateNumber()+1)/2;
    DoubleMatrix1D result = DoubleFactory1D.dense.make(size);
    DetermineDiagonalAndRightHandSide ddrhInstance = new Linear.DDRH_default();
    for(Entry<CmpVertex,Map<String,CmpVertex>> entryA:gr.transitionMatrix.entrySet())
    {
      // Now iterate through states
      Iterator<Entry<CmpVertex,Map<String,CmpVertex>>> stateB_It = gr.transitionMatrix.entrySet().iterator();
      while(stateB_It.hasNext())
      {
        Entry<CmpVertex,Map<String,CmpVertex>> stateB = stateB_It.next();

        int currentStatePair = gr.wmethod.vertexToIntNR(stateB.getKey(),entryA.getKey());
        ddrhInstance.compute(entryA.getValue(),stateB.getValue());
        result.setQuick(currentStatePair, ddrhInstance.getRightHandSide());
        if (stateB.getKey().equals(entryA.getKey())) break; // we only process a triangular subset.
      }
    }
   
    return result;
View Full Code Here


  {
    LSolver solver = gr.linear.buildMatrix(ThreadNumber);
    DoubleMatrix2D Ax=solver.toDoubleMatrix2D();
    Assert.assertEquals(getExpectedMatrix2DSlowly(gr),Ax);
    if (expectedAx != null) Assert.assertEquals(expectedAx, Ax);
    DoubleMatrix1D b=solver.toDoubleMatrix1D();
    if (expectedB != null) Assert.assertEquals(expectedB, b);Assert.assertEquals(getExpectedMatrix1DSlowly(gr),b);
    solver.solveExternally();// check if we have a solution, just in case it fails.

    // Now check consistency.
    gr.config.setAttenuationK_testOnly(1);DoubleMatrix2D Ax1 = gr.linear.buildMatrix(ThreadNumber).toDoubleMatrix2D();
    gr.config.setAttenuationK(0);DoubleMatrix2D Ax0 = gr.linear.buildMatrix(ThreadNumber).toDoubleMatrix2D();
    DoubleMatrix1D one = DoubleFactory1D.dense.make(Ax1.rows(), 1), a=DoubleFactory1D.dense.make(Ax.rows(), 0);
   
    // check A(1)*one >= 0
    Ax1.zMult(one, a);for(int i=0;i<a.size();++i) Assert.assertTrue(a.getQuick(i)>=0);
   
    // check (A(1)-A(0))*one = b
View Full Code Here

  public final void testBuildMatrix3()
  {
    LearnerGraph gr=new LearnerGraph(buildGraph("A-a->B\nC-a->D""testAddToBuffer2"),config);
    final int size =4*5/2;
    DoubleMatrix2D matrix=DoubleFactory2D.sparse.identity(size);
    DoubleMatrix1D row=DoubleFactory1D.dense.make(size, 0);
    matrix.setQuick(0, 2, -k);row.setQuick(0, 1);
    matrix.setQuick(5, 9, -k);row.setQuick(5, 1);
    matrix.setQuick(3, 7, -k);row.setQuick(3, 1);
   
    checkBuildMatrix(gr,matrix,row);
  }
View Full Code Here

  public final void testBuildMatrix4()
  {
    LearnerGraph gr=new LearnerGraph(buildGraph("A-a->B\nA-b->C\nD-a->C""testAddToBuffer3"),config);
    final int size =4*5/2;
    DoubleMatrix2D matrix=DoubleFactory2D.sparse.identity(size);
    DoubleMatrix1D row=DoubleFactory1D.dense.make(size, 0);
    matrix.setQuick(0, 0, 2);matrix.setQuick(0, 2, -k);matrix.setQuick(0, 5, -k);row.setQuick(0, 2);
    matrix.setQuick(1, 1, 2);
    matrix.setQuick(3, 3, 2);
    matrix.setQuick(6, 6, 2);
 
    matrix.setQuick(9, 5, -k);row.setQuick(9, 1);
    matrix.setQuick(6, 4, -k);row.setQuick(6, 1);
       
    checkBuildMatrix(gr,matrix,row);
  }
View Full Code Here

  public final void testBuildMatrix5()
  {
    LearnerGraph gr=new LearnerGraph(buildGraph("A-a->B\nA-b->C\nD-a->C\nD-b->C","testAddToBuffer4"),config);
    final int size =4*5/2;
    DoubleMatrix2D matrix=DoubleFactory2D.sparse.identity(size);
    DoubleMatrix1D row=DoubleFactory1D.dense.make(size, 0);
    matrix.setQuick(0, 0, 2);matrix.setQuick(0, 2, -k);matrix.setQuick(0, 5, -k);row.setQuick(0, 2);
    matrix.setQuick(1, 1, 2);
    matrix.setQuick(3, 3, 2);
    matrix.setQuick(6, 6, 2);
 
    matrix.setQuick(9, 5, -2*k);row.setQuick(9, 2);
    matrix.setQuick(6, 4, -k);matrix.setQuick(6, 5, -k);row.setQuick(6, 2);
    matrix.setQuick(9, 9, 2);
    matrix.setQuick(8, 8, 2);
    matrix.setQuick(7, 7, 2);
    matrix.setQuick(6, 6, 2);
   
View Full Code Here

  public final void testBuildMatrix6()
  {
    LearnerGraph gr=new LearnerGraph(buildGraph("A-a->B\nA-b->C\nD-a->C\nD-b->C\nD-c->A","testAddToBuffer5"),config);
    final int size =4*5/2;
    DoubleMatrix2D matrix=DoubleFactory2D.sparse.identity(size);
    DoubleMatrix1D row=DoubleFactory1D.dense.make(size, 0);
    matrix.setQuick(0, 0, 2);matrix.setQuick(0, 2, -k);matrix.setQuick(0, 5, -k);row.setQuick(0, 2);
    matrix.setQuick(1, 1, 2);
    matrix.setQuick(3, 3, 2);
    matrix.setQuick(6, 6, 2);
 
    matrix.setQuick(9, 5, -2*k);matrix.setQuick(9, 0, -k);row.setQuick(9, 3);
    matrix.setQuick(6, 4, -k);matrix.setQuick(6, 5, -k);row.setQuick(6, 2);
    matrix.setQuick(9, 9, 3);
    matrix.setQuick(8, 8, 3);
    matrix.setQuick(7, 7, 3);
    matrix.setQuick(6, 6, 3);
   
View Full Code Here

  {
    LearnerGraph gr=new LearnerGraph(buildGraph("A-a->B\nA-b->C\nA-c->C\nD-a->C\nD-b->C\nD-c->A","testAddToBuffer7"),config);

    final int size =4*5/2;
    DoubleMatrix2D matrix=DoubleFactory2D.sparse.identity(size);
    DoubleMatrix1D row=DoubleFactory1D.dense.make(size, 0);
    matrix.setQuick(0, 0, 3);matrix.setQuick(0, 2, -k);matrix.setQuick(0, 5, -2*k);row.setQuick(0, 3);
    matrix.setQuick(1, 1, 3);
    matrix.setQuick(3, 3, 3);
    matrix.setQuick(6, 6, 3);
 
    matrix.setQuick(9, 5, -2*k);matrix.setQuick(9, 0, -k);row.setQuick(9, 3);
    matrix.setQuick(6, 4, -k);matrix.setQuick(6, 5, -k);matrix.setQuick(6, 3, -k);row.setQuick(6, 3);
    matrix.setQuick(9, 9, 3);
    matrix.setQuick(8, 8, 3);
    matrix.setQuick(7, 7, 3);
    matrix.setQuick(6, 6, 3);
   
View Full Code Here

  {
    LearnerGraph gr=new LearnerGraph(buildGraph("A-a->B\nA-b->C\nA-c->C\nD-a->C\nD-b->C\nD-d->C\nD-c->A","testAddToBuffer8"),config);

    final int size =4*5/2;
    DoubleMatrix2D matrix=DoubleFactory2D.sparse.identity(size);
    DoubleMatrix1D row=DoubleFactory1D.dense.make(size, 0);
    matrix.setQuick(0, 0, 3);matrix.setQuick(0, 2, -k);matrix.setQuick(0, 5, -2*k);row.setQuick(0, 3);
    matrix.setQuick(1, 1, 3);
    matrix.setQuick(3, 3, 3);
    matrix.setQuick(6, 6, 3);
 
    matrix.setQuick(9, 5, -3*k);matrix.setQuick(9, 0, -k);row.setQuick(9, 4);
    matrix.setQuick(6, 4, -k);matrix.setQuick(6, 5, -k);matrix.setQuick(6, 3, -k);row.setQuick(6, 3);
    matrix.setQuick(9, 9, 4);
    matrix.setQuick(8, 8, 4);
    matrix.setQuick(7, 7, 4);
    matrix.setQuick(6, 6, 4);
   
View Full Code Here

  @Test
  public final void TestComputeStateCompatibility1a()
  {
    LearnerGraph gr = new LearnerGraph(buildGraph(machineCompatibility1,"TestComputeStateCompatibility1"),config);
    DoubleMatrix1D result = DoubleFactory1D.dense.make(gr.linear.computeStateCompatibility(ThreadNumber,null));
    Assert.assertTrue(DoubleFactory1D.dense.make(new double[]{1+k,1,1,0,0,0}).equals(result));
 
View Full Code Here

  @Test
  public final void TestComputeStateCompatibility2a()
  {
    LearnerGraph gr=new LearnerGraph(buildGraph(machineCompatibility2,"TestComputeStateCompatibility1"),config);
    DoubleMatrix1D result = DoubleFactory1D.dense.make(gr.linear.computeStateCompatibility(ThreadNumber,null));
    DoubleMatrix1D expected=DoubleFactory1D.dense.make(new double[]{1+k,PAIR_INCOMPATIBLE,1,PAIR_INCOMPATIBLE,PAIR_INCOMPATIBLE,0});
    Assert.assertTrue(expected.equals(result));
  }
View Full Code Here

TOP

Related Classes of cern.colt.matrix.DoubleMatrix1D

Copyright © 2018 www.massapicom. All rights reserved.
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.