Re: [ViennaCL-devel] shared_ptr and compressed_matrix

2018-06-15 Thread Karl Rupp
Hi, I have discovered something rather odd. If I run an minimal example (as shown in eigen-with-viennacl.cpp) all runs fine where I pass a compressed_matrix directly to the copy e.g. Eigen::SparseMatrix spAm; ... code to fill spAm; viennacl::matrix A =

Re: [ViennaCL-devel] Failing to find platforms with pocl

2017-08-11 Thread Karl Rupp
Hi Charles, I have another curious situation. I have installed pocl 0.14 on a ubuntu 14.04 system. I can install and run clinfo without any problems. However, when I compile and run my context.cpp file (https://github.com/cdeterman/gpuR/blob/develop/src/context.cpp) and try to run the

Re: [ViennaCL-devel] Efficient unary operation

2016-11-23 Thread Karl Rupp
Hi Charles, > Right now, if I want to take the negative of every element in a matrix I > end up doing the following: > > // previously assigned > viennacl::matrix vcl_A; > > // matrix of zeros to subtract from > viennacl::matrix vcl_Z = > viennacl::zero_matrix(vcl_A.size1(),vcl_A.size2()); > > //

Re: [ViennaCL-devel] vector vs vector_base classes

2016-11-10 Thread Karl Rupp
Hi, > What exactly is the distinction between vector and vector_base classes? vector_base is the base class for dense vectors (viennacl::vector) as well as vector proxy objects (vector_range, vector_slice). Since some of the constructors of vector_base are tricky, the recommendation is to use

Re: [ViennaCL-devel] Fw: Matrix * Vector CL_OUT_OF_RESOURCES error

2016-10-21 Thread Karl Rupp
. There is NO warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. Thanks, Andy *From:* Karl Rupp <r...@iue.tuwien.ac.at> *Sent:* Tuesday, October 11, 2016 10:31 AM *To:* Andrew Palumbo; viennacl

Re: [ViennaCL-devel] Fw: Matrix * Vector CL_OUT_OF_RESOURCES error

2016-10-12 Thread Karl Rupp
Hi Andy, > As I mentioned before, I'm using libviennacl-dev version 1.7.1 installed > from the ubunto repo. > > > > When I run your attached code, I get do get an error: > > > andrew@michael:~/Downloads$ g++ DenseVectorMmul.cpp > -I/usr/include/viennacl/ -lOpenCL -o denseVec > >

Re: [ViennaCL-devel] Fw: Matrix * Vector CL_OUT_OF_RESOURCES error

2016-09-15 Thread Karl Rupp
gt; double* B_values_ptr = B_values; > > // this is currently the constructor that we're using through > scala/javacpp. > viennacl::vector B_vec(B_values_ptr, oclCtx.memory_type(), n, > 0, 1); > > // perform multiplication and pass result to a vector constructor > viennacl::ve

Re: [ViennaCL-devel] Matrix * Vector CL_OUT_OF_RESOURCES error

2016-09-15 Thread Karl Rupp
rform multiplication and pass result to a vector constructor > viennacl::vector C_vec(viennacl::linalg::prod(A_dense_matrix , > B_vec)); > > // print out vec > std::cout << "ViennaCL: " << C_vec << std::endl; > > > // just exit with success for now if

Re: [ViennaCL-devel] Calculating matrix inverse

2016-09-08 Thread Karl Rupp
Hi Charles, lu_factorize factors the matrix A into a lower triangular matrix L (with unit diagonal) and an upper triangular matrix U. The values in A are overwritten with these values. If you want to obtain the inverse, you have to call viennacl::linalg::lu_substitute(vcl_A, vcl_B); where

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-17 Thread Karl Rupp
Hi Charles, > Just adding my opinion here as I have been following this thread. Would > it be possible to have both the .so library and header only options > available or is it a strictly 'this-or-that' scenario? a header-only version should be possible, but exposes the user to all kinds of

Re: [ViennaCL-devel] Interfacing with clMAGMA?

2016-08-17 Thread Karl Rupp
Hi Charles, > There is a fair amount of output, hopefully something here provides a > clue that you can understand. Ok, so let me explain the relevant messages: > ViennaCL: Setting handle kernel argument 0xbac6690 at pos 0 for kernel > assign_cpu (...) > ViennaCL: Setting handle kernel argument

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-17 Thread Karl Rupp
Hi Dmitriy, > We could (and probably should?) add such a convenience header file > at the expense of increased compilation times (and reduced > encapsulation of source code against compiler issues). > > > +1 on single header! :) thanks for the feedback:

Re: [ViennaCL-devel] Interfacing with clMAGMA?

2016-08-12 Thread Karl Rupp
Hi Charles, call .handle()/.handle1()/.handle2() to get the abstract memory buffers, and call .opencl_handle() on them to get the cl_mem handles: A.handle().opencl_handle() Similarly, the command queue is obtained with viennacl::ocl::get_queue().handle().get() Unfortunately it's not

