Increasing lsqcurvefit/nlinfit speed

4 visualizaciones (últimos 30 días)
Raghav Poddar
Raghav Poddar el 6 de Mzo. de 2024
Comentada: Matt J el 6 de Mzo. de 2024
Hi,
I'm looking to increase the speed with which lsqcurvefit is running in my code. The "fitting function" attempts to find a solution by using a lookup table, who's axes are some of my free paramters. I suspect speed is taking a hit because my table is almost 3GB. Right now, I'm sending the table to the anonymous function instead of loading it within the function each iteration (see below - "crossSecSimC2H4" is the database).
[estimCustomCSec(i,:),resnorm,residCsec,~,~,~,Jac] = lsqcurvefit(@(freeParamsC,waveExperiment)fitCustomCrossSecData...
(freeParamsC,waveExperiment,numDens(i),L,absLaser3(:,i,1),bCoeffC2H4,bCoeffCH4,xC2H4range,xCH4range,crossSecSimC2H4...
,crossSecSimCH4,waveExperiment,waveDatabase),freeParamsC,waveExperiment,absLaser3(:,i,1)',[],[],options);
Am I losing speed because even "sending the file" to the function takes long? Or is there no copying of data, in which case the fitting process itself is computationally expensive?

Respuesta aceptada

Matt J
Matt J el 6 de Mzo. de 2024
Editada: Matt J el 6 de Mzo. de 2024
It would definitely be faster to attach the table as an external variable to an anonymous or nested function, than to freshly load it from disk every iteration.
The slow-down might be coming from the time required for table-lookup. A direct computation can sometimes be done entirely within the CPU and its cache, whereas a lookup process, especially in a table that size, would require frequent RAM access.
  2 comentarios
Raghav Poddar
Raghav Poddar el 6 de Mzo. de 2024
I see, that makes sense. Thanks!
Matt J
Matt J el 6 de Mzo. de 2024
You're welcome, but if your question is resolved, please Accept-click the answer.

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre App Building en Help Center y File Exchange.

Etiquetas

Productos


Versión

R2019b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by