- Join over
**1.5M+ people** - Join over
**100K+ communities** - Free
**without limits** - Create
**your own community**

- May 14 2019 01:33poulson closed #276
- May 14 2019 01:33poulson commented #276
- May 14 2019 01:30
poulson on master

Update README.md (compare)

- May 14 2019 01:24
poulson on master

added logo (compare)

- May 10 2019 21:37jedbrown commented #276
- May 10 2019 21:25poulson commented #276
- May 10 2019 21:21jedbrown commented #276
- May 10 2019 17:08poulson commented #276
- May 10 2019 17:03tesch1 opened #276
- Apr 10 2019 13:18JM1 commented #275
- Apr 10 2019 13:16JM1 synchronize #275
- Apr 10 2019 08:56JM1 opened #275
- Mar 06 2019 03:47Raviteja1996 closed #274
- Mar 05 2019 05:46Raviteja1996 opened #274
- Feb 11 2019 21:53BenBrock commented #228
- Feb 11 2019 21:52BenBrock commented #228
- Feb 11 2019 21:51poulson commented #228
- Feb 11 2019 21:50poulson commented #228
- Feb 11 2019 21:45BenBrock commented #228
- Jan 23 2019 23:57adambaskerville commented #273

It includes some random projections and could be used as a template for contributing randomized QR

I need to use functions because my matrix vector multiplication actually is a integral operator.

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.

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`

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);

C_copy=C;

C_copy.Resize(C_copy.Height()*C_copy.Width(),1);

The

`CIRC,CIRC`

distribution puts the entire matrix on one process
i started working on it last year but didn't find much material online.

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!

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
The python files were not installed in the right folder as a package under python3.5/site-packages but in the top-level folder

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.

@AidanGG Sorry, I hadn't checked this gitter in a while for temporary reasons. The

`COMPACT_SVD`

is indeed meant to have that behavior: the finite-precision reduction to bidiagonal form introduces a perturbation of magnitude similar to `eps || A ||_2 max(m,n)`

, and so one should ignore singular values below this amount.
any cases of SVD not converging should be a bug: do you have a reproducing example?

there are several SVD algorithms (e.g., Divide and Conquer and QR), so it is unlikely that multiple of them are not converging for the same matrix

and you are correct that ScaLAPACK does not support 64-bit integers (and the problem is in the library itself, not Elemental)

I've been working on getting the LP solvers to work on all of the http://www.netlib.org/lp/data/ examples using a symmetric-indefinite solver and hadn't been checking in as much as I should have

Is there any reason why we can’t further generalize elastic net to penalize l1*||Gx||_1 + l2*||Gx||_2 for a square matrix G.

or even two matrices G,H (where in my case G=H)

Looking at the EN code, it looks like I need to just patch up the corresponding bits of the objective function