Ticket #868: test.tol

Set triplet = {Set Include( "triplet.oza" )[1]};

Set uniqueTriple = Classify( triplet, Real( Set a, Set b ) {
    Real cmp1 = Compare(a[1],b[1]);
    If( cmp1, cmp1, Compare( a[2],b[2] ) )
 } );

Real check_r =  284;
Real check_c = 137;

Set EvalSet( triplet, Real( Set cell ) {
    Real If( cell[1] == check_r & cell[2] == check_c, {
        WriteLn( "check cell in triplet = " << cell[ 3 ] );
        0
      }, 0 )
  } );

Real M = SetMax( EvalSet( triplet, Real( Set t ) { t[1] } ) );
Real N = SetMax( EvalSet( triplet, Real( Set t ) { t[2] } ) );

Real check1 = Eq( Card( uniqueTriple ), Card( triplet ) );

Real check2 = If( check1, {
    VMatrix A1 = Convert( Triplet( triplet, M, N ), "Cholmod.R.Sparse" );
    Real cell1 = VMatDat( A1, check_r, check_c );
    WriteLn( "check cell in A1 = " << cell1 );
    VMatrix A2 = Convert( TripletUnique( SetMat(triplet), M, N ), "Cholmod.R.Sparse" );
    Real cell2 = VMatDat( A2, check_r, check_c );
    WriteLn( "check cell in A2 = " << cell2 );
    VMatrix diff = A2 - A1;
    Real VMatMax( Abs( diff ) )
  }, 0 );

WriteLn( "Max(Abs(diff)) = " << check2 );