A parallel reduced-space sequential-quadratic programming algorithm for frequency-domain small animal optical tomography

Gu, Xuejun; Kim, Hyun Koo; Masciotti, James M.; Hielscher, Andreas H.

doi:10.1117/12.809484

Cited by 4 publications

(3 citation statements)

References 8 publications

Supporting

Mentioning

Contrasting

Unclassified

Order By: Relevance

“…Na biblioteca NLOpt, há algoritmos de dois tipos diferentes dentro dessa classe de algoritmos, sendo os que não se baseiam em condições derivativas, rotulados nos testes e na biblioteca NLOpt como LN, e os que são baseados em métodos dos gradientes para otimização local, rotulados por LD. Há inúmeros algoritmos listados nessas duas categorias e disponibilizados pela NLOpt como os algoritmos de Otimização de Restrições por Aproximações Lineares (Contrained Optimzation by Linear Approximations -COBYLA) [ [38] e três variações dos métodos pré-condicionados truncados de Newton [39], [40]. Como a classe e o tipo de algoritmo são diferentes do ESCH aqui apresentado, optou-se em utilizar o algoritmo NLOPT_LN_COBYLA, por apresentar aproximações lineares que podem ser comparadas a parte dos recursos do DIRECT-L e também por proceder otimizações em que as informações das derivadas do espaço de busca não são [41] informadas com restrições de desigualdades não lineares e igualdades automaticamente modificadas pela combinação de pares de igualdades [40], [41], [42].…”

Section: Resultsunclassified

Evolutionary Strategy Algorithm in a Complex Photonic Coupler Device Optimization

Silva-Sant

Gonçalves

Hernandez

2018

IEEE Latin Am. Trans.

View full text Add to dashboard Cite

show abstract

Section: Resultsunclassified

Evolutionary Strategy Algorithm in a Complex Photonic Coupler Device Optimization

Silva-Sant

Gonçalves

Hernandez

2018

IEEE Latin Am. Trans.

View full text Add to dashboard Cite

show abstract

“…On the other hand, RTE-based DOT codes are based on multiple right hand sides, which correspond to multiple light sources illumination. Traditional methods to solve multiple right hand sides are to solve each right hand side separately or solve multiple right hand sides simultaneously in parallel [24, 25]. However, the extensive computational power those parallel solvers required are not always available; hence we focus here on the numerical method for solving multiple right hand sides simultaneously on a single thread or limited threads (thread number is less than source number).…”

Section: Introductionmentioning

confidence: 99%

Fast linear solver for radiative transport equation with multiple right hand sides in diffuse optical tomography

Jia

Kim

Hielscher

2015

Journal of Quantitative Spectroscopy and Radiative Transfer

Self Cite

View full text Add to dashboard Cite

It is well known that radiative transfer equation (RTE) provides more accurate tomographic results than its diffusion approximation (DA). However, RTE-based tomographic reconstruction codes have limited applicability in practice due to their high computational cost. In this article, we propose a new efficient method for solving the RTE forward problem with multiple light sources in an all-at-once manner instead of solving it for each source separately. To this end, we introduce here a novel linear solver called block biconjugate gradient stabilized method (block BiCGStab) that makes full use of the shared information between different right hand sides to accelerate solution convergence. Two parallelized block BiCGStab methods are proposed for additional acceleration under limited threads situation. We evaluate the performance of this algorithm with numerical simulation studies involving the Delta-Eddington approximation to the scattering phase function. The results show that the single threading block RTE solver proposed here reduces computation time by a factor of 1.5~3 as compared to the traditional sequential solution method and the parallel block solver by a factor of 1.5 as compared to the traditional parallel sequential method. This block linear solver is, moreover, independent of discretization schemes and preconditioners used; thus further acceleration and higher accuracy can be expected when combined with other existing discretization schemes or preconditioners.

show abstract

“…The adaptive finite element method uses local mesh evolution strategy not only to accelerate the simulation with fewer discretized points but also to improve upon the discretization errors, forward solution accuracy and therefore reconstruction quality (Joshi et al 2004, Lv et al 2006). When the simulation region is decomposed into several subdomains, parallel simulation implementation can effectively accelerate the adaptive FEM simulation (Gu et al 2009, Lu and Chatziioannou 2009), enabling fully parallel adaptive finite element photon migration simulation with high-order RTE approximations that are suitable for use within inverse image reconstruction algorithms (Lu et al 2009b). …”

Section: Introductionmentioning

confidence: 99%

A parallel adaptive finite element simplified spherical harmonics approximation solver for frequency domain fluorescence molecular imaging

Zhu

Shen

et al. 2010

Phys. Med. Biol.

View full text Add to dashboard Cite

Fluorescence molecular imaging/tomography may play an important future role in preclinical research and clinical diagnostics. Time-and frequency-domain fluorescence imaging can acquire more measurement information than the continuous wave (CW) counterpart, improving the image quality of fluorescence molecular tomography. Although diffusion approximation (DA) theory has been extensively applied in optical molecular imaging, high-order photon migration models need to be further investigated to match quantitation provided by nuclear imaging. In this paper, a frequencydomain parallel adaptive finite element solver is developed with simplified spherical harmonics (SP N ) approximations. To fully evaluate the performance of the SP N approximations, a fast timeresolved tetrahedron-based Monte Carlo fluorescence simulator suitable for complex heterogeneous geometries is developed using a convolution strategy to realize the simulation of the fluorescence excitation and emission. The validation results show that high-order SP N can effectively correct the modeling errors of the diffusion equation, especially when the tissues have high absorption characteristics or when high modulation frequency measurements are used. Furthermore, the parallel adaptive mesh evolution strategy improves the modeling precision and the simulation speed significantly on a realistic digital mouse phantom. This solver is a promising platform for fluorescence molecular tomography using high-order approximations to the radiative transfer equation.

show abstract

A parallel reduced-space sequential-quadratic programming algorithm for frequency-domain small animal optical tomography

Cited by 4 publications

References 8 publications

Evolutionary Strategy Algorithm in a Complex Photonic Coupler Device Optimization

Evolutionary Strategy Algorithm in a Complex Photonic Coupler Device Optimization

Fast linear solver for radiative transport equation with multiple right hand sides in diffuse optical tomography

A parallel adaptive finite element simplified spherical harmonics approximation solver for frequency domain fluorescence molecular imaging

Contact Info

Product

Resources

About