Re: [ViennaCL-devel] Does ViennaCL supports opencl 1.1 EP

2016-08-11 Thread Karl Rupp
Hi Tanguy, > we are working with an arm 7 on an imx6q board and wonder if ViennaCL > supports the opencl 1.1 Embedded Profile. No, the embedded profile is not supported. Since the embedded profile is not part of OpenCL 2.0 and since newer mobile GPUs tend to have full OpenCL support, it is

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-07 Thread Karl Rupp
Hi Andy, the relevant tests for sparse matrices times dense matrices are in tests/spmdm.cpp. In particular, I recreated a test case based on your description and couldn't find any issues: viennacl::compressed_matrix compressed_A; viennacl::matrix B1(std_A.size(),

Re: [ViennaCL-devel] compressed_matrix %*% matrix_Base

2016-08-04 Thread Karl Rupp
> ---- > *From:* Andrew Palumbo <ap@outlook.com> > *Sent:* Wednesday, August 3, 2016 6:44:10 PM > *To:* Karl Rupp; viennacl-devel@lists.sourceforge.net > *Subject:* Re: [ViennaCL-devel] compressed_ma

Re: [ViennaCL-devel] Removing Boost elements

2016-08-03 Thread Karl Rupp
Hi Charles, > I know the intent is to remove boost from viennacl. As such, I was > looking at some files I could contribute to that removal. That said, > before I start blindly making changes, is there a plan for how to > replace elements such as `boost::numeric::ublas::prod`? Well, the

Re: [ViennaCL-devel] boost/ublass requirements

2016-07-26 Thread Karl Rupp
Hi Andy, > I was just looking at looking at svd.hpp, qr.hpp and fspai.hpp, and I > notice that there are boost requirements and includes in these files. > > > Is there any way of running svd, qr or cholesky_solve without > instlalling boost? currently this is tied to uBLAS, but we want to get

Re: [ViennaCL-devel] Copying Values out of a compressed_matrix

2016-07-23 Thread Karl Rupp
Hi, > yes. this seems to be the case. if i force out-of-order CSR into > in-order CSR everything seems to work. Can't see the documentation > explicitly mentioning this if this is the case indeed. > > Karl, can you please confirm only in-order CSRs are supported? Thanks! out-of-order CSR works

Re: [ViennaCL-devel] Copying Values out of a compressed_matrix

2016-07-23 Thread Karl Rupp
Hi, > PS > (4) column indices admit out-of-order placements of elements within each > row. Column indices *have* to be in ascending order for sparse matrix-matrix multiplication. Best regards, Karli > > Thank you. > -Dmitriy > > On Fri, Jul 22, 2016 at 12:56 PM, Dmitriy Lyubimov

Re: [ViennaCL-devel] Copying Values out of a compressed_matrix

2016-07-23 Thread Karl Rupp
gt; Thanks very much for your time, Karl- much appreciated. > > > Andy > > ---- > *From:* Karl Rupp <r...@iue.tuwien.ac.at <mailto:r...@iue.tuwien.ac.at>> > *Sent:* Friday, July 22, 2016

Re: [ViennaCL-devel] Copying Values out of a compressed_matrix

2016-07-22 Thread Karl Rupp
ead of product_size_row * x) > > > I've attached the corrected file. > > > Thanks > > > Andy > > > *From:* Andrew Palumbo <ap@outlook.com> > *Sent:* Thursday, July 21, 2016 11:03

Re: [ViennaCL-devel] Copying Values out of a compressed_matrix

2016-07-20 Thread Karl Rupp
ers similarly with the same type > of results. I know that the use of Javacpp obfuscates what the problem > may be. But I believe the Memorry is properly allocated. > > > > Sorry for the mistake. > > > Thanks, > > > Andy > > > ---

Re: [ViennaCL-devel] Copying Values out of a compressed_matrix

2016-07-20 Thread Karl Rupp
Hi Andy, instead of viennacl::backend::memory_copy(), you want to use viennacl::backend::memory_read(), which directly transfers the data into your buffer(s). If you *know* that your handles are in host memory, you can even grab the values directly via

Re: [ViennaCL-devel] Initializing matrices via direct serialization (row major or CCS)

2016-07-18 Thread Karl Rupp
Hi, > Is DGEMM your performance-critical operation? Are there any other > performance-critical operations? > > > For now we are only looking at (especially sparse) blas3 and > decompositions. Basically, your normal R base functionality for > in-memory sparse algebra. Sparse

Re: [ViennaCL-devel] Initializing matrices via direct serialization (row major or CCS)

