Examples of LSolver

MacOS X - specific instructions

It may be advisable to build gcc first (gmp web page suggests it will not work when built with XCode so using a recent gcc may be advisable), It is worth noting that the full path of the interface is encoded in the shared libraries hence moving this interface is only possible if libraries are rebuilt. For this reason, LSolver is stuck in the rpnicore rather than in learner/linear. Instructions for building OpenBLAS on Linux, Uncompress the .zip, then from the main directory run 
 make PREFIX=/usr/local/soft/OpenBLAS-be853da DYNAMIC_ARCH=1 UTEST_CHECK=1 NO_LAPACK=1	NO_AFFINITY=1
The first parameter is responsible for building in support for all CPU architectures, otherwise only the current one will This will attempt to download CUnit and then fail. 
 cd utest;make
should build and install CUnit into utest. Subsequently doing the following command from the main directory (the one just above utest) 
  ln -s utest/CUnit-2.1-2/include/CUnit .
would fix unit tests that ship with OpenBLAS and they will run when the above make command is run. After doing 
  make PREFIX=/usr/local/soft/OpenBLAS-be853da DYNAMIC_ARCH=1 UTEST_CHECK=1 NO_LAPACK=1 install
OpenBLAS is available. SuiteSparse_config.mk should have 
 CC=x86_64-w64-mingw32-gcc AR=x86_64-w64-mingw32-ar RANLIB = x86_64-w64-mingw32-ranlib F77 = x86_64-w64-mingw32-gfortran BLAS = /usr/local/soft/OpenBLAS-be853da/lib/libopenblas.lib INSTALL_LIB=/usr/local/soft/umfpack-5.6.1-openblas/lib INSTALL_INCLUDE=/usr/local/soft/umfpack-5.6.1-openblas/include
OpenBLAS on Windows7-x86_64 First, build OpenBLAS with both options for Unix (above) and those mentioned in quickbuild.64bit.

When building UMFPACK, most options are the same as for Unix but there are a few changes,

Finally, Makefile in the UMFPACK/Demo directory has the "run" target, where it runs demos and records the results. Every time it runs a demo, the output has to be piped through d2u to correct line endings on Windows and through sed 's/e\([+-]\)0\([0-9]\+\)/e\1\2/g' to convert the differences in which exponents are output on Windows. Hence the extra content is like 
 ... _demo | d2u |  sed  's/e\([+-]\)0\([0-9]\+\)/e\1\2/g' > my_ ... .out

Examples of statechum.analysis.learning.rpnicore.LSolver

156
157
158
159
160
161
162
163
164
      }
    }, IllegalArgumentException.class,"inconsistent dimension");   
    checkForCorrectException(new whatToRun() {
      @SuppressWarnings("unused")
      public @Override void run() throws NumberFormatException {
        new LSolver(Ap, Ai, Ax, b, x);
      }
    }, IllegalArgumentException.class,"inconsistent dimension");   
  }
View Full Code Here

Examples of statechum.analysis.learning.rpnicore.LSolver

179
180
181
182
183
184
185
186
187
      }
    }, IllegalArgumentException.class,"too few");   
    checkForCorrectException(new whatToRun() {
      @SuppressWarnings("unused")
      public @Override void run() throws NumberFormatException {
        new LSolver(Ap, Ai, Ax, b, x);
      }
    }, IllegalArgumentException.class,"too few");   
  }
View Full Code Here

Examples of statechum.analysis.learning.rpnicore.LSolver

201
202
203
204
205
206
207
208
209
      }
    }, IllegalArgumentException.class,"Ap[0] should be 0");   
    checkForCorrectException(new whatToRun() {
      @SuppressWarnings("unused")
      public @Override void run() throws NumberFormatException {
        new LSolver(Ap, Ai, Ax, b, x);
      }
    }, IllegalArgumentException.class,"Ap[0] should be 0");   
  }
View Full Code Here

Examples of statechum.analysis.learning.rpnicore.LSolver

223
224
225
226
227
228
229
230
231
      }
    }, IllegalArgumentException.class,"inconsistent dimension");   
    checkForCorrectException(new whatToRun() {
      @SuppressWarnings("unused")
      public @Override void run() throws NumberFormatException {
        new LSolver(Ap, Ai, Ax, b, x);
      }
    }, IllegalArgumentException.class,"inconsistent dimension");   
  }
View Full Code Here

Examples of statechum.analysis.learning.rpnicore.LSolver

245
246
247
248
249
250
251
252
253
      }
    }, IllegalArgumentException.class,"inconsistent dimension");   
    checkForCorrectException(new whatToRun() {
      @SuppressWarnings("unused")
      public @Override void run() throws NumberFormatException {
        new LSolver(Ap, Ai, Ax, b, x);
      }
    }, IllegalArgumentException.class,"inconsistent dimension");   
  }
