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

[Stéphane Mottelet]

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