2016-07-14 Thread Karl Rupp
tiplies (even though we've got some improvements in a pull request), so for this particular operation you will get better performance from MKL, OpenBLAS, or libflame. Best regards, Karli > On Tue, Jul 12, 2016 at 9:27 AM, Karl Rupp <r...@iue.tuwien.ac.at > <mailto:r...@iue.tuwien.ac.at>> wr

Re: [ViennaCL-devel] Initializing matrices via direct serialization (row major or CCS)

2016-07-12 Thread Karl Rupp
he padded dimensions, use internal_size1() and internal_size2() for the internal number of rows and columns. http://viennacl.sourceforge.net/doc/manual-types.html#manual-types-matrix Best regards, Karli > > On Tue, Jul 12, 2016 at 5:53 AM, Karl Rupp <r...@iue.tuwien.ac.at > <mail

Re: [ViennaCL-devel] Initializing matrices via direct serialization (row major or CCS)

2016-07-12 Thread Karl Rupp
Hi Dmitriy, On 07/12/2016 07:17 AM, Dmitriy Lyubimov wrote: > Hi, > > I am trying to create some elementary wrappers for VCL in javacpp. > > Everything goes fine, except i really would rather not use those "cpu" > types (std::map, > std::vector) and rather initialize matrices directly by feeding

Re: [ViennaCL-devel] Custom OpenCL kernel confusion

2016-06-10 Thread Karl Rupp
eger gemm kernel for square or rectangular matrices. I > really didn't think it would be difficult but I think I have > fallen in a rabbit hole at this point and likely just confusing > myself. > > Regards, > Charles > > &

Re: [ViennaCL-devel] Custom OpenCL kernel confusion

2016-06-10 Thread Karl Rupp
Hi Charles, > Here is the current kernel > with all the different attempts commented out (where MdimPad and PdimPad > or the padded dimensions). where is the third dimension? Are you assuming C to be M-by-M? > If I don't have a size condition check, the > device quickly runs out of resources

Re: [ViennaCL-devel] Custom OpenCL kernel confusion

2016-05-23 Thread Karl Rupp
Hey, > Ah yes, thanks Karl. I remember that now. With that said, are there > recommendations on how kernels should be written to address the padded > columns? I am imagining some if/else or loop limits on indices but > thought I would ask here before I start trying to do that. I am trying >

Re: [ViennaCL-devel] Custom OpenCL kernel confusion

2016-05-23 Thread Karl Rupp
Hi, On 05/23/2016 05:38 PM, Charles Determan wrote: > I am experimenting with the custom OpenCL kernel functionality, > specifically a naive matrix multiplication as an example. > > My OpenCL Kernel: > __kernel void iMatMult(const int Mdim, const int Pdim, > __global const

Re: [ViennaCL-devel] gpuR support

2016-05-16 Thread Karl Rupp
Hi Charles, > In light of the recent question coming here regarding gpuR I was > wondering if I should alter the error function within RViennaCL. Do you > mind the users may come here if a function throws a ViennaCL error? As long as the ViennaCL mailing list don't get swamped by gpuR-related

Re: [ViennaCL-devel] Forking multiple OpenCL contexts

2016-05-11 Thread Karl Rupp
ctx); Best regards, Karli > > On Tue, May 10, 2016 at 10:54 AM, Charles Determan > <cdeterma...@gmail.com <mailto:cdeterma...@gmail.com>> wrote: > > That does it, didn't think to look for an explicit context.hpp > file. Makes sense that the it is with the ma

Re: [ViennaCL-devel] Forking multiple OpenCL contexts

2016-05-09 Thread Karl Rupp
Hi Charles, setContext() is not thread-safe, so if mclapply() is executing in parallel, there will be a race. MPI works across processes, so globals are not shared (and hence setContext() has no problems). If you want to run multithreaded, you should manage the contexts explicitly. That is,

Re: [ViennaCL-devel] Setting context platforms and devices

2016-04-28 Thread Karl Rupp
e very helpful to know how important (vs. just "nice to have") such 'native' GPU support is for you and users of gpuR. Best regards, Karli > On Thu, Apr 28, 2016 at 11:12 AM, Karl Rupp <r...@iue.tuwien.ac.at > <mailto:r...@iue.tuwien.ac.at>> wrote: > > H

Re: [ViennaCL-devel] Handling a matrix from OpenCL

2016-03-30 Thread Karl Rupp
2016-03-29 22:24 GMT+02:00 Karl Rupp <r...@iue.tuwien.ac.at > <mailto:r...@iue.tuwien.ac.at>>: > > Hi Jose, > > > I have been looking for documentation regarding how to handle a matrix > > (sparse) in an OpenCL kernel. My main objective is to can

Re: [ViennaCL-devel] Handling a matrix from OpenCL

2016-03-29 Thread Karl Rupp
Hi Jose, > I have been looking for documentation regarding how to handle a matrix > (sparse) in an OpenCL kernel. My main objective is to can fill a sparse > matrix using an OpenCL kernel. > > Would it be possible?? > > P.S. I can eventually precompute the matrix and indexes arrays sizes. Yes,

Re: [ViennaCL-devel] FFT-2D

2016-03-02 Thread Karl Rupp
ve to ensure that the data is > dyadic (or to be zero-padded to make it dyadic)? You can pass any data size, but you will only see good performance with a power of 2. Best regards, Karli > > ---- &

Re: [ViennaCL-devel] Examples of importing raw buffers

2016-03-01 Thread Karl Rupp
Hi Sumit, sorry, I overlooked this email in my mailers threaded view. > I have a question here: > a.) Suppose I have three raw buffers (let us assume float); I would like > to do something like this: > v = (I- a) / (W-a+EPS), where I, W and a are the three buffers and EPS > is a scalar to

