{"id":58222,"date":"2025-08-16T06:49:50","date_gmt":"2025-08-15T22:49:50","guid":{"rendered":"https:\/\/www.wsisp.com\/helps\/58222.html"},"modified":"2025-08-16T06:49:50","modified_gmt":"2025-08-15T22:49:50","slug":"gcc-c%e5%ae%9e%e7%8e%b0matlab%e7%9f%a9%e9%98%b5%e8%ae%a1%e7%ae%97%e5%92%8c%e6%95%b0%e5%ad%a6%e5%87%bd%e6%95%b0%e5%8a%9f%e8%83%bd","status":"publish","type":"post","link":"https:\/\/www.wsisp.com\/helps\/58222.html","title":{"rendered":"GCC C++\u5b9e\u73b0Matlab\u77e9\u9635\u8ba1\u7b97\u548c\u6570\u5b66\u51fd\u6570\u529f\u80fd"},"content":{"rendered":"

\u5b8c\u6574\u7684GCC C++\u4ee3\u7801\u5b9e\u73b0，\u6574\u5408\u4e86MATIO\u5e93\u3001Eigen\u77e9\u9635\u8ba1\u7b97\u3001\u591a\u7ebf\u7a0b\u5904\u7406\u53ca\u65e5\u5fd7\u7cfb\u7edf，\u5e76\u652f\u6301MATLAB\u98ce\u683c\u7684\u77e9\u9635\u8fd0\u7b97\u548c\u57fa\u672c\u6570\u5b66\u51fd\u6570\u3002<\/p>\n

#<\/span>include<\/span> <iostream><\/span><\/span>
\n#<\/span>include<\/span> <fstream><\/span><\/span>
\n#<\/span>include<\/span> <vector><\/span><\/span>
\n#<\/span>include<\/span> <map><\/span><\/span>
\n#<\/span>include<\/span> <regex><\/span><\/span>
\n#<\/span>include<\/span> <thread><\/span><\/span>
\n#<\/span>include<\/span> <future><\/span><\/span>
\n#<\/span>include<\/span> <queue><\/span><\/span>
\n#<\/span>include<\/span> <mutex><\/span><\/span>
\n#<\/span>include<\/span> <condition_variable><\/span><\/span>
\n#<\/span>include<\/span> <filesystem><\/span><\/span>
\n#<\/span>include<\/span> <sstream><\/span><\/span>
\n#<\/span>include<\/span> <cctype><\/span><\/span>
\n#<\/span>include<\/span> <stdexcept><\/span><\/span>
\n#<\/span>include<\/span> <chrono><\/span><\/span>
\n#<\/span>include<\/span> <iomanip><\/span><\/span>
\n#<\/span>include<\/span> <backward.hpp><\/span><\/span>
\n#<\/span>include<\/span> <matio.h><\/span><\/span>
\n#<\/span>include<\/span> <Eigen\/Dense><\/span><\/span>
\n#<\/span>include<\/span> <cmath><\/span><\/span>
\n#<\/span>include<\/span> <limits><\/span><\/span>
\n#<\/span>include<\/span> <spdlog\/spdlog.h><\/span><\/span>
\n#<\/span>include<\/span> <spdlog\/sinks\/stdout_color_sinks.h><\/span><\/span>
\n#<\/span>include<\/span> <spdlog\/sinks\/basic_file_sink.h><\/span><\/span><\/p>\n

namespace<\/span> fs =<\/span> std::<\/span>filesystem;<\/span>
\nnamespace<\/span> backward {<\/span> backward::<\/span>SignalHandling sh;<\/span> }<\/span><\/p>\n

using<\/span> Eigen::<\/span>MatrixXd;<\/span>
\nusing<\/span> Eigen::<\/span>VectorXd;<\/span>
\nusing<\/span> EigenMap =<\/span> std::<\/span>map<<\/span>std::<\/span>string,<\/span> MatrixXd><\/span>;<\/span><\/p>\n

\/\/ \u5168\u5c40\u65e5\u5fd7\u5668\u521d\u59cb\u5316<\/span>
\nstd::<\/span>shared_ptr<<\/span>spdlog::<\/span>logger><\/span> init_logger<\/span>(<\/span>)<\/span> {<\/span>
\n auto<\/span> console_sink =<\/span> std::<\/span>make_shared<\/span><<\/span>spdlog::<\/span>sinks::<\/span>stdout_color_sink_mt><\/span><\/span><\/span>(<\/span>)<\/span>;<\/span>
\n console_sink-><\/span>set_pattern<\/span>(<\/span>"%^[%Y-%m-%d %H:%M:%S.%e] [%l]%$ %v"<\/span>)<\/span>;<\/span><\/p>\n

auto<\/span> file_sink =<\/span> std::<\/span>make_shared<\/span><<\/span>spdlog::<\/span>sinks::<\/span>basic_file_sink_mt><\/span><\/span><\/span>(<\/span>"matrix_calc.log"<\/span>,<\/span> true<\/span>)<\/span>;<\/span>
\n file_sink-><\/span>set_pattern<\/span>(<\/span>"[%Y-%m-%d %H:%M:%S.%e] [%l] %v"<\/span>)<\/span>;<\/span><\/p>\n

std::<\/span>vector<<\/span>spdlog::<\/span>sink_ptr><\/span> sinks{<\/span>console_sink,<\/span> file_sink}<\/span>;<\/span>
\n auto<\/span> logger =<\/span> std::<\/span>make_shared<\/span><<\/span>spdlog::<\/span>logger><\/span><\/span><\/span>(<\/span>"mat_exec"<\/span>,<\/span> sinks.<\/span>begin<\/span>(<\/span>)<\/span>,<\/span> sinks.<\/span>end<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n logger-><\/span>set_level<\/span>(<\/span>spdlog::<\/span>level::<\/span>debug)<\/span>;<\/span>
\n spdlog::<\/span>register_logger<\/span>(<\/span>logger)<\/span>;<\/span>
\n return<\/span> logger;<\/span>
\n}<\/span><\/p>\n

auto<\/span> logger =<\/span> init_logger<\/span>(<\/span>)<\/span>;<\/span><\/p>\n

\/\/ \u5f02\u5e38\u5904\u7406\u5b8f\u5b9a\u4e49<\/span>
\n#<\/span>define<\/span> THROW_EXCEPTION<\/span>(<\/span>msg)<\/span> <\/span>\\\\<\/span>
\n do<\/span> {<\/span> <\/span>\\\\<\/span>
\n std::<\/span>stringstream ss;<\/span> <\/span>\\\\<\/span>
\n ss <<<\/span> msg <<<\/span> <\/span>" [\u6587\u4ef6: "<\/span> <<<\/span> __FILE__<\/span> <<<\/span> <\/span>" \u884c: "<\/span> <<<\/span> __LINE__<\/span> <<<\/span> <\/span>"]"<\/span>;<\/span> <\/span>\\\\<\/span>
\n logger-><\/span>error<\/span>(<\/span>ss.<\/span>str<\/span>(<\/span>)<\/span>)<\/span>;<\/span> <\/span>\\\\<\/span>
\n throw<\/span> std::<\/span>runtime_error<\/span>(<\/span>ss.<\/span>str<\/span>(<\/span>)<\/span>)<\/span>;<\/span> <\/span>\\\\<\/span>
\n }<\/span> while<\/span>(<\/span>0<\/span>)<\/span><\/span><\/span><\/p>\n

\/\/ \u77e9\u9635\u8fd0\u7b97\u5de5\u5177\u7c7b<\/span>
\nnamespace<\/span> MatrixOps {<\/span>
\n \/\/ \u57fa\u672c\u77e9\u9635\u8fd0\u7b97<\/span>
\n MatrixXd add<\/span>(<\/span>const<\/span> MatrixXd&<\/span> a,<\/span> const<\/span> MatrixXd&<\/span> b)<\/span> {<\/span>
\n if<\/span> (<\/span>a.<\/span>rows<\/span>(<\/span>)<\/span> !=<\/span> b.<\/span>rows<\/span>(<\/span>)<\/span> ||<\/span> a.<\/span>cols<\/span>(<\/span>)<\/span> !=<\/span> b.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"\u77e9\u9635\u52a0\u6cd5\u7ef4\u5ea6\u4e0d\u5339\u914d"<\/span>)<\/span>;<\/span>
\n return<\/span> a +<\/span> b;<\/span>
\n }<\/span><\/p>\n

MatrixXd subtract<\/span>(<\/span>const<\/span> MatrixXd&<\/span> a,<\/span> const<\/span> MatrixXd&<\/span> b)<\/span> {<\/span>
\n if<\/span> (<\/span>a.<\/span>rows<\/span>(<\/span>)<\/span> !=<\/span> b.<\/span>rows<\/span>(<\/span>)<\/span> ||<\/span> a.<\/span>cols<\/span>(<\/span>)<\/span> !=<\/span> b.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"\u77e9\u9635\u51cf\u6cd5\u7ef4\u5ea6\u4e0d\u5339\u914d"<\/span>)<\/span>;<\/span>
\n return<\/span> a –<\/span> b;<\/span>
\n }<\/span><\/p>\n

MatrixXd multiply<\/span>(<\/span>const<\/span> MatrixXd&<\/span> a,<\/span> const<\/span> MatrixXd&<\/span> b)<\/span> {<\/span>
\n if<\/span> (<\/span>a.<\/span>cols<\/span>(<\/span>)<\/span> !=<\/span> b.<\/span>rows<\/span>(<\/span>)<\/span>)<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"\u77e9\u9635\u4e58\u6cd5\u7ef4\u5ea6\u4e0d\u5339\u914d"<\/span>)<\/span>;<\/span>
\n return<\/span> a *<\/span> b;<\/span>
\n }<\/span><\/p>\n

MatrixXd elementwise_multiply<\/span>(<\/span>const<\/span> MatrixXd&<\/span> a,<\/span> const<\/span> MatrixXd&<\/span> b)<\/span> {<\/span>
\n if<\/span> (<\/span>a.<\/span>rows<\/span>(<\/span>)<\/span> !=<\/span> b.<\/span>rows<\/span>(<\/span>)<\/span> ||<\/span> a.<\/span>cols<\/span>(<\/span>)<\/span> !=<\/span> b.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"\u77e9\u9635\u9010\u5143\u7d20\u4e58\u6cd5\u7ef4\u5ea6\u4e0d\u5339\u914d"<\/span>)<\/span>;<\/span>
\n return<\/span> a.<\/span>cwiseProduct<\/span>(<\/span>b)<\/span>;<\/span>
\n }<\/span><\/p>\n

MatrixXd divide<\/span>(<\/span>const<\/span> MatrixXd&<\/span> a,<\/span> const<\/span> MatrixXd&<\/span> b)<\/span> {<\/span>
\n if<\/span> (<\/span>a.<\/span>rows<\/span>(<\/span>)<\/span> !=<\/span> b.<\/span>rows<\/span>(<\/span>)<\/span> ||<\/span> a.<\/span>cols<\/span>(<\/span>)<\/span> !=<\/span> b.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"\u77e9\u9635\u9010\u5143\u7d20\u9664\u6cd5\u7ef4\u5ea6\u4e0d\u5339\u914d"<\/span>)<\/span>;<\/span>
\n return<\/span> a.<\/span>cwiseQuotient<\/span>(<\/span>b)<\/span>;<\/span>
\n }<\/span><\/p>\n

MatrixXd inv<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span>
\n if<\/span> (<\/span>m.<\/span>rows<\/span>(<\/span>)<\/span> !=<\/span> m.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span> THROW_EXCEPTION<\/span>(<\/span>"\u6c42\u9006\u8981\u6c42\u65b9\u9635"<\/span>)<\/span>;<\/span>
\n return<\/span> m.<\/span>inverse<\/span>(<\/span>)<\/span>;<\/span>
\n }<\/span><\/p>\n

MatrixXd pinv<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m,<\/span> double<\/span> tol =<\/span> 1e-6<\/span>)<\/span> {<\/span>
\n Eigen::<\/span>JacobiSVD<<\/span>MatrixXd><\/span> svd<\/span>(<\/span>m,<\/span> Eigen::<\/span>ComputeThinU |<\/span> Eigen::<\/span>ComputeThinV)<\/span>;<\/span>
\n VectorXd sv =<\/span> svd.<\/span>singularValues<\/span>(<\/span>)<\/span>;<\/span>
\n MatrixXd S_inv =<\/span> MatrixXd<\/span>::<\/span>Zero<\/span>(<\/span>m.<\/span>cols<\/span>(<\/span>)<\/span>,<\/span> m.<\/span>rows<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n for<\/span> (<\/span>int<\/span> i =<\/span> 0<\/span>;<\/span> i <<\/span> sv.<\/span>size<\/span>(<\/span>)<\/span>;<\/span> ++<\/span>i)<\/span> {<\/span>
\n if<\/span> (<\/span>sv<\/span>(<\/span>i)<\/span> ><\/span> tol)<\/span> S_inv<\/span>(<\/span>i,<\/span> i)<\/span> =<\/span> 1.0<\/span> \/<\/span> sv<\/span>(<\/span>i)<\/span>;<\/span>
\n }<\/span>
\n return<\/span> svd.<\/span>matrixV<\/span>(<\/span>)<\/span> *<\/span> S_inv *<\/span> svd.<\/span>matrixU<\/span>(<\/span>)<\/span>.<\/span>transpose<\/span>(<\/span>)<\/span>;<\/span>
\n }<\/span><\/p>\n

MatrixXd eig<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span>
\n if<\/span> (<\/span>m.<\/span>rows<\/span>(<\/span>)<\/span> !=<\/span> m.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span> THROW_EXCEPTION<\/span>(<\/span>"\u7279\u5f81\u503c\u8981\u6c42\u65b9\u9635"<\/span>)<\/span>;<\/span>
\n Eigen::<\/span>EigenSolver<<\/span>MatrixXd><\/span> solver<\/span>(<\/span>m)<\/span>;<\/span>
\n return<\/span> solver.<\/span>eigenvalues<\/span>(<\/span>)<\/span>.<\/span>real<\/span>(<\/span>)<\/span>;<\/span>
\n }<\/span><\/p>\n

MatrixXd svd<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span>
\n Eigen::<\/span>JacobiSVD<<\/span>MatrixXd><\/span> svd<\/span>(<\/span>m,<\/span> Eigen::<\/span>ComputeThinU |<\/span> Eigen::<\/span>ComputeThinV)<\/span>;<\/span>
\n return<\/span> svd.<\/span>singularValues<\/span>(<\/span>)<\/span>;<\/span>
\n }<\/span><\/p>\n

MatrixXd chol<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span>
\n if<\/span> (<\/span>m.<\/span>rows<\/span>(<\/span>)<\/span> !=<\/span> m.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span> THROW_EXCEPTION<\/span>(<\/span>"Cholesky\u5206\u89e3\u8981\u6c42\u65b9\u9635"<\/span>)<\/span>;<\/span>
\n Eigen::<\/span>LLT<<\/span>MatrixXd><\/span> llt<\/span>(<\/span>m)<\/span>;<\/span>
\n if<\/span> (<\/span>llt.<\/span>info<\/span>(<\/span>)<\/span> ==<\/span> Eigen::<\/span>NumericalIssue)<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"\u77e9\u9635\u975e\u6b63\u5b9a，\u65e0\u6cd5Cholesky\u5206\u89e3"<\/span>)<\/span>;<\/span>
\n return<\/span> llt.<\/span>matrixL<\/span>(<\/span>)<\/span>;<\/span>
\n }<\/span><\/p>\n

\/\/ \u57fa\u672c\u6570\u5b66\u51fd\u6570<\/span>
\n MatrixXd sin<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>sin<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd cos<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>cos<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd tan<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>tan<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd asin<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>asin<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd acos<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>acos<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd atan<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>atan<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd sinh<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>sinh<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd cosh<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>cosh<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd tanh<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>tanh<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd asinh<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>asinh<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd acosh<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>acosh<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd atanh<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>atanh<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd exp<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>exp<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd log<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>log<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd log10<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>log10<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd log2<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span>
\n return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>log<\/span>(<\/span>)<\/span> \/<\/span> std::<\/span>log<\/span>(<\/span>2<\/span>)<\/span>;<\/span>
\n }<\/span>
\n MatrixXd sqrt<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>sqrt<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd abs<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>abs<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd sign<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span>
\n return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>unaryExpr<\/span>(<\/span>[<\/span>]<\/span>(<\/span>double<\/span> x)<\/span> {<\/span>
\n return<\/span> (<\/span>x ><\/span> 0<\/span>)<\/span> ?<\/span> 1<\/span> :<\/span> (<\/span>(<\/span>x <<\/span> 0<\/span>)<\/span> ?<\/span> –<\/span>1<\/span> :<\/span> 0<\/span>)<\/span>;<\/span>
\n }<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span>
\n }<\/span>
\n MatrixXd round<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>round<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd floor<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>floor<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd ceil<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>ceil<\/span>(<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span> }<\/span>
\n MatrixXd fix<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span>
\n return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>unaryExpr<\/span>(<\/span>[<\/span>]<\/span>(<\/span>double<\/span> x)<\/span> {<\/span>
\n return<\/span> (<\/span>x ><\/span> 0<\/span>)<\/span> ?<\/span> std::<\/span>floor<\/span>(<\/span>x)<\/span> :<\/span> std::<\/span>ceil<\/span>(<\/span>x)<\/span>;<\/span>
\n }<\/span>)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span>
\n }<\/span>
\n MatrixXd pow<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m,<\/span> double<\/span> exponent)<\/span> {<\/span>
\n return<\/span> m.<\/span>array<\/span>(<\/span>)<\/span>.<\/span>pow<\/span>(<\/span>exponent)<\/span>.<\/span>matrix<\/span>(<\/span>)<\/span>;<\/span>
\n }<\/span>
\n MatrixXd atan2<\/span>(<\/span>const<\/span> MatrixXd&<\/span> y,<\/span> const<\/span> MatrixXd&<\/span> x)<\/span> {<\/span>
\n if<\/span> (<\/span>y.<\/span>rows<\/span>(<\/span>)<\/span> !=<\/span> x.<\/span>rows<\/span>(<\/span>)<\/span> ||<\/span> y.<\/span>cols<\/span>(<\/span>)<\/span> !=<\/span> x.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span> {<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"atan2: \u77e9\u9635\u7ef4\u5ea6\u4e0d\u5339\u914d"<\/span>)<\/span>;<\/span>
\n }<\/span>
\n MatrixXd result<\/span>(<\/span>y.<\/span>rows<\/span>(<\/span>)<\/span>,<\/span> y.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n for<\/span> (<\/span>int<\/span> i =<\/span> 0<\/span>;<\/span> i <<\/span> y.<\/span>size<\/span>(<\/span>)<\/span>;<\/span> ++<\/span>i)<\/span> {<\/span>
\n result<\/span>(<\/span>i)<\/span> =<\/span> std::<\/span>atan2<\/span>(<\/span>y<\/span>(<\/span>i)<\/span>,<\/span> x<\/span>(<\/span>i)<\/span>)<\/span>;<\/span>
\n }<\/span>
\n return<\/span> result;<\/span>
\n }<\/span>
\n MatrixXd mod<\/span>(<\/span>const<\/span> MatrixXd&<\/span> x,<\/span> const<\/span> MatrixXd&<\/span> y)<\/span> {<\/span>
\n if<\/span> (<\/span>x.<\/span>rows<\/span>(<\/span>)<\/span> !=<\/span> y.<\/span>rows<\/span>(<\/span>)<\/span> ||<\/span> x.<\/span>cols<\/span>(<\/span>)<\/span> !=<\/span> y.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span> {<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"mod: \u77e9\u9635\u7ef4\u5ea6\u4e0d\u5339\u914d"<\/span>)<\/span>;<\/span>
\n }<\/span>
\n MatrixXd result<\/span>(<\/span>x.<\/span>rows<\/span>(<\/span>)<\/span>,<\/span> x.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n for<\/span> (<\/span>int<\/span> i =<\/span> 0<\/span>;<\/span> i <<\/span> x.<\/span>size<\/span>(<\/span>)<\/span>;<\/span> ++<\/span>i)<\/span> {<\/span>
\n result<\/span>(<\/span>i)<\/span> =<\/span> std::<\/span>fmod<\/span>(<\/span>x<\/span>(<\/span>i)<\/span>,<\/span> y<\/span>(<\/span>i)<\/span>)<\/span>;<\/span>
\n }<\/span>
\n return<\/span> result;<\/span>
\n }<\/span>
\n MatrixXd rem<\/span>(<\/span>const<\/span> MatrixXd&<\/span> x,<\/span> const<\/span> MatrixXd&<\/span> y)<\/span> {<\/span>
\n if<\/span> (<\/span>x.<\/span>rows<\/span>(<\/span>)<\/span> !=<\/span> y.<\/span>rows<\/span>(<\/span>)<\/span> ||<\/span> x.<\/span>cols<\/span>(<\/span>)<\/span> !=<\/span> y.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span> {<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"rem: \u77e9\u9635\u7ef4\u5ea6\u4e0d\u5339\u914d"<\/span>)<\/span>;<\/span>
\n }<\/span>
\n MatrixXd result<\/span>(<\/span>x.<\/span>rows<\/span>(<\/span>)<\/span>,<\/span> x.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n for<\/span> (<\/span>int<\/span> i =<\/span> 0<\/span>;<\/span> i <<\/span> x.<\/span>size<\/span>(<\/span>)<\/span>;<\/span> ++<\/span>i)<\/span> {<\/span>
\n double<\/span> a =<\/span> x<\/span>(<\/span>i)<\/span>;<\/span>
\n double<\/span> b =<\/span> y<\/span>(<\/span>i)<\/span>;<\/span>
\n result<\/span>(<\/span>i)<\/span> =<\/span> (<\/span>b ==<\/span> 0<\/span>)<\/span> ?<\/span> std::<\/span>numeric_limits<\/span><<\/span>double<\/span>><\/span>::<\/span>quiet_NaN<\/span>(<\/span>)<\/span> :<\/span> a –<\/span> b *<\/span> std::<\/span>floor<\/span>(<\/span>a \/<\/span> b)<\/span>;<\/span>
\n }<\/span>
\n return<\/span> result;<\/span>
\n }<\/span>
\n MatrixXd max<\/span>(<\/span>const<\/span> MatrixXd&<\/span> a,<\/span> const<\/span> MatrixXd&<\/span> b)<\/span> {<\/span>
\n if<\/span> (<\/span>a.<\/span>rows<\/span>(<\/span>)<\/span> !=<\/span> b.<\/span>rows<\/span>(<\/span>)<\/span> ||<\/span> a.<\/span>cols<\/span>(<\/span>)<\/span> !=<\/span> b.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span> {<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"max: \u77e9\u9635\u7ef4\u5ea6\u4e0d\u5339\u914d"<\/span>)<\/span>;<\/span>
\n }<\/span>
\n MatrixXd result<\/span>(<\/span>a.<\/span>rows<\/span>(<\/span>)<\/span>,<\/span> a.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n for<\/span> (<\/span>int<\/span> i =<\/span> 0<\/span>;<\/span> i <<\/span> a.<\/span>size<\/span>(<\/span>)<\/span>;<\/span> ++<\/span>i)<\/span> {<\/span>
\n result<\/span>(<\/span>i)<\/span> =<\/span> std::<\/span>max<\/span>(<\/span>a<\/span>(<\/span>i)<\/span>,<\/span> b<\/span>(<\/span>i)<\/span>)<\/span>;<\/span>
\n }<\/span>
\n return<\/span> result;<\/span>
\n }<\/span>
\n MatrixXd min<\/span>(<\/span>const<\/span> MatrixXd&<\/span> a,<\/span> const<\/span> MatrixXd&<\/span> b)<\/span> {<\/span>
\n if<\/span> (<\/span>a.<\/span>rows<\/span>(<\/span>)<\/span> !=<\/span> b.<\/span>rows<\/span>(<\/span>)<\/span> ||<\/span> a.<\/span>cols<\/span>(<\/span>)<\/span> !=<\/span> b.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span> {<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"min: \u77e9\u9635\u7ef4\u5ea6\u4e0d\u5339\u914d"<\/span>)<\/span>;<\/span>
\n }<\/span>
\n MatrixXd result<\/span>(<\/span>a.<\/span>rows<\/span>(<\/span>)<\/span>,<\/span> a.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n for<\/span> (<\/span>int<\/span> i =<\/span> 0<\/span>;<\/span> i <<\/span> a.<\/span>size<\/span>(<\/span>)<\/span>;<\/span> ++<\/span>i)<\/span> {<\/span>
\n result<\/span>(<\/span>i)<\/span> =<\/span> std::<\/span>min<\/span>(<\/span>a<\/span>(<\/span>i)<\/span>,<\/span> b<\/span>(<\/span>i)<\/span>)<\/span>;<\/span>
\n }<\/span>
\n return<\/span> result;<\/span>
\n }<\/span>
\n MatrixXd sum<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span>
\n return<\/span> MatrixXd<\/span>::<\/span>Constant<\/span>(<\/span>1<\/span>,<\/span> 1<\/span>,<\/span> m.<\/span>sum<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n }<\/span>
\n MatrixXd mean<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span>
\n return<\/span> MatrixXd<\/span>::<\/span>Constant<\/span>(<\/span>1<\/span>,<\/span> 1<\/span>,<\/span> m.<\/span>mean<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n }<\/span>
\n MatrixXd prod<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span>
\n double<\/span> product =<\/span> 1.0<\/span>;<\/span>
\n for<\/span> (<\/span>int<\/span> i =<\/span> 0<\/span>;<\/span> i <<\/span> m.<\/span>size<\/span>(<\/span>)<\/span>;<\/span> ++<\/span>i)<\/span> {<\/span>
\n product *=<\/span> m<\/span>(<\/span>i)<\/span>;<\/span>
\n }<\/span>
\n return<\/span> MatrixXd<\/span>::<\/span>Constant<\/span>(<\/span>1<\/span>,<\/span> 1<\/span>,<\/span> product)<\/span>;<\/span>
\n }<\/span>
\n MatrixXd max_val<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span>
\n return<\/span> MatrixXd<\/span>::<\/span>Constant<\/span>(<\/span>1<\/span>,<\/span> 1<\/span>,<\/span> m.<\/span>maxCoeff<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n }<\/span>
\n MatrixXd min_val<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span>
\n return<\/span> MatrixXd<\/span>::<\/span>Constant<\/span>(<\/span>1<\/span>,<\/span> 1<\/span>,<\/span> m.<\/span>minCoeff<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n }<\/span>
\n}<\/span><\/p>\n

\/\/ MAT\u6587\u4ef6\u5904\u7406\u7c7b<\/span>
\nclass<\/span> MatFileHandler<\/span> {<\/span>
\npublic<\/span>:<\/span>
\n static<\/span> EigenMap load<\/span>(<\/span>const<\/span> std::<\/span>string&<\/span> path)<\/span> {<\/span>
\n if<\/span> (<\/span>!<\/span>fs::<\/span>exists<\/span>(<\/span>path)<\/span>)<\/span> THROW_EXCEPTION<\/span>(<\/span>"\u6587\u4ef6\u4e0d\u5b58\u5728: "<\/span> +<\/span> path)<\/span>;<\/span><\/p>\n

mat_t*<\/span> mat =<\/span> Mat_Open<\/span>(<\/span>path.<\/span>c_str<\/span>(<\/span>)<\/span>,<\/span> MAT_ACC_RDONLY)<\/span>;<\/span>
\n if<\/span> (<\/span>!<\/span>mat)<\/span> THROW_EXCEPTION<\/span>(<\/span>"\u65e0\u6cd5\u6253\u5f00MAT\u6587\u4ef6: "<\/span> +<\/span> path)<\/span>;<\/span><\/p>\n

EigenMap variables;<\/span>
\n matvar_t*<\/span> var;<\/span>
\n while<\/span> (<\/span>(<\/span>var =<\/span> Mat_VarReadNext<\/span>(<\/span>mat)<\/span>)<\/span> !=<\/span> nullptr<\/span>)<\/span> {<\/span>
\n if<\/span> (<\/span>var-><\/span>data_type ==<\/span> MAT_T_DOUBLE &&<\/span> var-><\/span>rank ==<\/span> 2<\/span>)<\/span> {<\/span>
\n MatrixXd matrix<\/span>(<\/span>var-><\/span>dims[<\/span>0<\/span>]<\/span>,<\/span> var-><\/span>dims[<\/span>1<\/span>]<\/span>)<\/span>;<\/span>
\n memcpy<\/span>(<\/span>matrix.<\/span>data<\/span>(<\/span>)<\/span>,<\/span> var-><\/span>data,<\/span> var-><\/span>nbytes)<\/span>;<\/span>
\n variables[<\/span>var-><\/span>name]<\/span> =<\/span> matrix;<\/span>
\n logger-><\/span>debug<\/span>(<\/span>"\u52a0\u8f7d\u53d8\u91cf: {} [{}x{}]"<\/span>,<\/span> var-><\/span>name,<\/span>
\n var-><\/span>dims[<\/span>0<\/span>]<\/span>,<\/span> var-><\/span>dims[<\/span>1<\/span>]<\/span>)<\/span>;<\/span>
\n }<\/span>
\n Mat_VarFree<\/span>(<\/span>var)<\/span>;<\/span>
\n }<\/span>
\n Mat_Close<\/span>(<\/span>mat)<\/span>;<\/span>
\n return<\/span> variables;<\/span>
\n }<\/span><\/p>\n

static<\/span> void<\/span> save<\/span>(<\/span>const<\/span> std::<\/span>string&<\/span> path,<\/span> const<\/span> EigenMap&<\/span> variables)<\/span> {<\/span>
\n mat_t*<\/span> mat =<\/span> Mat_CreateVer<\/span>(<\/span>path.<\/span>c_str<\/span>(<\/span>)<\/span>,<\/span> nullptr<\/span>,<\/span> MAT_FT_MAT5)<\/span>;<\/span>
\n if<\/span> (<\/span>!<\/span>mat)<\/span> THROW_EXCEPTION<\/span>(<\/span>"\u65e0\u6cd5\u521b\u5efaMAT\u6587\u4ef6: "<\/span> +<\/span> path)<\/span>;<\/span><\/p>\n

for<\/span> (<\/span>const<\/span> auto<\/span>&<\/span> [<\/span>name,<\/span> matrix]<\/span> :<\/span> variables)<\/span> {<\/span>
\n size_t dims[<\/span>2<\/span>]<\/span> =<\/span> {<\/span>
\n static_cast<\/span><<\/span>size_t><\/span><\/span><\/span>(<\/span>matrix.<\/span>rows<\/span>(<\/span>)<\/span>)<\/span>,<\/span>
\n static_cast<\/span><<\/span>size_t><\/span><\/span><\/span>(<\/span>matrix.<\/span>cols<\/span>(<\/span>)<\/span>)<\/span>
\n }<\/span>;<\/span>
\n matvar_t*<\/span> var =<\/span> Mat_VarCreate<\/span>(<\/span>
\n name.<\/span>c_str<\/span>(<\/span>)<\/span>,<\/span> MAT_C_DOUBLE,<\/span> MAT_T_DOUBLE,<\/span> 2<\/span>,<\/span> dims,<\/span>
\n (<\/span>void<\/span>*<\/span>)<\/span>matrix.<\/span>data<\/span>(<\/span>)<\/span>,<\/span> MAT_F_DONT_COPY_DATA
\n )<\/span>;<\/span>
\n if<\/span> (<\/span>!<\/span>var ||<\/span> Mat_VarWrite<\/span>(<\/span>mat,<\/span> var,<\/span> MAT_COMPRESSION_NONE)<\/span> !=<\/span> 0<\/span>)<\/span> {<\/span>
\n Mat_VarFree<\/span>(<\/span>var)<\/span>;<\/span>
\n Mat_Close<\/span>(<\/span>mat)<\/span>;<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"\u5199\u5165\u53d8\u91cf\u5931\u8d25: "<\/span> +<\/span> name)<\/span>;<\/span>
\n }<\/span>
\n Mat_VarFree<\/span>(<\/span>var)<\/span>;<\/span>
\n logger-><\/span>debug<\/span>(<\/span>"\u4fdd\u5b58\u53d8\u91cf: {} [{}x{}]"<\/span>,<\/span> name,<\/span> dims[<\/span>0<\/span>]<\/span>,<\/span> dims[<\/span>1<\/span>]<\/span>)<\/span>;<\/span>
\n }<\/span>
\n Mat_Close<\/span>(<\/span>mat)<\/span>;<\/span>
\n }<\/span>
\n}<\/span>;<\/span><\/p>\n

\/\/ \u5176\u4ed6\u683c\u5f0f\u8f93\u51fa\u7c7b<\/span>
\nclass<\/span> FormatExporter<\/span> {<\/span>
\npublic<\/span>:<\/span>
\n static<\/span> void<\/span> save_csv<\/span>(<\/span>const<\/span> std::<\/span>string&<\/span> path,<\/span> const<\/span> EigenMap&<\/span> variables)<\/span> {<\/span>
\n std::<\/span>ofstream file<\/span>(<\/span>path)<\/span>;<\/span>
\n if<\/span> (<\/span>!<\/span>file.<\/span>is_open<\/span>(<\/span>)<\/span>)<\/span> THROW_EXCEPTION<\/span>(<\/span>"\u65e0\u6cd5\u521b\u5efaCSV\u6587\u4ef6: "<\/span> +<\/span> path)<\/span>;<\/span><\/p>\n

for<\/span> (<\/span>const<\/span> auto<\/span>&<\/span> [<\/span>name,<\/span> matrix]<\/span> :<\/span> variables)<\/span> {<\/span>
\n file <<<\/span> "=== "<\/span> <<<\/span> name <<<\/span> " ("<\/span> <<<\/span> matrix.<\/span>rows<\/span>(<\/span>)<\/span> <<<\/span> "x"<\/span> <<<\/span> matrix.<\/span>cols<\/span>(<\/span>)<\/span> <<<\/span> ") ==="<\/span> <<<\/span> std::<\/span>endl;<\/span>
\n for<\/span> (<\/span>int<\/span> i =<\/span> 0<\/span>;<\/span> i <<\/span> matrix.<\/span>rows<\/span>(<\/span>)<\/span>;<\/span> ++<\/span>i)<\/span> {<\/span>
\n for<\/span> (<\/span>int<\/span> j =<\/span> 0<\/span>;<\/span> j <<\/span> matrix.<\/span>cols<\/span>(<\/span>)<\/span>;<\/span> ++<\/span>j)<\/span> {<\/span>
\n file <<<\/span> matrix<\/span>(<\/span>i,<\/span> j)<\/span>;<\/span>
\n if<\/span> (<\/span>j <<\/span> matrix.<\/span>cols<\/span>(<\/span>)<\/span> –<\/span> 1<\/span>)<\/span> file <<<\/span> ","<\/span>;<\/span>
\n }<\/span>
\n file <<<\/span> std::<\/span>endl;<\/span>
\n }<\/span>
\n file <<<\/span> std::<\/span>endl;<\/span>
\n }<\/span>
\n }<\/span><\/p>\n

static<\/span> void<\/span> save_excel<\/span>(<\/span>const<\/span> std::<\/span>string&<\/span> path,<\/span> const<\/span> EigenMap&<\/span> variables)<\/span> {<\/span>
\n logger-><\/span>warn<\/span>(<\/span>"Excel\u683c\u5f0f\u6682\u672a\u5b9e\u73b0: {}"<\/span>,<\/span> path)<\/span>;<\/span>
\n }<\/span><\/p>\n

static<\/span> void<\/span> save_parquet<\/span>(<\/span>const<\/span> std::<\/span>string&<\/span> path,<\/span> const<\/span> EigenMap&<\/span> variables)<\/span> {<\/span>
\n logger-><\/span>warn<\/span>(<\/span>"Parquet\u683c\u5f0f\u6682\u672a\u5b9e\u73b0: {}"<\/span>,<\/span> path)<\/span>;<\/span>
\n }<\/span>
\n}<\/span>;<\/span><\/p>\n

\/\/ M\u4ee3\u7801\u89e3\u6790\u5668<\/span>
\nclass<\/span> MCodeParser<\/span> {<\/span>
\npublic<\/span>:<\/span>
\n explicit<\/span> MCodeParser<\/span>(<\/span>const<\/span> std::<\/span>string&<\/span> script_path)<\/span> {<\/span>
\n parse_script<\/span>(<\/span>script_path)<\/span>;<\/span>
\n }<\/span><\/p>\n

void<\/span> execute<\/span>(<\/span>EigenMap&<\/span> variables)<\/span> const<\/span> {<\/span>
\n for<\/span> (<\/span>size_t i =<\/span> 0<\/span>;<\/span> i <<\/span> instructions.<\/span>size<\/span>(<\/span>)<\/span>;<\/span> ++<\/span>i)<\/span> {<\/span>
\n try<\/span> {<\/span>
\n const<\/span> auto<\/span>&<\/span> [<\/span>lhs,<\/span> rhs]<\/span> =<\/span> instructions[<\/span>i]<\/span>;<\/span>
\n MatrixXd result =<\/span> evaluate_expression<\/span>(<\/span>rhs,<\/span> variables)<\/span>;<\/span>
\n variables[<\/span>lhs]<\/span> =<\/span> result;<\/span>
\n logger-><\/span>info<\/span>(<\/span>"\u6267\u884c\u6210\u529f: {} = [{}x{}\u77e9\u9635] (\u884c {})"<\/span>,<\/span>
\n lhs,<\/span> result.<\/span>rows<\/span>(<\/span>)<\/span>,<\/span> result.<\/span>cols<\/span>(<\/span>)<\/span>,<\/span> i +<\/span> 1<\/span>)<\/span>;<\/span>
\n }<\/span> catch<\/span> (<\/span>const<\/span> std::<\/span>exception&<\/span> e)<\/span> {<\/span>
\n logger-><\/span>error<\/span>(<\/span>"\u6267\u884c\u5931\u8d25 (\u884c {}): {}\\\\n{}"<\/span>,<\/span>
\n i +<\/span> 1<\/span>,<\/span> instructions[<\/span>i]<\/span>.<\/span>first +<\/span> " = "<\/span> +<\/span> instructions[<\/span>i]<\/span>.<\/span>second,<\/span> e.<\/span>what<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n }<\/span>
\n }<\/span>
\n }<\/span><\/p>\n

private<\/span>:<\/span>
\n std::<\/span>vector<<\/span>std::<\/span>pair<<\/span>std::<\/span>string,<\/span> std::<\/span>string>><\/span> instructions;<\/span><\/p>\n

void<\/span> parse_script<\/span>(<\/span>const<\/span> std::<\/span>string&<\/span> path)<\/span> {<\/span>
\n std::<\/span>ifstream file<\/span>(<\/span>path)<\/span>;<\/span>
\n if<\/span> (<\/span>!<\/span>file.<\/span>is_open<\/span>(<\/span>)<\/span>)<\/span> THROW_EXCEPTION<\/span>(<\/span>"\u65e0\u6cd5\u6253\u5f00M\u6587\u4ef6: "<\/span> +<\/span> path)<\/span>;<\/span><\/p>\n

std::<\/span>string line;<\/span>
\n int<\/span> line_num =<\/span> 0<\/span>;<\/span>
\n while<\/span> (<\/span>std::<\/span>getline<\/span>(<\/span>file,<\/span> line)<\/span>)<\/span> {<\/span>
\n line_num++<\/span>;<\/span>
\n \/\/ \u79fb\u9664\u6ce8\u91ca\u548c\u7a7a\u767d<\/span>
\n size_t comment_pos =<\/span> line.<\/span>find<\/span>(<\/span>'%'<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>comment_pos !=<\/span> std::<\/span>string::<\/span>npos)<\/span> {<\/span>
\n line =<\/span> line.<\/span>substr<\/span>(<\/span>0<\/span>,<\/span> comment_pos)<\/span>;<\/span>
\n }<\/span>
\n line =<\/span> trim<\/span>(<\/span>line)<\/span>;<\/span>
\n if<\/span> (<\/span>line.<\/span>empty<\/span>(<\/span>)<\/span>)<\/span> continue<\/span>;<\/span><\/p>\n

\/\/ \u89e3\u6790\u8d4b\u503c\u8bed\u53e5<\/span>
\n size_t eq_pos =<\/span> line.<\/span>find<\/span>(<\/span>'='<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>eq_pos ==<\/span> std::<\/span>string::<\/span>npos)<\/span> {<\/span>
\n logger-><\/span>warn<\/span>(<\/span>"\u5ffd\u7565\u65e0\u6548\u8bed\u53e5 (\u884c {}): {}"<\/span>,<\/span> line_num,<\/span> line)<\/span>;<\/span>
\n continue<\/span>;<\/span>
\n }<\/span><\/p>\n

std::<\/span>string lhs =<\/span> trim<\/span>(<\/span>line.<\/span>substr<\/span>(<\/span>0<\/span>,<\/span> eq_pos)<\/span>)<\/span>;<\/span>
\n std::<\/span>string rhs =<\/span> trim<\/span>(<\/span>line.<\/span>substr<\/span>(<\/span>eq_pos +<\/span> 1<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>lhs.<\/span>empty<\/span>(<\/span>)<\/span> ||<\/span> rhs.<\/span>empty<\/span>(<\/span>)<\/span>)<\/span> {<\/span>
\n logger-><\/span>warn<\/span>(<\/span>"\u5ffd\u7565\u4e0d\u5b8c\u6574\u8d4b\u503c (\u884c {}): {}"<\/span>,<\/span> line_num,<\/span> line)<\/span>;<\/span>
\n continue<\/span>;<\/span>
\n }<\/span><\/p>\n

instructions.<\/span>emplace_back<\/span>(<\/span>lhs,<\/span> rhs)<\/span>;<\/span>
\n }<\/span>
\n logger-><\/span>info<\/span>(<\/span>"\u89e3\u6790\u5b8c\u6210: {} \u6761\u6709\u6548\u6307\u4ee4"<\/span>,<\/span> instructions.<\/span>size<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n }<\/span><\/p>\n

static<\/span> std::<\/span>string trim<\/span>(<\/span>const<\/span> std::<\/span>string&<\/span> s)<\/span> {<\/span>
\n auto<\/span> start =<\/span> std::<\/span>find_if_not<\/span>(<\/span>s.<\/span>begin<\/span>(<\/span>)<\/span>,<\/span> s.<\/span>end<\/span>(<\/span>)<\/span>,<\/span> ::<\/span>isspace)<\/span>;<\/span>
\n auto<\/span> end =<\/span> std::<\/span>find_if_not<\/span>(<\/span>s.<\/span>rbegin<\/span>(<\/span>)<\/span>,<\/span> s.<\/span>rend<\/span>(<\/span>)<\/span>,<\/span> ::<\/span>isspace)<\/span>.<\/span>base<\/span>(<\/span>)<\/span>;<\/span>
\n return<\/span> (<\/span>start <<\/span> end)<\/span> ?<\/span> std::<\/span>string<\/span>(<\/span>start,<\/span> end)<\/span> :<\/span> ""<\/span>;<\/span>
\n }<\/span><\/p>\n

MatrixXd evaluate_expression<\/span>(<\/span>const<\/span> std::<\/span>string&<\/span> expr,<\/span> const<\/span> EigenMap&<\/span> vars)<\/span> const<\/span> {<\/span>
\n \/\/ \u5904\u7406\u51fd\u6570\u8c03\u7528<\/span>
\n static<\/span> std::<\/span>regex func_re<\/span>(<\/span>R"((\\\\w+)\\\\s*\\\\(\\\\s*([^)]*)\\\\s*\\\\))"<\/span>)<\/span>;<\/span>
\n std::<\/span>smatch match;<\/span>
\n if<\/span> (<\/span>std::<\/span>regex_match<\/span>(<\/span>expr,<\/span> match,<\/span> func_re)<\/span>)<\/span> {<\/span>
\n std::<\/span>string func =<\/span> match[<\/span>1<\/span>]<\/span>;<\/span>
\n std::<\/span>vector<<\/span>std::<\/span>string><\/span> args =<\/span> split_arguments<\/span>(<\/span>match[<\/span>2<\/span>]<\/span>)<\/span>;<\/span>
\n return<\/span> call_function<\/span>(<\/span>func,<\/span> args,<\/span> vars)<\/span>;<\/span>
\n }<\/span><\/p>\n

\/\/ \u5904\u7406\u4e8c\u5143\u8fd0\u7b97<\/span>
\n static<\/span> std::<\/span>vector<<\/span>std::<\/span>pair<<\/span>std::<\/span>regex,<\/span> std::<\/span>function<<\/span>MatrixXd<\/span>(<\/span>const<\/span> MatrixXd&<\/span>,<\/span> const<\/span> MatrixXd&<\/span>)<\/span>>><\/span>><\/span> ops =<\/span> {<\/span>
\n {<\/span>std::<\/span>regex<\/span>(<\/span>R"(^\\\\s*(\\\\w+)\\\\s*\\\\+\\\\s*(\\\\w+)\\\\s*$)"<\/span>)<\/span>,<\/span> MatrixOps::<\/span>add}<\/span>,<\/span>
\n {<\/span>std::<\/span>regex<\/span>(<\/span>R"(^\\\\s*(\\\\w+)\\\\s*\\\\-\\\\s*(\\\\w+)\\\\s*$)"<\/span>)<\/span>,<\/span> MatrixOps::<\/span>subtract}<\/span>,<\/span>
\n {<\/span>std::<\/span>regex<\/span>(<\/span>R"(^\\\\s*(\\\\w+)\\\\s*\\\\*\\\\s*(\\\\w+)\\\\s*$)"<\/span>)<\/span>,<\/span> MatrixOps::<\/span>multiply}<\/span>,<\/span>
\n {<\/span>std::<\/span>regex<\/span>(<\/span>R"(^\\\\s*(\\\\w+)\\\\s*\\\\\/\\\\s*(\\\\w+)\\\\s*$)"<\/span>)<\/span>,<\/span> MatrixOps::<\/span>divide}<\/span>,<\/span>
\n {<\/span>std::<\/span>regex<\/span>(<\/span>R"(^\\\\s*(\\\\w+)\\\\s*\\\\*\\\\*\\\\s*(\\\\w+)\\\\s*$)"<\/span>)<\/span>,<\/span> [<\/span>]<\/span>(<\/span>const<\/span> MatrixXd&<\/span> a,<\/span> const<\/span> MatrixXd&<\/span> b)<\/span> {<\/span> return<\/span> MatrixOps<\/span>::<\/span>pow<\/span>(<\/span>a,<\/span> b<\/span>(<\/span>0<\/span>,<\/span>0<\/span>)<\/span>)<\/span>;<\/span> }<\/span>}<\/span>,<\/span>
\n }<\/span>;<\/span><\/p>\n

for<\/span> (<\/span>const<\/span> auto<\/span>&<\/span> [<\/span>op_re,<\/span> op_func]<\/span> :<\/span> ops)<\/span> {<\/span>
\n if<\/span> (<\/span>std::<\/span>regex_match<\/span>(<\/span>expr,<\/span> match,<\/span> op_re)<\/span>)<\/span> {<\/span>
\n std::<\/span>string a_name =<\/span> match[<\/span>1<\/span>]<\/span>;<\/span>
\n std::<\/span>string b_name =<\/span> match[<\/span>2<\/span>]<\/span>;<\/span><\/p>\n

if<\/span> (<\/span>!<\/span>vars.<\/span>count<\/span>(<\/span>a_name)<\/span>)<\/span> THROW_EXCEPTION<\/span>(<\/span>"\u672a\u5b9a\u4e49\u53d8\u91cf: "<\/span> +<\/span> a_name)<\/span>;<\/span>
\n if<\/span> (<\/span>!<\/span>vars.<\/span>count<\/span>(<\/span>b_name)<\/span>)<\/span> THROW_EXCEPTION<\/span>(<\/span>"\u672a\u5b9a\u4e49\u53d8\u91cf: "<\/span> +<\/span> b_name)<\/span>;<\/span><\/p>\n

return<\/span> op_func<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>a_name)<\/span>,<\/span> vars.<\/span>at<\/span>(<\/span>b_name)<\/span>)<\/span>;<\/span>
\n }<\/span>
\n }<\/span><\/p>\n

\/\/ \u5904\u7406\u4e00\u5143\u8fd0\u7b97<\/span>
\n static<\/span> std::<\/span>vector<<\/span>std::<\/span>pair<<\/span>std::<\/span>regex,<\/span> std::<\/span>function<<\/span>MatrixXd<\/span>(<\/span>const<\/span> MatrixXd&<\/span>)<\/span>>><\/span>><\/span> unary_ops =<\/span> {<\/span>
\n {<\/span>std::<\/span>regex<\/span>(<\/span>R"(^\\\\s*sin\\\\s*\\\\(\\\\s*(\\\\w+)\\\\s*\\\\)$)"<\/span>)<\/span>,<\/span> [<\/span>]<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> MatrixOps<\/span>::<\/span>sin<\/span>(<\/span>m)<\/span>;<\/span> }<\/span>}<\/span>,<\/span>
\n {<\/span>std::<\/span>regex<\/span>(<\/span>R"(^\\\\s*cos\\\\s*\\\\(\\\\s*(\\\\w+)\\\\s*\\\\)$)"<\/span>)<\/span>,<\/span> [<\/span>]<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> MatrixOps<\/span>::<\/span>cos<\/span>(<\/span>m)<\/span>;<\/span> }<\/span>}<\/span>,<\/span>
\n {<\/span>std::<\/span>regex<\/span>(<\/span>R"(^\\\\s*tan\\\\s*\\\\(\\\\s*(\\\\w+)\\\\s*\\\\)$)"<\/span>)<\/span>,<\/span> [<\/span>]<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> MatrixOps<\/span>::<\/span>tan<\/span>(<\/span>m)<\/span>;<\/span> }<\/span>}<\/span>,<\/span>
\n {<\/span>std::<\/span>regex<\/span>(<\/span>R"(^\\\\s*exp\\\\s*\\\\(\\\\s*(\\\\w+)\\\\s*\\\\)$)"<\/span>)<\/span>,<\/span> [<\/span>]<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> MatrixOps<\/span>::<\/span>exp<\/span>(<\/span>m)<\/span>;<\/span> }<\/span>}<\/span>,<\/span>
\n {<\/span>std::<\/span>regex<\/span>(<\/span>R"(^\\\\s*log\\\\s*\\\\(\\\\s*(\\\\w+)\\\\s*\\\\)$)"<\/span>)<\/span>,<\/span> [<\/span>]<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> MatrixOps<\/span>::<\/span>log<\/span>(<\/span>m)<\/span>;<\/span> }<\/span>}<\/span>,<\/span>
\n {<\/span>std::<\/span>regex<\/span>(<\/span>R"(^\\\\s*sqrt\\\\s*\\\\(\\\\s*(\\\\w+)\\\\s*\\\\)$)"<\/span>)<\/span>,<\/span> [<\/span>]<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> MatrixOps<\/span>::<\/span>sqrt<\/span>(<\/span>m)<\/span>;<\/span> }<\/span>}<\/span>,<\/span>
\n {<\/span>std::<\/span>regex<\/span>(<\/span>R"(^\\\\s*abs\\\\s*\\\\(\\\\s*(\\\\w+)\\\\s*\\\\)$)"<\/span>)<\/span>,<\/span> [<\/span>]<\/span>(<\/span>const<\/span> MatrixXd&<\/span> m)<\/span> {<\/span> return<\/span> MatrixOps<\/span>::<\/span>abs<\/span>(<\/span>m)<\/span>;<\/span> }<\/span>}<\/span>,<\/span>
\n }<\/span>;<\/span><\/p>\n

