{"title":"大规模一般形式正则化的混合 LSMR 算法","authors":"Yanfei Yang","doi":"arxiv-2409.09104","DOIUrl":null,"url":null,"abstract":"The hybrid LSMR algorithm is proposed for large-scale general-form\nregularization. It is based on a Krylov subspace projection method where the\nmatrix $A$ is first projected onto a subspace, typically a Krylov subspace,\nwhich is implemented via the Golub-Kahan bidiagonalization process applied to\n$A$, with starting vector $b$. Then a regularization term is employed to the\nprojections. Finally, an iterative algorithm is exploited to solve a least\nsquares problem with constraints. The resulting algorithms are called the\n{hybrid LSMR algorithm}. At every step, we exploit LSQR algorithm to solve the\ninner least squares problem, which is proven to become better conditioned as\nthe number of $k$ increases, so that the LSQR algorithm converges faster. We\nprove how to select the stopping tolerances for LSQR in order to guarantee that\nthe regularized solution obtained by iteratively computing the inner least\nsquares problems and the one obtained by exactly computing the inner least\nsquares problems have the same accuracy. Numerical experiments illustrate that\nthe best regularized solution by the hybrid LSMR algorithm is as accurate as\nthat by JBDQR which is a joint bidiagonalization based algorithm.","PeriodicalId":501162,"journal":{"name":"arXiv - MATH - Numerical Analysis","volume":"101 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2024-09-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Hybrid LSMR algorithms for large-scale general-form regularization\",\"authors\":\"Yanfei Yang\",\"doi\":\"arxiv-2409.09104\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"The hybrid LSMR algorithm is proposed for large-scale general-form\\nregularization. It is based on a Krylov subspace projection method where the\\nmatrix $A$ is first projected onto a subspace, typically a Krylov subspace,\\nwhich is implemented via the Golub-Kahan bidiagonalization process applied to\\n$A$, with starting vector $b$. Then a regularization term is employed to the\\nprojections. Finally, an iterative algorithm is exploited to solve a least\\nsquares problem with constraints. The resulting algorithms are called the\\n{hybrid LSMR algorithm}. At every step, we exploit LSQR algorithm to solve the\\ninner least squares problem, which is proven to become better conditioned as\\nthe number of $k$ increases, so that the LSQR algorithm converges faster. We\\nprove how to select the stopping tolerances for LSQR in order to guarantee that\\nthe regularized solution obtained by iteratively computing the inner least\\nsquares problems and the one obtained by exactly computing the inner least\\nsquares problems have the same accuracy. Numerical experiments illustrate that\\nthe best regularized solution by the hybrid LSMR algorithm is as accurate as\\nthat by JBDQR which is a joint bidiagonalization based algorithm.\",\"PeriodicalId\":501162,\"journal\":{\"name\":\"arXiv - MATH - Numerical Analysis\",\"volume\":\"101 1\",\"pages\":\"\"},\"PeriodicalIF\":0.0000,\"publicationDate\":\"2024-09-13\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"arXiv - MATH - Numerical Analysis\",\"FirstCategoryId\":\"1085\",\"ListUrlMain\":\"https://doi.org/arxiv-2409.09104\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"\",\"JCRName\":\"\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - MATH - Numerical Analysis","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2409.09104","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
Hybrid LSMR algorithms for large-scale general-form regularization
The hybrid LSMR algorithm is proposed for large-scale general-form
regularization. It is based on a Krylov subspace projection method where the
matrix $A$ is first projected onto a subspace, typically a Krylov subspace,
which is implemented via the Golub-Kahan bidiagonalization process applied to
$A$, with starting vector $b$. Then a regularization term is employed to the
projections. Finally, an iterative algorithm is exploited to solve a least
squares problem with constraints. The resulting algorithms are called the
{hybrid LSMR algorithm}. At every step, we exploit LSQR algorithm to solve the
inner least squares problem, which is proven to become better conditioned as
the number of $k$ increases, so that the LSQR algorithm converges faster. We
prove how to select the stopping tolerances for LSQR in order to guarantee that
the regularized solution obtained by iteratively computing the inner least
squares problems and the one obtained by exactly computing the inner least
squares problems have the same accuracy. Numerical experiments illustrate that
the best regularized solution by the hybrid LSMR algorithm is as accurate as
that by JBDQR which is a joint bidiagonalization based algorithm.