[Scilab-Dev] Question about sparse(ij,v,mn) internals

[Stéphane Mottelet]

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.



-- 
Département de Génie Informatique
EA 4297 Transformations Intégrées de la Matière Renouvelable
Université de Technologie de Compiègne -  CS 60319
60203 Compiègne cedex