for<\/span> (<\/span>const<\/span> auto<\/span>&<\/span> [<\/span>op_re,<\/span> op_func]<\/span> :<\/span> unary_ops)<\/span> {<\/span>
\n if<\/span> (<\/span>std::<\/span>regex_match<\/span>(<\/span>expr,<\/span> match,<\/span> op_re)<\/span>)<\/span> {<\/span>
\n std::<\/span>string var_name =<\/span> match[<\/span>1<\/span>]<\/span>;<\/span>
\n if<\/span> (<\/span>!<\/span>vars.<\/span>count<\/span>(<\/span>var_name)<\/span>)<\/span> THROW_EXCEPTION<\/span>(<\/span>"\u672a\u5b9a\u4e49\u53d8\u91cf: "<\/span> +<\/span> var_name)<\/span>;<\/span>
\n return<\/span> op_func<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>var_name)<\/span>)<\/span>;<\/span>
\n }<\/span>
\n }<\/span><\/p>\n

\/\/ \u5904\u7406\u53d8\u91cf\u5f15\u7528<\/span>
\n if<\/span> (<\/span>vars.<\/span>count<\/span>(<\/span>expr)<\/span>)<\/span> {<\/span>
\n return<\/span> vars.<\/span>at<\/span>(<\/span>expr)<\/span>;<\/span>
\n }<\/span><\/p>\n

