Yummy ! This old-style code using the stack really sucks, now I am not so sure I will make it :-(

 S.

Le 11/05/2015 15:21, Vincent COUVERT a écrit :
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.




_______________________________________________
dev mailing list
dev@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/dev


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

_______________________________________________
dev mailing list
dev@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/dev

Reply via email to