প্রশ্ন ট্যাগ «linear-solver»

সমাধান করার জন্য আমার কাছে বেশ কয়েকটি চ্যালেঞ্জিং অ-উত্তেজক গ্লোবাল অপ্টিমাইজেশান সমস্যা রয়েছে। বর্তমানে আমি ম্যাটল্যাবের অপটিমাইজেশন টুলবক্সfmincon()'sqp' ব্যবহার করছি (বিশেষত, অ্যালগরিদম = সহ ), যা বেশ কার্যকর । তবে আমার বেশিরভাগ কোডটি পাইথনে রয়েছে এবং আমি পাইথনেও অপ্টিমাইজেশন করতে পছন্দ করি। পাইথন বাইন্ডিংগুলির সাথে প্রতিযোগিতা করতে পারে এমন কোনও …

77 python  optimization  nonlinear-programming  linear-algebra  linear-solver  sparse-matrix  ode  automatic-differentiation  efficiency  matrix  symbolic-computation  pde  multiphysics  operator-splitting  fluid-dynamics  pde  multigrid  pde  hyperbolic-pde  quadrature  precision  numerical-analysis  fluid-dynamics  interpolation  c  ode  testing  data-sets  pde  fluid-dynamics  hyperbolic-pde  pde  finite-element  fluid-dynamics  stability  incompressible  pde  hyperbolic-pde  fluid-dynamics  pde  precision  floating-point  convex-optimization  conditioning  pde  hyperbolic-pde  stability  implicit-methods  explicit-methods  fluid-dynamics  accuracy  pde  hyperbolic-pde  petsc  linear-algebra 

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

49 linear-algebra  linear-solver  sparse-matrix 