THROW_EXCEPTION<\/span>(<\/span>"\u65e0\u6cd5\u89e3\u6790\u8868\u8fbe\u5f0f: "<\/span> +<\/span> expr)<\/span>;<\/span>
\n }<\/span><\/p>\n

static<\/span> std::<\/span>vector<<\/span>std::<\/span>string><\/span> split_arguments<\/span>(<\/span>const<\/span> std::<\/span>string&<\/span> args_str)<\/span> {<\/span>
\n std::<\/span>vector<<\/span>std::<\/span>string><\/span> args;<\/span>
\n std::<\/span>istringstream iss<\/span>(<\/span>args_str)<\/span>;<\/span>
\n std::<\/span>string arg;<\/span>
\n while<\/span> (<\/span>std::<\/span>getline<\/span>(<\/span>iss,<\/span> arg,<\/span> ','<\/span>)<\/span>)<\/span> {<\/span>
\n arg =<\/span> trim<\/span>(<\/span>arg)<\/span>;<\/span>
\n if<\/span> (<\/span>!<\/span>arg.<\/span>empty<\/span>(<\/span>)<\/span>)<\/span> args.<\/span>push_back<\/span>(<\/span>arg)<\/span>;<\/span>
\n }<\/span>
\n return<\/span> args;<\/span>
\n }<\/span><\/p>\n

