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
they removed all open source links
Carlos Borges
@cecborges_twitter
Hi Jack, are there any plans to include GMRES or any other iterative solvers in Elemental?
Jack Poulson
@poulson
there is a GMRES in Elemental
it just isn't documented yet I think
it is heavily used within Elemental's Interior Point Methods
see include/El/solve/FGMRES.hpp for the flexible GMRES implementation
and include/El/solve/LGMRES.hpp for the left-preconditioned standard GMRES
with that said, the iterative solvers in the library have only existed to the point that they were needed for the more "direct" techniques like the IPMs
not that I am opposed to the iterative solvers improving
Jack Poulson
@poulson
Hmm, it seems that the Clang bug response time is not so impressive :-( https://llvm.org/bugs/show_bug.cgi?id=31249
Carlos Borges
@cecborges_twitter
Thanks for the response. I didn't see it in any of the available documentation. It is good to know. Best
Jack Poulson
@poulson
@cecborges_twitter The documentation is currently very behind due to many recent changes to the IPMs, sparse-direct solvers, etc.
Carlos Borges
@cecborges_twitter
@poulson Thanks for letting me. I noticed that you have LSQR implemented using matrices as input, do you have plans to implement it using functions as input?
Jack Poulson
@poulson
I don't recall committing LSQR
do you have a link?
Carlos Borges
@cecborges_twitter
isn't your least squares function implementing LSQR? I thought I read that somewhere. Maybe I am wrong.
Jack Poulson
@poulson
LeastSquares uses Flexible GMRES with the preconditioner being an iteratively-refined application of a sparse-direct solution to a nearby quasi-semidefinite problem
there is a note that Michael Saunders tends to use LSQR for this rather than FGMRES
I like GMRES(k) and FGMRES(k)
Carlos Borges
@cecborges_twitter
ah ok. my bad, I saw lsqr there and I thought that the implementation was using it. Sorry about that. But anyway, do you plan to implement LSQR using functions as input?
Jack Poulson
@poulson
it would be nice but isn't in my critical path right now
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.