در levmar یک پیاده سازی از غیر خطی لونبرگ-مارکوارت حداقل الگوریتم مربع در C / C ++ می باشد
در
از معمول lmder از Minpack، اجرا در '80s اولیه در آزمایشگاه ملی آرگون، شاید اجرای رایگان الگوریتم LM به طور گسترده استفاده می شود. lmder در FORTRAN77 نوشته شده است و در طول سال ها ثابت کرده است که یک قطعه قابل اعتماد از نرم افزار است. با توجه به اینکه فورترن را می توان از C / C ++ نامیده می شود، ممکن است در مورد انگیزه برای نوشتن یک نسخه از LM در C. خوب تعجب، مشکل این است که زمانی که FORTRAN از C نامیده می شود، برنامه نویس باید آگاه باشند (و مطابق با ) قوانین متعددی در ارتباط با نام mangling، عبور استدلال، طرح حافظه آرایه چند بعدی، همایش ها ارتباط، و غیره، که غیر طبیعی هستند، نسبت به قوانین عادی C. دلیل دوم این است که این روش طول می کشد برای مسلم است که یک کامپایلر فرترن برای محیط زیست برنامه نویسی هدف در دسترس است، که ممکن است لزوما مورد باشد. یکی دیگر از دلایل است که با عدم درک عملکرد داخلی اجرای FORTRAN: گهگاه، زمانی که لازم است به دقت درک آنچه کد فرترن می کند، قطعه خاصی از آن ممکن است بدون هیچ گونه دانش فرترن به نظر می رسد غیر قابل درک به برنامه نویسان. فرترن به صورت خودکار به مترجمان C (به عنوان مثال f2c) مشکل را حل نمی کند چرا که کد C تولید شده است بسیار ناخوانا به انسان "uninitiated". علاوه بر این، مستندات توصیف ریاضیات که بر اساس آن پیاده سازی شده است بر اساس ممکن است مبهم یا غیر قابل دسترس. تاریخ و زمان آخرین اما نه کم، پیاده سازی LM نامزد در C باید آزاد و فنی صدا باشد. به عنوان مثال، نوع C از الگوریتم LM ارائه شده در کتاب "دستور عددی" (به عنوان مثال mrqmin)، است که همیشه یک انتخاب عملی: علاوه بر آن که حق کپی رایت، آن است که مشهور به عدم استحکام
در
از. به دلایل بالا، من بسته levmar که شامل پیاده سازی C از طعم LM که همچنین قابل استفاده با C ++ هستند توسعه یافته اند. levmar شامل پیاده سازی LM دقت دو جداره و تک، هر دو با اختلاف تحلیلی و محدود تقریبی Jacobians. این ارائه شده است رایگان، تحت شرایط مجوز GNU است. نظریه ریاضی پشت levmar نامحدود است که در جزئیات در یادداشت های سخنرانی مواد و روش ها تحت عنوان غیر خطی حداقل مربعات مشکلات توصیف، توسط K. مدسن، HB نیلسن و O. Tingleff، دانشگاه فنی دانمارک. پیاده سازی نرم افزار Matlab از الگوریتم های ارائه شده در یادداشت های سخنرانی نیز در دسترس هستند. توجه داشته باشید با این حال که تدوین مشکل به حداقل رساندن اینجا به تصویب رسید کمی متفاوت از که شرح داده شده در یادداشت های سخنرانی است استفاده
در
از تابع: کنید.
levmar ارائه می دهد چندین کاربر صدازدنی توابع اطاعت از زیر نامگذاری: در این نامه برای اولین بار (D یا بازدید کنندگان) مشخص با دقت مضاعف و یا تک و پسوند (_der یا _dif) نشان دهنده ژاکوبین تحلیلی یا تقریبی. اگر حال حاضر، LEC، سال قبل از میلاد و اجزای blec دلالت معادله خطی، جعبه و جعبه به طور همزمان و محدودیت های معادله خطی بود. بیشتر به طور خاص، levmar شامل توابع زیر:
در
در بهینه سازی نامحدود
در
در dlevmar_der (): دقت مضاعف، ژاکوبین تحلیلی
در dlevmar_dif (): دقت مضاعف، تفاضل محدود تقریبی slevmar_der ژاکوبین
از (): دقت تک، تحلیلی slevmar_dif
از ژاکوبین (): دقت تک، تفاضل محدود تقریبی ژاکوبین
در بهتر
بهینه سازی مقید
در
در dlevmar_lec_der () : دقت مضاعف، محدودیت های معادله خطی، ژاکوبین تحلیلی
در dlevmar_lec_dif (): دقت مضاعف، محدودیت های معادله خطی، تفاضل محدود تقریبی ژاکوبین
در slevmar_lec_der (): دقت تک، محدودیت های معادله خطی،
از ژاکوبین تحلیلی slevmar_lec_dif (): دقت تک، محدودیت های معادله خطی، تفاضل محدود تقریبی ژاکوبین
در
در dlevmar_bc_der (): دقت مضاعف، جعبه محدودیت، ژاکوبین تحلیلی
در dlevmar_bc_dif (): دقت مضاعف، جعبه محدودیت، محدود تفاوت تقریبی ژاکوبین
در slevmar_bc_der (): دقت تک، جعبه محدودیت، تحلیلی slevmar_bc_dif
از ژاکوبین (): دقت تک، جعبه محدودیت، تفاضل محدود تقریبی ژاکوبین
در
در dlevmar_blec_der (): دو دقت، جعبه و محدودیت های معادله خطی، تحلیلی dlevmar_blec_dif
از ژاکوبین (): دقت مضاعف، جعبه و محدودیت های معادله خطی، تفاضل محدود تقریبی ژاکوبین
در slevmar_blec_der (): دقت تک، جعبه و معادله خطی محدودیت، ژاکوبین تحلیلی
در slevmar_blec_dif (): دقت تک، جعبه و معادله خطی محدودیت، تفاضل محدود تقریبی ژاکوبین
در
در توجه داشته باشید که با استفاده از اختلاف محدود به تقریب نتایج ژاکوبین در ارزیابی تکراری از تابع به نصب می شود. به منظور کاهش تعداد کل این ارزیابی، توابع xxxxxxx_dif اجرای تقریب قاطع به ژاکوبین با استفاده از رتبه یک به روز رسانی Broyden است. همه توابع حل مشکل همان، یعنی آنها به دنبال P بردار پارامتر که به بهترین وجه توصیف (از نظر هنجار L2) اندازه گیری بردار x. بطور دقیقتر، با توجه به یک تابع بردار F: R ^ متر -> R ^ N با N> = M، آنها محاسبه AP به طوری که F (P) ~ = X، یعنی هنجار مربع || E || ^ 2 = | | XF (ص) || ^ 2 به حداقل رسیده است. همچنین، جعبه محدودیت از پوند فرم های [i]
چه جدید در این نسخه است:
در این نسخه اضافه می کند پشتیبانی از پوسته پوسته شدن قطر به xlevmar_bc_der ()، که می تواند همگرایی زمانی که بهبود نقطه شروع است به دور از minimizer درست است.
سیستم های حل کننده خطی حمایت از تجزیه موازی Cholesky با پلاسما، کتابخانه جبر خطی برای پردازنده های چند هسته ای.
حل کننده خطی اند ثابت شده است به طوری که آنها در مثلث کمتر از ماتریس متقارن، که در نتیجه در عملکرد حافظه نهان بهتر عمل می کنند.
فایل های پیکربندی CMake برای ساخت و ساز پروژه تجدید نظر شده.
چند تغییرات جزئی دیگر ساخته شده اند.
احتیاط برای نسخه های کاربران قبلی::
جدید در نسخه 2.5 است که برای هماهنگی بهتر و اجتناب از درگیری، نام فایل هدر levmar از lm.h به levmar تغییر یافت. ساعت. لطفا بر این اساس به روز رسانی فایل های منبع خود را.
اضافه شدن پشتیبانی برای به حداقل رساندن تحت جعبه به طور همزمان، معادله خطی و محدودیت های نابرابری (توابع xlevmar_bleic_der را ببینید () و xlevmar_bleic_dif ()).
فراغ راحتی اجرا به xlevmar_bleic_der () و xlevmar_bleic_dif () برخورد با موارد ساده تر محدود، به عنوان مثال xlevmar_blic_der () و xlevmar_blic_dif () به حداقل رساندن پشتیبانی زیر فقط جعبه و نابرابری خطی محدودیت.
اضافه شده سیستم های حل کننده خطی بر اساس تجزیه UDUt (Cholesky به عنوان مثال، SQRT-رایگان).
حذف برخی از کپی کردن حافظه های غیر ضروری را از بیشتر حل خطی.
اضافه شده مشکل آزمون داده اتصالات جدید (آزبورن).
نویسنده چند تغییرات جزئی دیگر.
احتیاط برای نسخه های قبلی اعضا:: اندازه استدلال اطلاعات (به عنوان مثال، LM_INFO_SZ) توسط یکی افزایش یافته بود به جای یک میدان بازگشت جدید (به عنوان مثال
جدید در نسخه 2.4 است ، اطلاعات [9]) که به تعداد کل سیستمهای خطی حل در طول این دوره از به حداقل رساندن مربوط. شاخص های زمینه قبلا بازگشت در آرایه اطلاعات بدون تغییر باقی می ماند.
یک طرح کش-کارآمد تر برای محاسبه بلند J تقریبی ^ T * J و J ^ T * E برای مشکلات به حداقل رساندن کوچک اجرا شود.
اطمینان می داد که هر حافظه کاری حفظ بین فراخوانی از حل کننده خطی است بر فسخ روال levmar منتشر شد.
اطمینان می داد که حل کننده خطی استخدام حداقل مقدار مورد نیاز از حافظه کمکی، اجتناب از گاه به گاه بیش از تخصیص.
یک طرح أنوندنج دستگاه حلقه برای بالا بردن سرعت محاسبه E = X-HX اجرا شده است.
یک زن و شوهر از مسائل با هم ترازی حافظه در سیستم های 64 بیتی ثابت.
توابع اضافه شده برای محاسبه ضریب تعیین.
سر و کار داشتن با مسائل چند جزئی.
نظر یافت نشد