MatrixXd call_function<\/span>(<\/span>const<\/span> std::<\/span>string&<\/span> name,<\/span> const<\/span> std::<\/span>vector<<\/span>std::<\/span>string><\/span>&<\/span> args,<\/span> const<\/span> EigenMap&<\/span> vars)<\/span> const<\/span> {<\/span>
\n if<\/span> (<\/span>args.<\/span>empty<\/span>(<\/span>)<\/span>)<\/span> THROW_EXCEPTION<\/span>(<\/span>"\u51fd\u6570\u8c03\u7528\u7f3a\u5c11\u53c2\u6570: "<\/span> +<\/span> name)<\/span>;<\/span><\/p>\n

\/\/ \u68c0\u67e5\u6240\u6709\u53c2\u6570\u662f\u5426\u5b58\u5728<\/span>
\n for<\/span> (<\/span>const<\/span> auto<\/span>&<\/span> arg :<\/span> args)<\/span> {<\/span>
\n if<\/span> (<\/span>!<\/span>vars.<\/span>count<\/span>(<\/span>arg)<\/span>)<\/span> THROW_EXCEPTION<\/span>(<\/span>"\u51fd\u6570\u53c2\u6570\u672a\u5b9a\u4e49: "<\/span> +<\/span> arg)<\/span>;<\/span>
\n }<\/span><\/p>\n

\/\/ \u8c03\u7528\u5bf9\u5e94\u77e9\u9635\u51fd\u6570<\/span>
\n if<\/span> (<\/span>name ==<\/span> "inv"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>inv<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "pinv"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>pinv<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "eig"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>eig<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "svd"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>svd<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "chol"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>chol<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "sin"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>sin<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "cos"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>cos<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "tan"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>tan<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "exp"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>exp<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "log"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>log<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "log10"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>log10<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "log2"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>log2<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "sqrt"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>sqrt<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "abs"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>abs<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "sign"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>sign<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "round"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>round<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "floor"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>floor<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "ceil"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>ceil<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "fix"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>fix<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "sum"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>sum<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "mean"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>mean<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "prod"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>prod<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "max"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>max_val<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "min"<\/span> &&<\/span> args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 1<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>min_val<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>)<\/span>;<\/span><\/p>\n

\/\/ \u4e8c\u5143\u51fd\u6570<\/span>
\n if<\/span> (<\/span>args.<\/span>size<\/span>(<\/span>)<\/span> ==<\/span> 2<\/span>)<\/span> {<\/span>
\n if<\/span> (<\/span>name ==<\/span> "atan2"<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>atan2<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>,<\/span> vars.<\/span>at<\/span>(<\/span>args[<\/span>1<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "mod"<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>mod<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>,<\/span> vars.<\/span>at<\/span>(<\/span>args[<\/span>1<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "rem"<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>rem<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>,<\/span> vars.<\/span>at<\/span>(<\/span>args[<\/span>1<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "max"<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>max<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>,<\/span> vars.<\/span>at<\/span>(<\/span>args[<\/span>1<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>name ==<\/span> "min"<\/span>)<\/span> return<\/span> MatrixOps<\/span>::<\/span>min<\/span>(<\/span>vars.<\/span>at<\/span>(<\/span>args[<\/span>0<\/span>]<\/span>)<\/span>,<\/span> vars.<\/span>at<\/span>(<\/span>args[<\/span>1<\/span>]<\/span>)<\/span>)<\/span>;<\/span>
\n }<\/span><\/p>\n

THROW_EXCEPTION<\/span>(<\/span>"\u4e0d\u652f\u6301\u7684\u51fd\u6570\u6216\u53c2\u6570\u6570\u91cf\u9519\u8bef: "<\/span> +<\/span> name)<\/span>;<\/span>
\n }<\/span>
\n}<\/span>;<\/span><\/p>\n

\/\/ \u7ebf\u7a0b\u6c60\u5904\u7406\u5668<\/span>
\nclass<\/span> ThreadPool<\/span> {<\/span>
\npublic<\/span>:<\/span>
\n explicit<\/span> ThreadPool<\/span>(<\/span>size_t threads =<\/span> std::<\/span>thread::<\/span>hardware_concurrency<\/span>(<\/span>)<\/span>)<\/span>
\n :<\/span> stop<\/span>(<\/span>false<\/span>)<\/span> {<\/span>
\n for<\/span> (<\/span>size_t i =<\/span> 0<\/span>;<\/span> i <<\/span> threads;<\/span> ++<\/span>i)<\/span> {<\/span>
\n workers.<\/span>emplace_back<\/span>(<\/span>[<\/span>this<\/span>]<\/span> {<\/span>
\n while<\/span> (<\/span>true<\/span>)<\/span> {<\/span>
\n std::<\/span>function<<\/span>void<\/span>(<\/span>)<\/span>><\/span> task;<\/span>
\n {<\/span>
\n std::<\/span>unique_lock<<\/span>std::<\/span>mutex><\/span> lock<\/span>(<\/span>queue_mutex)<\/span>;<\/span>
\n condition.<\/span>wait<\/span>(<\/span>lock,<\/span> [<\/span>this<\/span>]<\/span> {<\/span> return<\/span> stop ||<\/span> !<\/span>tasks.<\/span>empty<\/span>(<\/span>)<\/span>;<\/span> }<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>stop &&<\/span> tasks.<\/span>empty<\/span>(<\/span>)<\/span>)<\/span> return<\/span>;<\/span>
\n task =<\/span> std::<\/span>move<\/span>(<\/span>tasks.<\/span>front<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n tasks.<\/span>pop<\/span>(<\/span>)<\/span>;<\/span>
\n }<\/span>
\n task<\/span>(<\/span>)<\/span>;<\/span>
\n }<\/span>
\n }<\/span>)<\/span>;<\/span>
\n }<\/span>
\n logger-><\/span>info<\/span>(<\/span>"\u7ebf\u7a0b\u6c60\u521d\u59cb\u5316\u5b8c\u6210，\u7ebf\u7a0b\u6570: {}"<\/span>,<\/span> workers.<\/span>size<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n }<\/span><\/p>\n

template<\/span><<\/span>class<\/span> F<\/span>><\/span>
\n void<\/span> enqueue<\/span>(<\/span>F&&<\/span> f)<\/span> {<\/span>
\n {<\/span>
\n std::<\/span>unique_lock<<\/span>std::<\/span>mutex><\/span> lock<\/span>(<\/span>queue_mutex)<\/span>;<\/span>
\n if<\/span> (<\/span>stop)<\/span> throw<\/span> std::<\/span>runtime_error<\/span>(<\/span>"\u7ebf\u7a0b\u6c60\u5df2\u505c\u6b62，\u65e0\u6cd5\u6dfb\u52a0\u4efb\u52a1"<\/span>)<\/span>;<\/span>
\n tasks.<\/span>emplace<\/span>(<\/span>std::<\/span>forward<\/span><<\/span>F><\/span><\/span><\/span>(<\/span>f)<\/span>)<\/span>;<\/span>
\n }<\/span>
\n condition.<\/span>notify_one<\/span>(<\/span>)<\/span>;<\/span>
\n }<\/span><\/p>\n

~<\/span>ThreadPool<\/span>(<\/span>)<\/span> {<\/span>
\n {<\/span>
\n std::<\/span>unique_lock<<\/span>std::<\/span>mutex><\/span> lock<\/span>(<\/span>queue_mutex)<\/span>;<\/span>
\n stop =<\/span> true<\/span>;<\/span>
\n }<\/span>
\n condition.<\/span>notify_all<\/span>(<\/span>)<\/span>;<\/span>
\n for<\/span> (<\/span>std::<\/span>thread&<\/span> worker :<\/span> workers)<\/span> {<\/span>
\n if<\/span> (<\/span>worker.<\/span>joinable<\/span>(<\/span>)<\/span>)<\/span> worker.<\/span>join<\/span>(<\/span>)<\/span>;<\/span>
\n }<\/span>
\n }<\/span><\/p>\n

private<\/span>:<\/span>
\n std::<\/span>vector<<\/span>std::<\/span>thread><\/span> workers;<\/span>
\n std::<\/span>queue<<\/span>std::<\/span>function<<\/span>void<\/span>(<\/span>)<\/span>>><\/span> tasks;<\/span>
\n std::<\/span>mutex queue_mutex;<\/span>
\n std::<\/span>condition_variable condition;<\/span>
\n bool<\/span> stop;<\/span>
\n}<\/span>;<\/span><\/p>\n

\/\/ \u6587\u4ef6\u5904\u7406\u51fd\u6570<\/span>
\nvoid<\/span> process_file<\/span>(<\/span>const<\/span> std::<\/span>string&<\/span> mat_path,<\/span>
\n const<\/span> std::<\/span>string&<\/span> m_script,<\/span>
\n const<\/span> std::<\/span>vector<<\/span>std::<\/span>string><\/span>&<\/span> output_formats)<\/span> {<\/span>
\n try<\/span> {<\/span>
\n logger-><\/span>info<\/span>(<\/span>"\u5f00\u59cb\u5904\u7406: {}"<\/span>,<\/span> mat_path)<\/span>;<\/span><\/p>\n

\/\/ \u52a0\u8f7dMAT\u6587\u4ef6\u53d8\u91cf<\/span>
\n EigenMap variables =<\/span> MatFileHandler<\/span>::<\/span>load<\/span>(<\/span>mat_path)<\/span>;<\/span>
\n logger-><\/span>info<\/span>(<\/span>"\u52a0\u8f7d\u53d8\u91cf\u6570\u91cf: {}"<\/span>,<\/span> variables.<\/span>size<\/span>(<\/span>)<\/span>)<\/span>;<\/span><\/p>\n

\/\/ \u89e3\u6790\u5e76\u6267\u884cM\u4ee3\u7801<\/span>
\n MCodeParser parser<\/span>(<\/span>m_script)<\/span>;<\/span>
\n parser.<\/span>execute<\/span>(<\/span>variables)<\/span>;<\/span><\/p>\n

\/\/ \u751f\u6210\u8f93\u51fa\u6587\u4ef6\u540d（\u79fb\u9664\u539f\u6269\u5c55\u540d）<\/span>
\n std::<\/span>string base_name =<\/span> fs::<\/span>path<\/span>(<\/span>mat_path)<\/span>.<\/span>stem<\/span>(<\/span>)<\/span>.<\/span>string<\/span>(<\/span>)<\/span>;<\/span><\/p>\n

\/\/ \u4fdd\u5b58\u7ed3\u679c\u5230\u6307\u5b9a\u683c\u5f0f<\/span>
\n for<\/span> (<\/span>const<\/span> auto<\/span>&<\/span> format :<\/span> output_formats)<\/span> {<\/span>
\n std::<\/span>string output_path =<\/span> base_name +<\/span> "."<\/span> +<\/span> format;<\/span>
\n if<\/span> (<\/span>format ==<\/span> "mat"<\/span>)<\/span> {<\/span>
\n MatFileHandler<\/span>::<\/span>save<\/span>(<\/span>output_path,<\/span> variables)<\/span>;<\/span>
\n }<\/span> else<\/span> if<\/span> (<\/span>format ==<\/span> "csv"<\/span>)<\/span> {<\/span>
\n FormatExporter<\/span>::<\/span>save_csv<\/span>(<\/span>output_path,<\/span> variables)<\/span>;<\/span>
\n }<\/span> else<\/span> if<\/span> (<\/span>format ==<\/span> "xlsx"<\/span>)<\/span> {<\/span>
\n FormatExporter<\/span>::<\/span>save_excel<\/span>(<\/span>output_path,<\/span> variables)<\/span>;<\/span>
\n }<\/span> else<\/span> if<\/span> (<\/span>format ==<\/span> "parquet"<\/span>)<\/span> {<\/span>
\n FormatExporter<\/span>::<\/span>save_parquet<\/span>(<\/span>output_path,<\/span> variables)<\/span>;<\/span>
\n }<\/span> else<\/span> {<\/span>
\n logger-><\/span>warn<\/span>(<\/span>"\u4e0d\u652f\u6301\u7684\u8f93\u51fa\u683c\u5f0f: {}"<\/span>,<\/span> format)<\/span>;<\/span>
\n continue<\/span>;<\/span>
\n }<\/span>
\n logger-><\/span>info<\/span>(<\/span>"\u5df2\u4fdd\u5b58\u5230: {}"<\/span>,<\/span> output_path)<\/span>;<\/span>
\n }<\/span><\/p>\n

logger-><\/span>info<\/span>(<\/span>"\u5904\u7406\u5b8c\u6210: {}"<\/span>,<\/span> mat_path)<\/span>;<\/span>
\n }<\/span> catch<\/span> (<\/span>const<\/span> std::<\/span>exception&<\/span> e)<\/span> {<\/span>
\n logger-><\/span>error<\/span>(<\/span>"\u5904\u7406\u6587\u4ef6\u5931\u8d25: {}\\\\n{}"<\/span>,<\/span> mat_path,<\/span> e.<\/span>what<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n }<\/span>
\n}<\/span><\/p>\n

\/\/ \u76ee\u5f55\u5904\u7406\u51fd\u6570<\/span>
\nvoid<\/span> process_directory<\/span>(<\/span>const<\/span> std::<\/span>string&<\/span> dir_path,<\/span>
\n const<\/span> std::<\/span>string&<\/span> m_script,<\/span>
\n const<\/span> std::<\/span>vector<<\/span>std::<\/span>string><\/span>&<\/span> output_formats,<\/span>
\n ThreadPool&<\/span> pool)<\/span> {<\/span>
\n if<\/span> (<\/span>!<\/span>fs::<\/span>is_directory<\/span>(<\/span>dir_path)<\/span>)<\/span> {<\/span>
\n logger-><\/span>error<\/span>(<\/span>"\u4e0d\u662f\u6709\u6548\u76ee\u5f55: {}"<\/span>,<\/span> dir_path)<\/span>;<\/span>
\n return<\/span>;<\/span>
\n }<\/span><\/p>\n

for<\/span> (<\/span>const<\/span> auto<\/span>&<\/span> entry :<\/span> fs::<\/span>directory_iterator<\/span>(<\/span>dir_path)<\/span>)<\/span> {<\/span>
\n if<\/span> (<\/span>entry.<\/span>is_regular_file<\/span>(<\/span>)<\/span> &&<\/span> entry.<\/span>path<\/span>(<\/span>)<\/span>.<\/span>extension<\/span>(<\/span>)<\/span> ==<\/span> ".mat"<\/span>)<\/span> {<\/span>
\n pool.<\/span>enqueue<\/span>(<\/span>[<\/span>=<\/span>]<\/span> {<\/span>
\n process_file<\/span>(<\/span>entry.<\/span>path<\/span>(<\/span>)<\/span>.<\/span>string<\/span>(<\/span>)<\/span>,<\/span> m_script,<\/span> output_formats)<\/span>;<\/span>
\n }<\/span>)<\/span>;<\/span>
\n }<\/span>
\n }<\/span>
\n logger-><\/span>info<\/span>(<\/span>"\u76ee\u5f55\u4efb\u52a1\u5df2\u5168\u90e8\u52a0\u5165\u961f\u5217: {}"<\/span>,<\/span> dir_path)<\/span>;<\/span>
\n}<\/span><\/p>\n

\/\/ \u547d\u4ee4\u884c\u53c2\u6570\u89e3\u6790<\/span>
\nstruct<\/span> CmdArgs<\/span> {<\/span>
\n std::<\/span>string input_path;<\/span>
\n std::<\/span>string script_path;<\/span>
\n std::<\/span>vector<<\/span>std::<\/span>string><\/span> output_formats;<\/span>
\n size_t threads =<\/span> std::<\/span>thread::<\/span>hardware_concurrency<\/span>(<\/span>)<\/span>;<\/span>
\n}<\/span>;<\/span><\/p>\n

CmdArgs parse_args<\/span>(<\/span>int<\/span> argc,<\/span> char<\/span>*<\/span> argv[<\/span>]<\/span>)<\/span> {<\/span>
\n if<\/span> (<\/span>argc <<\/span> 4<\/span>)<\/span> {<\/span>
\n std::<\/span>cerr <<<\/span> "\u7528\u6cd5: "<\/span> <<<\/span> argv[<\/span>0<\/span>]<\/span>
\n <<<\/span> " <\u8f93\u5165\u8def\u5f84(.mat\u6216\u76ee\u5f55)> <M\u811a\u672c\u8def\u5f84> <\u8f93\u51fa\u683c\u5f0f(\u7528\u9017\u53f7\u5206\u9694)> [\u7ebf\u7a0b\u6570]\\\\n"<\/span>
\n <<<\/span> "\u793a\u4f8b: "<\/span> <<<\/span> argv[<\/span>0<\/span>]<\/span>
\n <<<\/span> " data.mat calc.m mat,csv 4\\\\n"<\/span>
\n <<<\/span> "\u652f\u6301\u683c\u5f0f: mat, csv, xlsx(\u5f85\u5b9e\u73b0), parquet(\u5f85\u5b9e\u73b0)\\\\n"<\/span>;<\/span>
\n exit<\/span>(<\/span>1<\/span>)<\/span>;<\/span>
\n }<\/span><\/p>\n

CmdArgs args;<\/span>
\n args.<\/span>input_path =<\/span> argv[<\/span>1<\/span>]<\/span>;<\/span>
\n args.<\/span>script_path =<\/span> argv[<\/span>2<\/span>]<\/span>;<\/span><\/p>\n

\/\/ \u89e3\u6790\u8f93\u51fa\u683c\u5f0f<\/span>
\n std::<\/span>istringstream iss<\/span>(<\/span>argv[<\/span>3<\/span>]<\/span>)<\/span>;<\/span>
\n std::<\/span>string format;<\/span>
\n while<\/span> (<\/span>std::<\/span>getline<\/span>(<\/span>iss,<\/span> format,<\/span> ','<\/span>)<\/span>)<\/span> {<\/span>
\n if<\/span> (<\/span>!<\/span>format.<\/span>empty<\/span>(<\/span>)<\/span>)<\/span> args.<\/span>output_formats.<\/span>push_back<\/span>(<\/span>format)<\/span>;<\/span>
\n }<\/span><\/p>\n

\/\/ \u89e3\u6790\u7ebf\u7a0b\u6570<\/span>
\n if<\/span> (<\/span>argc >=<\/span> 5<\/span>)<\/span> {<\/span>
\n try<\/span> {<\/span>
\n args.<\/span>threads =<\/span> std::<\/span>stoul<\/span>(<\/span>argv[<\/span>4<\/span>]<\/span>)<\/span>;<\/span>
\n if<\/span> (<\/span>args.<\/span>threads ==<\/span> 0<\/span>)<\/span> args.<\/span>threads =<\/span> 1<\/span>;<\/span>
\n }<\/span> catch<\/span> (<\/span>.<\/span>.<\/span>.<\/span>)<\/span> {<\/span>
\n logger-><\/span>warn<\/span>(<\/span>"\u65e0\u6548\u7ebf\u7a0b\u6570，\u4f7f\u7528\u9ed8\u8ba4\u503c: {}"<\/span>,<\/span> args.<\/span>threads)<\/span>;<\/span>
\n }<\/span>
\n }<\/span><\/p>\n

\/\/ \u9a8c\u8bc1\u8f93\u5165<\/span>
\n if<\/span> (<\/span>!<\/span>fs::<\/span>exists<\/span>(<\/span>args.<\/span>input_path)<\/span>)<\/span> {<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"\u8f93\u5165\u8def\u5f84\u4e0d\u5b58\u5728: "<\/span> +<\/span> args.<\/span>input_path)<\/span>;<\/span>
\n }<\/span>
\n if<\/span> (<\/span>!<\/span>fs::<\/span>exists<\/span>(<\/span>args.<\/span>script_path)<\/span> ||<\/span> fs::<\/span>path<\/span>(<\/span>args.<\/span>script_path)<\/span>.<\/span>extension<\/span>(<\/span>)<\/span> !=<\/span> ".m"<\/span>)<\/span> {<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"\u65e0\u6548M\u811a\u672c\u6587\u4ef6: "<\/span> +<\/span> args.<\/span>script_path)<\/span>;<\/span>
\n }<\/span>
\n if<\/span> (<\/span>args.<\/span>output_formats.<\/span>empty<\/span>(<\/span>)<\/span>)<\/span> {<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"\u672a\u6307\u5b9a\u8f93\u51fa\u683c\u5f0f"<\/span>)<\/span>;<\/span>
\n }<\/span><\/p>\n

return<\/span> args;<\/span>
\n}<\/span><\/p>\n

int<\/span> main<\/span>(<\/span>int<\/span> argc,<\/span> char<\/span>*<\/span> argv[<\/span>]<\/span>)<\/span> {<\/span>
\n try<\/span> {<\/span>
\n logger-><\/span>info<\/span>(<\/span>"\u7a0b\u5e8f\u542f\u52a8"<\/span>)<\/span>;<\/span><\/p>\n

\/\/ \u89e3\u6790\u547d\u4ee4\u884c\u53c2\u6570<\/span>
\n CmdArgs args =<\/span> parse_args<\/span>(<\/span>argc,<\/span> argv)<\/span>;<\/span>
\n logger-><\/span>info<\/span>(<\/span>"\u8f93\u5165\u8def\u5f84: {}"<\/span>,<\/span> args.<\/span>input_path)<\/span>;<\/span>
\n logger-><\/span>info<\/span>(<\/span>"M\u811a\u672c: {}"<\/span>,<\/span> args.<\/span>script_path)<\/span>;<\/span>
\n logger-><\/span>info<\/span>(<\/span>"\u8f93\u51fa\u683c\u5f0f: {}"<\/span>,<\/span> fmt::<\/span>format<\/span>(<\/span>"{}"<\/span>,<\/span> fmt::<\/span>join<\/span>(<\/span>args.<\/span>output_formats,<\/span> ","<\/span>)<\/span>)<\/span>)<\/span>;<\/span>
\n logger-><\/span>info<\/span>(<\/span>"\u7ebf\u7a0b\u6570: {}"<\/span>,<\/span> args.<\/span>threads)<\/span>;<\/span><\/p>\n

\/\/ \u521d\u59cb\u5316\u7ebf\u7a0b\u6c60<\/span>
\n ThreadPool pool<\/span>(<\/span>args.<\/span>threads)<\/span>;<\/span><\/p>\n

\/\/ \u5904\u7406\u8f93\u5165（\u6587\u4ef6\u6216\u76ee\u5f55）<\/span>
\n if<\/span> (<\/span>fs::<\/span>is_regular_file<\/span>(<\/span>args.<\/span>input_path)<\/span> &&<\/span>
\n fs::<\/span>path<\/span>(<\/span>args.<\/span>input_path)<\/span>.<\/span>extension<\/span>(<\/span>)<\/span> ==<\/span> ".mat"<\/span>)<\/span> {<\/span>
\n \/\/ \u5904\u7406\u5355\u4e2aMAT\u6587\u4ef6<\/span>
\n pool.<\/span>enqueue<\/span>(<\/span>[<\/span>&<\/span>]<\/span> {<\/span>
\n process_file<\/span>(<\/span>args.<\/span>input_path,<\/span> args.<\/span>script_path,<\/span> args.<\/span>output_formats)<\/span>;<\/span>
\n }<\/span>)<\/span>;<\/span>
\n }<\/span> else<\/span> if<\/span> (<\/span>fs::<\/span>is_directory<\/span>(<\/span>args.<\/span>input_path)<\/span>)<\/span> {<\/span>
\n \/\/ \u5904\u7406\u76ee\u5f55\u4e0b\u6240\u6709MAT\u6587\u4ef6<\/span>
\n process_directory<\/span>(<\/span>args.<\/span>input_path,<\/span> args.<\/span>script_path,<\/span> args.<\/span>output_formats,<\/span> pool)<\/span>;<\/span>
\n }<\/span> else<\/span> {<\/span>
\n THROW_EXCEPTION<\/span>(<\/span>"\u8f93\u5165\u5fc5\u987b\u662f.mat\u6587\u4ef6\u6216\u76ee\u5f55: "<\/span> +<\/span> args.<\/span>input_path)<\/span>;<\/span>
\n }<\/span><\/p>\n

\/\/ \u7b49\u5f85\u6240\u6709\u4efb\u52a1\u5b8c\u6210（\u901a\u8fc7\u7ebf\u7a0b\u6c60\u6790\u6784\u5b9e\u73b0）<\/span>
\n logger-><\/span>info<\/span>(<\/span>"\u7b49\u5f85\u6240\u6709\u4efb\u52a1\u5b8c\u6210…"<\/span>)<\/span>;<\/span>
\n }<\/span> catch<\/span> (<\/span>const<\/span> std::<\/span>exception&<\/span> e)<\/span> {<\/span>
\n logger-><\/span>critical<\/span>(<\/span>"\u7a0b\u5e8f\u5f02\u5e38\u7ec8\u6b62: {}"<\/span>,<\/span> e.<\/span>what<\/span>(<\/span>)<\/span>)<\/span>;<\/span>
\n return<\/span> 1<\/span>;<\/span>
\n }<\/span><\/p>\n

logger-><\/span>info<\/span>(<\/span>"\u7a0b\u5e8f\u6b63\u5e38\u7ed3\u675f"<\/span>)<\/span>;<\/span>
\n return<\/span> 0<\/span>;<\/span>
\n}<\/span><\/p>\n

\u8fd9\u4e2a\u5b8c\u6574\u7684GCC C++\u4ee3\u7801\u5b9e\u73b0\u4e86\u4ee5\u4e0b\u529f\u80fd：<\/p>\n

\u4e3b\u8981\u7279\u6027<\/h3>\n
  • \n

    \u77e9\u9635\u8fd0\u7b97\u652f\u6301：<\/p>\n