কোন কোয়াড-ডাবল পাটিগণিত স্পারস ম্যাট্রিক্স প্যাকেজ আছে কি?


10

আমি কিছু অসুস্থ-শর্তযুক্ত বিশাল স্পারস লিনিয়ার সিস্টেমের সমীকরণের জন্য কাজ করছি। আমি তাদের সমাধানের জন্য ডাবল-ডাবল পাটিগণিত বা কোয়াড-ডাবল পাটিগণিত ব্যবহার করতে চাই। আমি জানি যে এমপিএকেকে নামে একটি প্যাকেজ রয়েছে যা নাকাটা, মাহো দ্বারা বিকাশ করা হয়েছে, যা কোয়াড-ডাবল গাণিতিকের অধীনে সংখ্যিক লিনিয়ার বীজগণিত গণনা সম্পাদন করতে পারে। তবে এটি ঘন ম্যাট্রিক্সের জন্য ডিজাইন করা হয়েছে, স্পার্স ম্যাট্রিক্সের জন্য নয়। আপনি কি জানেন যে কোয়াড-ডাবল পাটিগণিত স্পার্স ম্যাট্রিক্স প্যাকেজ আছে কি না?


আপনি কোন ধরণের ম্যাট্রিক নিয়ে কাজ করছেন? এটি কি সিমেট্রিক, হার্মিটিয়ান, ইতিবাচক সুনির্দিষ্ট? আপনি কি একটি স্পারস এলইউ করতে চান, বা পুনরুক্তি পদ্ধতি ব্যবহার করতে চান?
ভিক্টর লিউ

উত্তর:


3

সংস্করণ ৩.২ অনুসারে , পিইটিএসসি জিসিসি / গফর্ট্রান ৪.6 এবং আরও নতুন ক্ষেত্রে স্পার্স কোয়াড-স্পষ্টতা গণনা সমর্থন করে।

আপনার একটি কোয়াড-নির্ভুলতা BLAS এবং ল্যাপাকের দরকার হবে যা পিইটিএসসি আপনাকে নিম্নলিখিত (আংশিক) কনফিগার কমান্ড দিয়ে (কোয়াড সমর্থন সহ) সরবরাহ করতে পারে:

./configure --with-precision=__float128 --download-f2cblaslapack

আরও কিছু তথ্যের জন্য FAQ দেখুন ।

এছাড়াও, আমি nOOb এর সাথে একমত , যদি সম্ভব হয় তবে কোয়াড-নির্ভুলতায় সরিয়ে নেওয়ার আগে সিস্টেমটিকে নিয়মিত করার চেষ্টা করুন।


1
এটি কোয়াড নির্ভুলতা (128-বিট রিয়েলস), কোয়াড-ডাবল নয় (256-বিট রিয়েলস)। এতে বলা হয়েছে, কোয়াড নির্ভুলতা সাধারণত ডাবল স্পষ্টতা সহকারে স্থিতিশীলতার সমস্যাগুলি বোঝার জন্য যথেষ্ট এবং আপনি সাধারণত সিস্টেমটি স্কেল করে বিচক্ষণ করতে চান যাতে ডাবল নির্ভুলতা উত্পাদনের জন্য যথেষ্ট।
জেদ ব্রাউন

2

সংখ্যাযুক্ত নির্ভুলতা বাড়ানোর চেয়ে নিয়মিতকরণের মাধ্যমে অসুস্থ শর্তযুক্ত সিস্টেমগুলি আরও ভাল সমাধান করা হয়। বেহাল বিবরণগুলির জন্য "নিয়ন্ত্রিত অসুস্থ-পোজ" অনুসন্ধান করুন।


2

আমি ট্রিলিনোস গ্রন্থাগারটি চেষ্টা করে দেখতে পারি। তারা টেপেট্রার অধীনে স্পার্স ম্যাট্রিক্স লাইব্রেরি টেম্পলেট করেছে (যা তাদের মূল স্পার্স ম্যাট্রিক্স লাইব্রেরি এপেট্রাকে প্রতিস্থাপন করবে)। আপনি দ্বিগুণ, জটিল, চতুষ্কোণ ইত্যাদি টেমপ্লেট করতে পারেন এবং তাদের সম্ভবত পিইটিএসসি-এর পাশের সলভারগুলির (সবচেয়ে সরাসরি এবং পুনরুক্তি উভয়) বৃহত্তম নির্বাচন রয়েছে।

সম্পাদনা: মন্তব্যগুলি পড়ার পরে, টেপেট্রার তাত্ক্ষণিক উপযোগিতা মনে হচ্ছে কিছুটা প্রশ্নবিদ্ধ আর্ট কোয়াড নির্ভুলতা ...


দুর্ভাগ্যক্রমে দুর্ভাগ্যক্রমে: বেশিরভাগ সলভার এখনও টেপেট্রার সাথে কাজ করেন না:
ওল্ফগ্যাং

এটা দুর্ভাগ্যজনক. টেপেরা বিকাশ কতটা এগিয়ে গেছে সে সম্পর্কে আমি আশাবাদী কিন্তু অনিশ্চিত ছিলাম (সুতরাং "এপেট্রাকে প্রতিস্থাপনের কথা": পি)। আমি ভেবেছিলাম কমপক্ষে বেলোস লাইব্রেরি (যেমন ত্রিলিনোস সলভারগুলি তৃতীয় পক্ষের কোডগুলির চারপাশে মোড়কের উপর নির্ভর করে না) যদিও টেপেট্রাকে সমর্থন করে?
জেসি চ্যান

আমি মনে করি একটি বেলোস 2 আছে। নিশ্চিতভাবেই, ত্রিলিনোস প্রকল্পটি তাদের উত্সগুলি টেপেট্রার পিছনে ফেলেছে এবং এটি ভবিষ্যতে মানক প্যাকেজ হবে। আমি মনে করি যদিও তারা এখনও সেখানে নেই।
ওল্ফগ্যাং ব্যাঙ্গারথ

এমনকি __float128টেপেরাও সমর্থিত হওয়ার বিষয়ে আমি অবগত নই , চতুর্থাংশ দ্বিগুণ। Tpetra Stand-Alone নয় এবং সব হেডার, এবং এমনকি যদি এটা ছিল, ভালো জিনিস std::complexএকমাত্র কাজ floatএবং double
জেড ব্রাউন

1

ম্যাটল্যাবের জন্য মাল্টিপ্রেসিশন কম্পিউটিং টুলবাক্সে বিচ্ছিন্ন ম্যাট্রিক্সের জন্য সমর্থন রয়েছে এবং চতুর্ভুজ নির্ভুলতার সাথে গণনার জন্য বিশেষত অনুকূলিত হয়েছে।

কোয়াড-স্পষ্টতা স্পার্স সলভারগুলির জন্য সময় বিবরণ এখানে: স্পার্স ম্যাট্রিক্সের জন্য সরাসরি সমাধানকারী

(আমি টুলবক্সের লেখক)।

আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.