[Scilab-Dev] Question about sparse(ij,v,mn) internals
Vincent COUVERT
vincent.couvert at scilab-enterprises.com
Mon May 11 15:21:38 CEST 2015
Hello Stéphane,
The source code (Fortran) is available here:
http://gitweb.scilab.org/?p=scilab.git;a=blob;f=scilab/modules/sparse/sci_gateway/fortran/sci_f_sparse.f;h=6172a4fe8c7dafd2ae7e44135db9c4c4be1c9794;hb=master
Regards.
Le 11/05/2015 15:03, Stéphane Mottelet a écrit :
> Hello Scilab devs/masters,
>
> For performance purposes, i need to implement a function which is
> missing to Scilab, allowing to set at once all non-zero terms of a
> sparse matrix. Some other languages have this feature, e.g. Julia
> allows to write such things.
>
> A.nzval=v
>
> of course this is interesting when the sparsity pattern does not
> change but actual element values do. I noticed that the use of
> "sparse" to reconstruct matrices in this context is the main
> bottleneck of my code... In Scilab the logical name of such a command
> should be "spset" and used like this
>
> spset(A,v)
>
> of course, the user has to be aware that Scilab sparse storage is line
> compressed and thus v should be correctly ordered, but this kind of
> feature is really missing. I think I can (try to) write such a
> function, so my question is : where is the source code of the C
> function which is interfaced as "sparse" macro in Scilab ?
>
> S.
>
>
>
More information about the dev
mailing list