Re: [ViennaCL-devel] Rank/Sort Functionality?

2016-02-24 Thread Karl Rupp
Hi Charles, > I read on a prior post on Karl Rupp's website regarding project ideas > for GSoC 2014 > (https://www.karlrupp.net/2014/02/mentored-project-ideas-for-gsoc-2014/). > One of those was to implement sort() methods. Was this ever > accomplished/implemented? I haven't found anything in

Re: [ViennaCL-devel] Current context index?

2016-02-01 Thread Karl Rupp
hing with the structure of these files that is > preventing this from working? > > Regards, > Charles > > On Fri, Jan 29, 2016 at 2:04 PM, Karl Rupp <r...@iue.tuwien.ac.at > <mailto:r...@iue.tuwien.ac.at>> wrote: > > Hi Charles, > > sorry for the

Re: [ViennaCL-devel] Current context index?

2016-01-29 Thread Karl Rupp
Hi Charles, sorry for the late response. Currently there is no way of getting the internal index. You can, however, edit viennacl/ocl/backend.hpp and add a public member function in viennacl::ocl::backend returning the index. Best regards, Karli On 01/29/2016 04:00 PM, Charles Determan

Re: [ViennaCL-devel] resizing columns to front of matrix?

2016-01-20 Thread Karl Rupp
rhead. Do you happen to have a suggestion for an intuitive interface? I can't think of a good one, so I'd rather leave it as-is with the explicit copying as suggested. Best regards, Karli > On Wed, Jan 20, 2016 at 11:08 AM, Karl Rupp <r...@iue.tuwien.ac.at > <mailto:r...@iue

[ViennaCL-devel] ViennaCL 1.7.1 released!

2016-01-20 Thread Karl Rupp
operations A full list of changes is available here: http://viennacl.sourceforge.net/doc/changelog.html Many thanks to the community for the precious input! Best regards, Karl Rupp -- Site24x7 APM Insight: Get Deep Visibility

Re: [ViennaCL-devel] Switching context device types

2016-01-15 Thread Karl Rupp
Hi Charles, > I know I can set a device type with: > > viennacl::ocl::set_context_device_type(id, viennacl::ocl::gpu_tag()); > > but is there a where for me to subsequently change that to `cpu_tag()` > within the same main() function? Depends on what you want to achieve. You can call

Re: [ViennaCL-devel] tests failed

