Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Activity
  • May 14 01:33
    poulson closed #276
  • May 14 01:33
    poulson commented #276
  • May 14 01:30

    poulson on master

    Update README.md (compare)

  • May 14 01:24

    poulson on master

    added logo (compare)

  • May 10 21:37
    jedbrown commented #276
  • May 10 21:25
    poulson commented #276
  • May 10 21:21
    jedbrown commented #276
  • May 10 17:08
    poulson commented #276
  • May 10 17:03
    tesch1 opened #276
  • Apr 10 13:18
    JM1 commented #275
  • Apr 10 13:16
    JM1 synchronize #275
  • Apr 10 08:56
    JM1 opened #275
  • Mar 06 03:47
    Raviteja1996 closed #274
  • Mar 05 05:46
    Raviteja1996 opened #274
  • Feb 11 21:53
    BenBrock commented #228
  • Feb 11 21:52
    BenBrock commented #228
  • Feb 11 21:51
    poulson commented #228
  • Feb 11 21:50
    poulson commented #228
  • Feb 11 21:45
    BenBrock commented #228
  • Jan 23 23:57
    adambaskerville commented #273
Carlos Borges
@cecborges_twitter
I am at ICES. I worked with George
Jack Poulson
@poulson
fair enough
Carlos Borges
@cecborges_twitter
I don't know if he wants to release it. I think he wants it. But I would not release it before asking him
Jack Poulson
@poulson
ah, fair enough if it is existing project
I would worry more about whether Elemental is the right place to release rather than if it should be released though
Carlos Borges
@cecborges_twitter
It is a existing project. We are implementing FAIMS(method for solving inverse scattering problems) at 3D with several processors. BTW, why are you shipping elemental to use with openblas?
I am writing the entire method with elemental instead of Petsc.
Jack Poulson
@poulson
What is wrong with OpenBLAS?
I assume you are implicitly asking why not BLIS?
Carlos Borges
@cecborges_twitter
yes
Jack Poulson
@poulson
(the answer is that OpenBLAS is in more package managers and BLIS is harder to configure)
Carlos Borges
@cecborges_twitter
Also, I don't know if you had this problem, but if I run elemental using more than one processor with multiple threads it is breaking(with openblas). I can actually send you the error output.
I can run using mip with several processors, but if I don't set openblas_num_threads=1 I get some memory problems. Any ideas of what it can be?
Sorry if I look like an idiot in this one. It is just that I can't get it to work in any way possible.
Jack Poulson
@poulson
What build configuration did you use?
There is no reason to feel ashamed.
Could you file a detailed GitHub issue?
Carlos Borges
@cecborges_twitter
I can. I am running some tests right now but I will do it tomorrow.
Jack Poulson
@poulson
No worries.
Carlos Borges
@cecborges_twitter
Thanks for the help.
Jack Poulson
@poulson
Btw, you may want to look at my randomized column pivoted QR
It includes some random projections and could be used as a template for contributing randomized QR
Carlos Borges
@cecborges_twitter
I will take a look. Is it committed to the website? The one that I implemented takes functions as input.
I need to use functions because my matrix vector multiplication actually is a integral operator.
Carlos Borges
@cecborges_twitter
I have also another question regarding the command resize for a a distributed matrix. I would expect it to work similarly to matlab, I mean if you have a 3 x3 matrix and try to resize to a 9x1 matrix, the command would stack the columns of the matrix. Is this what I should expect? I believe that in some of the col and row aligns it was not working like that.
Jack Poulson
@poulson
the Gemm calls can be replaced with std::function calls similar to the way I implemented FGMRES
ah, that would be reshape
which does not exist in Elemental
at least, I don't think I ever implemented it
err, actually I did implement Reshape
Carlos Borges
@cecborges_twitter
I implemented using bind, instead of function, but the code is very similar
Carlos Borges
@cecborges_twitter
Ah ok. I actually made a copy of the matrix using <T,CIRC,CIRC> and them I used resize on this copy and this gives me the same as reshape.
DistMatrix<T,El::CIRC,El::CIRC> C_copy(C.Height(),C.Width(),g);
C_copy=C;
C_copy.Resize(C_copy.Height()*C_copy.Width(),1);
Jack Poulson
@poulson
I would highly recommend calling Reshape instead
The CIRC,CIRC distribution puts the entire matrix on one process
Carlos Borges
@cecborges_twitter
Thanks I will do that, and it also saves me some trouble in having to change the distribution.
Jack Poulson
@poulson
please feel free to ask any questions if you run into anymore road blocks
Carlos Borges
@cecborges_twitter
thank you very much.
i started working on it last year but didn't find much material online.
Demian Wassermann
@demianw
Hi! I'm trying to install Elemental into an Anaconda environment python 3.5. After succesfully checking the python options building and installing with CMAKE, the Elemental python module seems to not be there. Any ideas on this? Thanks!
Jack Poulson
@poulson
This is good motivation to push out the Anaconda package being discussed at elemental/Elemental#207
What CMake configuration flags did you use?
the usual issue is that you need to make sure your PYTHONPATH includes the Python installation directory
Also DYLD_LIBRARY_PATH needs to be set: http://libelemental.org/documentation/dev/build.html#mac-os-x
Demian Wassermann
@demianw
Thanks for the response. I did set up the PYTHONPATH and I tried on a linux machine and OSX. In both cases the installation of the python package did not work
The python files were not installed in the right folder as a package under python3.5/site-packages but in the top-level folder
Ryan H. Lewis
@rhl-
You can control where most elements of the install go. Look at the redhat/elemental-spec for an example
Aidan Dang
@AidanGG
Hi Jack, thanks for merging my PR. I've got some issues with COMPACT_SVD. In https://github.com/elemental/Elemental/blob/master/src/lapack_like/spectral/SVD.cpp#L75, the default tol is used to cull 'low enough but not necessarily zero' singular values. Is this the intended behaviour? Lines 277, 469, 771 and 1126 of src/lapack_like/spectral/BidiagSVD.cpp seem to contradict this, and do not make use of the APosterioriThreshold function.