আমার জ্ঞান অনুসারে, রৈখিক সমীকরণের একটি পদ্ধতি সমাধানের জন্য 4 টি উপায় রয়েছে (যদি আরও থাকে তবে আমাকে সংশোধন করুন): যদি সিস্টেম ম্যাট্রিক্স একটি পূর্ণ-র‌্যাঙ্ক বর্গ ম্যাট্রিক্স হয় তবে আপনি ক্র্যামারস বিধি ব্যবহার করতে পারেন; সিস্টেম ম্যাট্রিক্সের বিপরীত বা সিউডোইনভার গণনা করুন; ম্যাট্রিক্স পচানোর পদ্ধতিগুলি ব্যবহার করুন (গাউসিয়ান বা গাউস-জর্ডান …

31 linear-algebra  linear-solver  iterative-method 

আমার এবং ম্যাট্রিক রয়েছে । বিক্ষিপ্ত হয় এবং সঙ্গে খুব বড় (কয়েক মিলিয়ন অনুক্রম হতে পারে।) একটি হল সঙ্গে লম্বা ম্যাট্রিক্স বরং ছোট ( ) এবং প্রতিটি কলামের করতে পারেন শুধুমাত্র একটি একক আছে বাকি হচ্ছে এন্ট্রি এর, যেমন যে । বিশাল, তাই এটি উল্টানো সত্যিই শক্ত, এবং আমি মতো …

22 linear-algebra  linear-solver  sparse-matrix  krylov-method 

এখানে অনেকগুলি বিভিন্ন লাইব্রেরি রয়েছে যা একটি বিচ্ছিন্ন রৈখিক সিস্টেমের সমীকরণের সমাধান করে, তবে আমি পার্থক্যগুলি কী তা নির্ধারণ করতে অসুবিধে করছি। আমি যতদূর বলতে পারি তিনটি বড় প্যাকেজ রয়েছে: ত্রিলিনোস , পিইটিএসসি এবং ইন্টেল এমকেএল । তারা সকলেই স্পার্স ম্যাট্রিক্স সলভ করতে পারে, তারা সবাই দ্রুত (যতদূর আমি বলতে …

21 linear-algebra  linear-solver  libraries 

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

17 linear-solver  conjugate-gradient  gmres 

মাল্টিলেভেল ইনভার্স-বেসড আইএলইউ পূর্বশর্তীদের বিশেষত ভিন্ন ভিন্ন জেনারেল হেলমহোল্টজের সিরিয়াল অভিনয়ে আমি অত্যন্ত অভিভূত , তবে কোনও উন্মুক্ত উত্স বাস্তবায়ন খুঁজে না পেয়ে আমি অবাক হয়েছি। বিশেষত, ILUPACK বাইনারিগুলি একাডেমিকদের জন্য অবাধে উপলভ্য করে তোলে তবে এটি তাদের উত্স কোড প্রকাশ করে না বলে মনে হয় না। সত্যিই কি কেস …

15 pde  algorithms  linear-solver  libraries  preconditioning 

মনে করুন আমাদের একটি লিনিয়ার সিস্টেম রয়েছে এবং আমরা এর কন্ডিশনার সম্পর্কে কিছুই জানি না এবং সমাধান সম্পর্কে প্রাথমিক কোনও তথ্য নেই। আমরা অন্ধভাবে গাউসিয়ান নির্মূল প্রয়োগ করি এবং কিছু সমাধান xxx । এই সমাধানটি ম্যাট্রিক্সের প্রাথমিক প্রাথমিক বিশ্লেষণ ছাড়াই বিশ্বাসযোগ্য (অর্থাত্ সিস্টেমটি ভালভাবে কন্ডিশনারযুক্ত) কিনা তা নির্ধারণ করা সম্ভব …

14 linear-solver  condition-number  conditioning 

মাল্টিগ্রিড (এমজি) একটি প্রাথমিক অনুমান x 0 নির্ধারণ করে এবং i = 0 , 1 এর জন্য নিম্নলিখিতটি পুনরাবৃত্তি করে লিনিয়ার সিস্টেম সমাধানের জন্য ব্যবহার করা যেতে পারে .. রূপান্তর পর্যন্ত:Ax=bAx=bAx=bx0x0x_0i=0,1..i=0,1..i=0,1.. অবশিষ্টাংশ গণনা কর i ri=b−Axiri=b−Axir_i = b-Ax_i একটি আনুমানিক পেতে একটি মাল্টিগ্রিড চক্র প্রয়োগ করুন , যেখানে A e …

13 linear-algebra  linear-solver  multigrid  krylov-method 

আমি ভাবছি যে টমাস অ্যালগরিদমটি অ্যালগোরিদমিক জটিলতার ক্ষেত্রে (LAPACK ইত্যাদির মতো বাস্তবায়ন প্যাকেজগুলির সন্ধান না করে) একসম্মত তির্যকভাবে স্পার্স ট্রাইডিজোনাল সিস্টেমকে প্রাধান্য দেওয়ার সবচেয়ে দ্রুততম উপায় (সম্ভবত)? আমি জানি যে টমাস অ্যালগরিদম এবং মাল্টিগ্রিড উভয়ই জটিলতা, তবে সম্ভবত মাল্টিগ্রিডের জন্য ধ্রুবক ফ্যাক্টরটি কম? মাল্টিগ্রিড দ্রুততর হতে পারে বলে আমার কাছে …

13 linear-solver  sparse-matrix  complexity  multigrid 

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

13 matrix  linear-solver  precision 

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

13 linear-algebra  linear-solver  preconditioning  krylov-method  gmres 

আমি থেকে সি ++ তে একটি বিদ্যমান কোড পোর্ট করছি এবং সমাধান করার জন্য একটি রৈখিক ব্যবস্থা আছে (আরও সাধারণ ফর্ম চেয়ে বেশি )xA=bxA=bxA=bAx=bAx=bAx=b ম্যাট্রিক্স ঘন, এবং সাধারণ ফর্ম, কিন্তু কোন 1000x1000 চেয়ে বড়। সুতরাং ম্যাটল্যাবে, সমাধানটি ফাংশন বা ফরোয়ার্ড-স্ল্যাশ স্বরলিপি দ্বারা পাওয়া যায়AAAmrdivide(b,A)x = b/A; BLAS এবং LAPACK রুটিন …

12 linear-solver  c++  lapack 

আমার অনেকগুলি ডান হাতের (300 থেকে 1000) সমান বিচ্ছিন্ন রৈখিক সিস্টেম (300x300 থেকে 1000x1000) সমাধান করতে হবে। এই প্রথম সমস্যাটি ছাড়াও, আমি বিভিন্ন সিস্টেমগুলি সমাধান করতে চাই, তবে একই অ-শূন্য উপাদানগুলির সাথে (কেবলমাত্র বিভিন্ন মান), এটি ধ্রুবক স্পারসিটি প্যাটার্ন সহ অনেকগুলি স্পার সিস্টেম। আমার ম্যাট্রিকগুলি অনির্দিষ্ট। ফ্যাক্টরাইজেশন এবং ইনিশিয়ালাইজেশনের পারফরম্যান্স …

12 linear-algebra  petsc  linear-solver  sparse-matrix 

আমি এমন একটি সমস্যা সমাধানের জন্য ম্যাটল্যাব ব্যবহার করছি যা প্রতিটি টাইমস্টেপে solving সমাধান করে, যেখানে। Mathbf time সময়ের সাথে পরিবর্তিত হয়এখনই, আমি ম্যাটল্যাব এর ব্যবহার করে এটি সম্পাদন করছি :খএ x = খAx=b\mathbf{A} \mathbf{x}=\mathbf{b}খb\mathbf{b}mldivide x = A\b আমার যতটা প্রয়োজন পূর্বরূপ তৈরি করার নমনীয়তা রয়েছে তাই আমি ভাবছি যে …

12 linear-algebra  matlab  linear-solver  matrix  linear-system