2015-12-15 Thread Karl Rupp
Hi, > I installed pocl-0.12 in Fedora-23, but tests are still failing. (...) Ok, this looks like a bug on our side, as the generator fails to pick a proper fallback. Is upgrading to ViennaCL 1.7.1 (which I'll roll out in the next few days and which will contain a fix) an option for you? Best

Re: [ViennaCL-devel] tests failed

2015-12-08 Thread Karl Rupp
Hi, these errors look a lot like your GPU driver isn't installed properly. Could you please try to reinstall the GPU driver? Best regards, Karli On 12/08/2015 08:01 PM, Ilya Gradina wrote: > Hi all, I am build ViennaCL for Fedora-Rawhide: >

Re: [ViennaCL-devel] Help in Interfacing another BLAS library?

2015-11-02 Thread Karl Rupp
Hi Oswin, >> FYI: ViennaCL 2.x will provide all functionality via a shared C >> library in order to make cases such as the one you describe much >> easier. This means that there will be stable function names for all >> operations. An optional C++ layer will sit on top, rather than the >> other

Re: [ViennaCL-devel] Help in Interfacing another BLAS library?

2015-10-31 Thread Karl Rupp
Hi Oswin, >> Unfortunately, C++ template metaprogramming is not very friendly to >> exchanging library interfaces at this point, free functions would >> be a lot easier. > I know, I have looked into the interfaces and figured that there might > be some rough edges to get around. My biggest

Re: [ViennaCL-devel] Help in Interfacing another BLAS library?

2015-10-30 Thread Karl Rupp
Hi Oswin, > I am currently investigating whether ViennaCL can be used as a possible > GPU backend for our linear Algebra as part of the Machine Learning > Library Shark. We started by using boost::ublas but then decided to > reimplement it in a better way(e.g. we dropped about 70% of ublas code,

Re: [ViennaCL-devel] HSA backend

2015-10-05 Thread Karl Rupp
Hi Vladimir, > As part of the masters thesis I've developed a HSA (hsafoundation.com > ) backend for viennacl library - > https://github.com/vpa1977/viennacl-dev/ . It is still work in progess, > but i was wondering if it would prove useful enough to be integrated >

Re: [ViennaCL-devel] Non-contiguous matrix subset?

2015-08-18 Thread Karl Rupp
Hi Charles, this is currently not implemented. I thought about gather/scatter routines such that you can extract arbitrary patterns from a matrix, but this is (again) not yet available. Best regards, Karli On 08/18/2015 09:13 PM, Charles Determan wrote: I have read in the documentation that

Re: [ViennaCL-devel] Matrix column or row index?

2015-08-17 Thread Karl Rupp
improvements, please let us know. :-) Best regards, Karli On Mon, Aug 17, 2015 at 8:15 AM, Karl Rupp r...@iue.tuwien.ac.at mailto:r...@iue.tuwien.ac.at wrote: Hi Charles, you can extract the row or column of a matrix using viennacl::row() and viennacl::column(), e.g

Re: [ViennaCL-devel] Matrix column or row index?

2015-08-17 Thread Karl Rupp
Hi Charles, you can extract the row or column of a matrix using viennacl::row() and viennacl::column(), e.g. viennacl::matrixT A(N, N); viennacl::vectorT x(N); x = viennacl::row(A, 2); x = viennacl::column(A, 3); The assignment of a vector to a matrix row or matrix column is more involved,

Re: [ViennaCL-devel] ViennaCL reductions

2015-08-05 Thread Karl Rupp
was the matrix size? Best regards, Karli *From:* Karl Rupp r...@iue.tuwien.ac.at *To:* Sumit Kumar dost_4_e...@yahoo.com *Cc:* viennacl-devel@lists.sourceforge.net viennacl-devel@lists.sourceforge.net *Sent:* Monday, August 3

Re: [ViennaCL-devel] ViennaCL Matrix multiplication

2015-08-05 Thread Karl Rupp
On 08/04/2015 08:21 AM, Sumit Kumar wrote: Hi Karl Based on all your suggestions, I came up with this implementation: // Copy Sparse Eigen matrices to Dense viennacl matrices typedef Eigen::MatrixScalarType, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor RMMatrix;

Re: [ViennaCL-devel] ViennaCL reductions

2015-08-05 Thread Karl Rupp
practical matrix sizes ;-) Either way, we are happy about any contributions. :-) Best regards, Karli *From:* Karl Rupp r...@iue.tuwien.ac.at *To:* Sumit Kumar dost_4_e...@yahoo.com *Cc:* viennacl-devel

Re: [ViennaCL-devel] ViennaCL fast_copy() matrices

2015-08-05 Thread Karl Rupp
grab a pointer to the internal data from an Eigen matrix with .data() and pass it to fast_copy(). Best regards, Karli *From:* Karl Rupp r...@iue.tuwien.ac.at *To:* Sumit Kumar dost_4_e...@yahoo.com *Cc:* viennacl-devel

Re: [ViennaCL-devel] CUDA slower than OpenCL in new R implementation?

2015-08-03 Thread Karl Rupp
On Mon, Aug 3, 2015 at 7:37 AM, Karl Rupp r...@iue.tuwien.ac.at mailto:r...@iue.tuwien.ac.at wrote: Hi Charles, I was benchmarking 4096x4096 matrices (again, with my R bindings). By 'slower' I mean that I am observing OpenCL at this size beating the OpenBLAS CPU

Re: [ViennaCL-devel] ViennaCL reductions

2015-08-02 Thread Karl Rupp
*From:* Karl Rupp r...@iue.tuwien.ac.at *To:* Sumit Kumar dost_4_e...@yahoo.com *Cc:* viennacl-devel@lists.sourceforge.net viennacl-devel@lists.sourceforge.net *Sent:* Friday, July 31, 2015 9:04 PM *Subject:* Re: [ViennaCL-devel] ViennaCL

Re: [ViennaCL-devel] CUDA slower than OpenCL in new R implementation?

2015-08-01 Thread Karl Rupp
Hi Charles, can you please quantify what you mean by 'slower'? How does 'slower' change as you increase the problem size? I would not be surprised if you see no performance gains below matrices of size 500-by-500. With the extra back-and-forth through PCI-Express you may even need matrices of

Re: [ViennaCL-devel] ViennaCL reductions

2015-07-31 Thread Karl Rupp
Hi Sumit, I am aware that Eigen can do it for its matrices and I am also aware that VCL cannot do it natively. My question was this: In your example of interfacing with Eigen, you have shown a VCL dense matrix interfacing with an Eigen dense matrix. Do you have any example of interfacing an

[ViennaCL-devel] ViennaCL 1.7.0 released!

