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
Jack Poulson
@poulson
I'm in the middle of refactoring the IPMs
a contribution would be welcome
Carlos Borges
@cecborges_twitter
I might implement that for my project. Unfortunately, I can't generate matrices. I implemented a randomized QR that can be changed to become a randomized SVD very easily that takes functions as input. I can send to you in case you want to clean it eventually and include in the code.
Jack Poulson
@poulson
the best mechanism is to submit a Pull Request via GitHub and we can take it from there
I assume you are referring to randomized low-rank QR?
Carlos Borges
@cecborges_twitter
sure. I will talk with the boss here and see if it is ok to send.
I am
Jack Poulson
@poulson
I thought you were in academia? ;-)
Carlos Borges
@cecborges_twitter
I am
Jack Poulson
@poulson
you have to ask for permission to release open source?
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.