View Full Code Here

Examples of statechum.analysis.learning.rpnicore.LSolver

267
268
269
270
271
272
273
274
275
      }
    }, IllegalArgumentException.class,"inconsistent dimension");   
    checkForCorrectException(new whatToRun() {
      @SuppressWarnings("unused")
      public @Override void run() throws NumberFormatException {
        new LSolver(Ap, Ai, Ax, b, x);
      }
    }, IllegalArgumentException.class,"inconsistent dimension");   
  }
View Full Code Here

Examples of statechum.analysis.learning.rpnicore.LSolver

289
290
291
292
293
294
295
296
297
      }
    }, IllegalArgumentException.class,"inconsistent dimension");   
    checkForCorrectException(new whatToRun() {
      @SuppressWarnings("unused")
      public @Override void run() throws NumberFormatException {
        new LSolver(Ap, Ai, Ax, b, x);
      }
    }, IllegalArgumentException.class,"inconsistent dimension");   
  }
View Full Code Here

Examples of statechum.analysis.learning.rpnicore.LSolver

299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
  @Test
  public final void testExternalSolver_fail6()
  {
    final int size = testMatrix.rows();
    final DoubleMatrix2D matrix = DoubleFactory2D.sparse.make(size,size);
    final LSolver solver = new LSolver(matrix,DoubleFactory1D.dense.make(new double[]{8., 45., -3., 3., 19.}));
    final DoubleMatrix1D x = DoubleFactory1D.dense.make(testMatrix.rows());
   
    for(int i=0;i<solver.j_b.length;++i) x.setQuick(i, solver.j_b[i]);
   
    checkForCorrectException(new whatToRun() {
      public @Override void run() throws NumberFormatException {
        solver.solveExternally(1);
      }
    }, IllegalArgumentException.class,"singular");   

    checkForCorrectException(new whatToRun() {
      public @Override void run() throws NumberFormatException {
        LUDecompositionQuick coltSolver = new LUDecompositionQuick();
        coltSolver.decompose(matrix);coltSolver.setLU(matrix);
        coltSolver.solve(x);
      }
    }, IllegalArgumentException.class,"singular");   

    checkForCorrectException(new whatToRun() {
      public @Override void run() throws NumberFormatException {
        for(int i=0;i<testMatrix.rows();++i) solver.j_x[i]=0;
        solver.solveUsingColt();
      }
    }, IllegalArgumentException.class,"singular");   
  }
View Full Code Here

Examples of statechum.analysis.learning.rpnicore.LSolver

341
342
343
344
345
346
347
348
349
350
351
352
353
354
3  0  4 0 6
0 -1 -3 2 0
0  0  1 0 0
0  4  2 0 1
     */
    LSolver s = new LSolver(testMatrix,DoubleFactory1D.dense.make(testMatrix.rows(),0));
    Assert.assertArrayEquals(new int[] {0, 2, 5, 9, 10, 12}, s.j_Ap);
    Assert.assertArrayEquals(new int[] {0,  1,  0,   2,  4,  1,  2,  3,   4,  2,  1,  4}, s.j_Ai);
    for(int i=0;i<s.j_Ap.length;++i)
      Assert.assertEquals(new double[]{2., 3., 3., -1., 4., 4., -3., 1., 2., 2., 6., 1.}[i], s.j_Ax[i],Configuration.fpAccuracy);
   
    DoubleMatrix2D mat = s.toDoubleMatrix2D();
    Assert.assertEquals(testMatrix, mat);
  }
View Full Code Here

Examples of statechum.analysis.learning.rpnicore.LSolver

375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
      matrix.setQuick(x, y, 0.5);
    }
    final DoubleMatrix1D b = DoubleFactory1D.dense.make(size);
    b.assign(randomGenerator);

    final LSolver solver = new LSolver(matrix,b);
    long tmStarted = new Date().getTime();
    LUDecompositionQuick coltSolver = new LUDecompositionQuick();
    coltSolver.decompose(matrix);coltSolver.setLU(matrix);
   
    DoubleMatrix1D vector = b.copy();
    coltSolver.solve(vector);
    long tmFinished = new Date().getTime();
    System.out.println(" time taken: "+((double)tmFinished-tmStarted)/1000);
    verifyAxb(matrix, b, vector);
   
    tmStarted = new Date().getTime();
    solver.solveExternally(1);
    tmFinished = new Date().getTime();
    System.out.println(" time taken: "+((double)tmFinished-tmStarted)/1000);   
    for(int i=0;i<matrix.rows();++i)
      Assert.assertEquals(solver.j_x[i], vector.getQuick(i),Configuration.fpAccuracy);
    verifyAxb(solver);
View Full Code Here
TOP
Copyright © 2018 www.massapi.com. 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.