2015-07-31 Thread Karl Rupp
. * Interface to/from the Armadillo library [4] A full list of changes is available here: http://viennacl.sourceforge.net/doc/changelog.html Updates to the benchmark section on the ViennaCL webpage in order to showcast the good performance will follow shortly. Best regards, Karl Rupp [1] https

Re: [ViennaCL-devel] ViennaCL reductions

2015-07-31 Thread Karl Rupp
out-of-the-box in ViennaCL to assign a ViennaCL sparse matrix to a ViennaCL dense matrix (and vice versa). Best regards, Karli *From:* Karl Rupp r...@iue.tuwien.ac.at *To:* Sumit Kumar dost_4_e...@yahoo.com *Cc

Re: [ViennaCL-devel] ViennaCL reductions

2015-07-29 Thread Karl Rupp
Hi Sumit, ViennaCL: Sparse Matrix-Matrix Multiplication http://www.iue.tuwien.ac.at/cse/index.php/gsoc/2014/ideas/167-viennacl-sparse-matrix-matrix-multiplication-2.html Description View on www.iue.tuwien.ac.at

Re: [ViennaCL-devel] ViennaCL reductions

2015-07-27 Thread Karl Rupp
any ViennaCL-objects): viennacl::ocl::setup_context(0, viennacl::ocl::platform().devices()[2]); Best regards, Karli *From:* Karl Rupp r...@iue.tuwien.ac.at *To:* Sumit Kumar dost_4_e...@yahoo.com *Cc:* viennacl-devel

Re: [ViennaCL-devel] ViennaCL reductions

2015-07-27 Thread Karl Rupp
Hi Sumit, Thanks for the update. I will check it out. As for the second one, indeed thats what I ended up doing: // Get some context information for OpenCL compatible GPU devices viennacl::ocl::platform pf; std::vectorviennacl::ocl::device devices =

Re: [ViennaCL-devel] ViennaCL reductions

2015-07-27 Thread Karl Rupp
*From:* Karl Rupp r...@iue.tuwien.ac.at *To:* Sumit Kumar dost_4_e...@yahoo.com *Cc:* viennacl-devel@lists.sourceforge.net viennacl-devel@lists.sourceforge.net *Sent:* Monday, July 27, 2015 7:03 PM *Subject:* Re: [ViennaCL-devel] ViennaCL reductions Hi Sumit

Re: [ViennaCL-devel] Column-wise kernels?

2015-07-27 Thread Karl Rupp
Hi Charles, I am working on writing some additional opencl kernels (potentially to incorporate in to viennacl) which involve column-wise reductions. A simple case would simply be the sum of each column of a matrix. However, I am having an extremely difficult time getting my kernel correct

Re: [ViennaCL-devel] Column-wise kernels?

2015-07-27 Thread Karl Rupp
routines are the performance bottleneck. Glad to hear that 1.7.0 is nearly completed. Does that mean we should expect a formal release soon? Yep. Expect the release on Wednesday. Best regards, Karli On Mon, Jul 27, 2015 at 9:57 AM, Karl Rupp r...@iue.tuwien.ac.at mailto:r

Re: [ViennaCL-devel] Column-wise kernels?

2015-07-27 Thread Karl Rupp
we should expect a formal release soon? Yep. Expect the release on Wednesday. Best regards, Karli On Mon, Jul 27, 2015 at 9:57 AM, Karl Rupp r...@iue.tuwien.ac.at mailto:r...@iue.tuwien.ac.at mailto:r...@iue.tuwien.ac.at mailto:r

Re: [ViennaCL-devel] ViennaCL reductions

2015-07-27 Thread Karl Rupp
, Karli *From:* Karl Rupp r...@iue.tuwien.ac.at *To:* Sumit Kumar dost_4_e...@yahoo.com *Cc:* viennacl-devel@lists.sourceforge.net viennacl-devel@lists.sourceforge.net *Sent:* Monday, July 27, 2015 9:13 PM *Subject:* Re

Re: [ViennaCL-devel] ViennaCL reductions

2015-07-26 Thread Karl Rupp
the matrix sparse. RMMatrix_Float is a typedef for a row-major Eigen Dense matrix; SparseMatrix being its sparse equivalenth Thanks and Regards Sumit *From:* Karl Rupp r...@iue.tuwien.ac.at *To:* Charles Determan cdeterma

Re: [ViennaCL-devel] ViennaCL reductions

2015-07-26 Thread Karl Rupp
Hi Sumit, Thanks for the information. I have a second problem that has cropped up! a.) Apparently, I cannot copy a Row-Major Eigen matrix to a vcl matrix! It comes up with weird errors. I was using the default statements from the examples folder. The moment I changed it to Column Major eigen

Re: [ViennaCL-devel] ViennaCL eigenvectors?

2015-07-14 Thread Karl Rupp
dependency), so it will migrate to master soon. :-) Best regards, Karli On Fri, Jul 10, 2015 at 4:20 AM, Karl Rupp r...@iue.tuwien.ac.at mailto:r...@iue.tuwien.ac.at wrote: Hi Charles, the interface for the power iteration has been updated to also return eigenvectors: https

