Because the random array we build from spartan sometimes will cause "minor
not positive definite error". I guess the reason is that the tiles
generated in each worker is almost the same ...


2014-04-09 15:05 GMT-04:00 Russell Power <notifications@github.com>:

> In tests/test_cholesky.py:
>
> >  #@test_common.with_ctx
> >  #def test_pr(ctx):
> >  def benchmark_cholesky(ctx, timer):
> >    print "#worker:", ctx.num_workers
> >    n = int(math.sqrt(ctx.num_workers))
> > -  ARRAY_SIZE = 2000 * n
> > -  A = expr.randn(ARRAY_SIZE, ARRAY_SIZE, tile_hint=(ARRAY_SIZE/n, ARRAY_SIZE/n))
> > -  A = expr.dot(A, expr.transpose(A)).force()
> > +  ARRAY_SIZE = 1600 * n
> > +
> > +  A = np.random.randn(ARRAY_SIZE, ARRAY_SIZE)
> > +  A = np.dot(A, A.T)
> > +
> > +  util.log_warn('prepare data!')
> > +  A = expr.force(from_numpy(A, tile_hint=(ARRAY_SIZE/n, ARRAY_SIZE/n)))
> > +
> > +  #A = expr.randn(ARRAY_SIZE, ARRAY_SIZE, tile_hint=(ARRAY_SIZE/n, ARRAY_SIZE/n))
> > +  #A = expr.dot(A, expr.transpose(A)).force()
>
> Why the change from Spartan to Numpy+copy to spartan here?
>
> --
> Reply to this email directly or view it on GitHub<https://github.com/spartan-array/spartan/pull/114/files#r11453505>
> .
>



-- 
Chen Qi
Department of Computer Science
School of EECS
Peking University