Re: [ViennaCL-devel] ViennaCL eigenvectors?

2015-07-10 Thread Karl Rupp
Hi Charles, the interface for the power iteration has been updated to also return eigenvectors: https://github.com/viennacl/viennacl-dev/commit/e80cc2141f266eb9b279dd45b7c4075b557bf558 Please let me know if you run into issues. Best regards, Karli On 07/08/2015 06:58 PM, Charles Determan

Re: [ViennaCL-devel] ViennaCL eigenvectors?

2015-07-08 Thread Karl Rupp
Hi Charles, I have seen that I can get all the eigenvalues with the lanczos algorithm in the lanczos.cpp example file but I don't see any documentation on eigenvectors. The only thing I have found is on the feature-improve-lanczos branch

Re: [ViennaCL-devel] NVCC compiler requirement for CUDA backend?

2015-05-26 Thread Karl Rupp
Hi Charles, I am new to the viennacl library after going through some of the documentation there was something I would like a bit of clarification. Could please confirm if the NVCC compiler is indeed required to use the CUDA backend or just for the makefiles of the examples? I ask this as I

Re: [ViennaCL-devel] NVCC compiler requirement for CUDA backend?

2015-05-26 Thread Karl Rupp
rather than C++ template stuff) instead? Best regards, Karli On Tue, May 26, 2015 at 7:47 AM, Karl Rupp r...@iue.tuwien.ac.at mailto:r...@iue.tuwien.ac.at wrote: Hi Charles, I am new to the viennacl library after going through some of the documentation there was something I

Re: [ViennaCL-devel] NVCC compiler requirement for CUDA backend?

2015-05-26 Thread Karl Rupp
the effort for the future. Any input on this is of course appreciated ;-) Best regards, Karli On Tue, May 26, 2015 at 8:14 AM, Karl Rupp r...@iue.tuwien.ac.at mailto:r...@iue.tuwien.ac.at wrote: Hi Charles, On 05/26/2015 03:08 PM, Charles Determan wrote: Thank you Karli

[ViennaCL-devel] IWOCL 2015

2015-03-31 Thread Karl Rupp
Hi, (slightly off-topic, but certainly of interest to users of the OpenCL-backend of ViennaCL) I'd like to take the opportunity to promote the 2015 edition of the Intl. Workshop on OpenCL (IWOCL): http://www.iwocl.org/ The big theme for this year is SyCL, a royalty-free, cross-platform C++

Re: [ViennaCL-devel] SIAM Conference - March 14-15, 2015

2015-03-03 Thread Karl Rupp
Hi Matt, On 03/03/2015 06:26 PM, Matthew Musto wrote: Folks, I see the conference is coming up very soon and I wanted to wish you all good luck on your presentation. Thanks. Since it's a poster presentation, it's more a one-to-one discussion rather than a 20-minute one-man show with

[ViennaCL-devel] ViennaCL 1.6.2 released!

2014-12-11 Thread Karl Rupp
with CUDA on Visual Studio, and further enhances the performance of pipelined iterative solvers. A full list of changes is available here: http://viennacl.sourceforge.net/doc/changelog.html Best regards, Karl Rupp

[ViennaCL-devel] ViennaCLBench 1.0.0 released

2014-12-01 Thread Karl Rupp
for Windows, Linux, and Mac OS. Builds from source are also reasonably convenient and documented, but more involved than just downloading the binaries. With best regards, Karl Rupp -- Download BIRT iHub F-Type - The Free

Re: [ViennaCL-devel] PyViennaCL 1.1.0 release?

2014-11-28 Thread Karl Rupp
Hi Toby, could you please provide a short update on your anticipated PyViennaCL 1.1.0 release? I'd like to know whether I should go ahead with promoting ViennaCL 1.6.1 plus the GUI, or whether it's worth waiting a few (less than 7) more days for PyViennaCL. I'm sorry, I've been quite

[ViennaCL-devel] PyViennaCL 1.1.0 release?

2014-11-26 Thread Karl Rupp
Hi Toby, could you please provide a short update on your anticipated PyViennaCL 1.1.0 release? I'd like to know whether I should go ahead with promoting ViennaCL 1.6.1 plus the GUI, or whether it's worth waiting a few (less than 7) more days for PyViennaCL. Best regards, Karli

Re: [ViennaCL-devel] ViennaCL Benchmark GUI 1.0.0 Release Candidate

2014-11-24 Thread Karl Rupp
Hi Phil, Worked well on my laptop :-) Good to know, thanks! :-) A couple of suggestions: - Maybe use layout N-T for GEMM, or perhaps it is already possible to chose? From my experience NT-col major (TN row major) always leads to higher performance on GEMM. Ah, good point. We can use

Re: [ViennaCL-devel] ViennaCL Benchmark GUI 1.0.0 Release Candidate

2014-11-24 Thread Karl Rupp
Hi Namik, It worked fine on my Windows 7 desktop machine; later I will test it on my Windows 8 laptop and report if anything is missing. I guess linking Qt dynamically is ok for the time being. But we should think about a static build for the long run. A static build is unfortunately not

Re: [ViennaCL-devel] ViennaCL Benchmark GUI 1.0.0 Release Candidate

2014-11-24 Thread Karl Rupp
on overflow. Alright, I'll try with a scrollbar then. Regards, Namik On Mon, Nov 24, 2014 at 11:35 AM, Karl Rupp r...@iue.tuwien.ac.at mailto:r...@iue.tuwien.ac.at wrote: Hi Namik, It worked fine on my Windows 7 desktop machine; later I will test

Re: [ViennaCL-devel] ViennaCL Benchmark GUI 1.0.0 Release Candidate

2014-11-24 Thread Karl Rupp
-Feature in the MatrixMarket :-) Looks like we are good for a nice release tomorrow. Further testing still appreciated ;-) Best regards, Karli On 11/23/2014 09:51 PM, Karl Rupp wrote: Hi guys, a release candidate for the benchmark GUI is available for download, I'd appreciate any testing

[ViennaCL-devel] ViennaCL Benchmark GUI 1.0.0 Release Candidate

2014-11-23 Thread Karl Rupp
Hi guys, a release candidate for the benchmark GUI is available for download, I'd appreciate any testing - particularly the Windows version: ** Windows ** http://viennaclbenchmark.sourceforge.net/ViennaCLBenchmark-1.0.0-RC.zip This is a self-contained package which is ready to launch after

Re: [ViennaCL-devel] Benchmark GUI - 1.0.0 TODOs

2014-11-22 Thread Karl Rupp
packager, we could use Qt Installer Framework http://qt-project.org/doc/qtinstallerframework-1.4/index.html. I've never used before, but considering the good documentation I don't think it'll be a problem. Regards, Namik On Tue, Nov 18, 2014 at 10:46 AM, Karl Rupp r...@iue.tuwien.ac.at

Re: [ViennaCL-devel] Release schedule and project hosting

2014-11-16 Thread Karl Rupp
Hi Namik, 2.) We currently use viennacl.sourceforge.net http://viennacl.sourceforge.net/ as a common domain. However, I wonder whether it better to use three separate domains for better managing the projects? What about using e.g. pyviennacl.sourceforge.net

Re: [ViennaCL-devel] Release schedule and project hosting

2014-11-16 Thread Karl Rupp
Hi Toby, 1.) Toby, how's the status with PyViennaCL? What is missing for a release? No code is missing, but some documentation is, mostly finishing off the new parts added during GSoC. I'm slowly making progress, but I will probably need another 7-10 days for this, as busy as I currently

[ViennaCL-devel] Release schedule and project hosting

2014-11-14 Thread Karl Rupp
Hi, we should push our next wave of releases out next week and have some coordinated promotion. The respective releases are: * PyViennaCL 1.1.0 (am I correct?) * ViennaCL 1.6.1 * ViennaCL Benchmark GUI 1.0.0 This brings up two important questions: 1.) Toby, how's the status with PyViennaCL?

Re: [ViennaCL-devel] Error on ViennaCL build with CUDA

2014-11-12 Thread Karl Rupp
Hi Aanchan, yes, this is a Boost-related problem. Since the error occurs only in the tests, I recommend disabling the tests by disabling ENABLE_TESTING in CMake (e.g. through the command line: cmake .. -DENABLE_TESTING=Off) A student once suggested to replace assert.h in Boost with the version

Re: [ViennaCL-devel] Error on ViennaCL build with CUDA

2014-11-12 Thread Karl Rupp
, Karli On Wed, Nov 12, 2014 at 3:00 PM, Karl Rupp r...@iue.tuwien.ac.at mailto:r...@iue.tuwien.ac.at wrote: Hi Aanchan, yes, this is a Boost-related problem. Since the error occurs only in the tests, I recommend disabling the tests by disabling ENABLE_TESTING in CMake

Re: [ViennaCL-devel] Roadmap update

2014-11-10 Thread Karl Rupp
Hey, I've updated our roadmap taking into account the latest release: https://github.com/viennacl/viennacl-dev/wiki/ViennaCL-Roadmap Feel free to add your topics and post your wishes :-) Awesome! Is it like a christmas present list? Can we post any wish? I'd like a pony,

Re: [ViennaCL-devel] More weird problems

2014-11-08 Thread Karl Rupp
Hi Toby, Thanks! The numerical errors with element-wise operations such as tan() or sin() look okay, that's just numerical noise. The following test cases deserve a closer look, though: test_matrix_matrix_trans_isub_C_float32 test_matrix_matrix_slice_trans_isub_C_float32

  1   2   3   >