গণনা বিজ্ঞান

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

12 linear-solver  c++  lapack 

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

12 multigrid 

আমি ভেরিয়েবল এবং এক্স 2 এর জন্য নীচের সমীকরণগুলির সিস্টেমটি সমাধান করার চেষ্টা করছি (অন্য সবগুলি ধ্রুবক):পি, এক্স1P,x1P,x_1এক্স2x2x_2 এ ( 1 - পি)2- কে1এক্স1= 0এ পি2- কে2এক্স2= 0( 1 - পি) ( আর1+ এক্স1)4এল1- পি( আর1+ এক্স2)4এল2= 0A(1−P)2−k1x1=0AP2−k2x2=0(1−P)(r1+x1)4L1−P(r1+x2)4L2=0\frac{A(1-P)}{2}-k_1x_1=0 \\ \frac{AP}{2}-k_2x_2=0 \\ \frac{(1-P)(r_1+x_1)^4}{L_1}-\frac{P(r_1+x_2)^4}{L_2}=0 আমি দেখতে পাচ্ছি যে আমি এই সিস্টেমের সমীকরণগুলিকে …

12 optimization  convergence  newton-method 

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

12 linear-algebra  petsc  linear-solver  sparse-matrix 

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

12 pde 

আমি এমন একটি সমস্যা সমাধানের জন্য ম্যাটল্যাব ব্যবহার করছি যা প্রতিটি টাইমস্টেপে 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 

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

12 fluid-dynamics  finite-volume  data-structures 

আমি একটি জটিল ইন্টিগ্রাল নিয়ে কাজ করছি যা শূন্যের কাছাকাছি নির্দিষ্ট মানগুলিতে NaN গুলি প্রদর্শন করে এবং এই মুহূর্তে আমি তাদের সাথে বেশ কুরুচিপূর্ণভাবে একটি ইসনান বিবৃতি ব্যবহার করছি যা ঘটনার সময় শূন্যের সাথে সংহতকরণ নির্ধারণ করে। আমি এটি ফরটারনে এনএমএস লাইব্রেরি দিয়ে (কিউ 1 ডি রুটিন - কিউ 1 …

12 quadrature 

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

12 parallel-computing  monte-carlo  gpu 

আমি এমন একটি লাইব্রেরি খুঁজছি যা সংখ্যার স্থিতিশীলতার জন্য আত্মবিশ্বাসের বৃহত স্পার্স ম্যাট্রিক্সে ম্যাট্রিক্স অপারেশন করে। ম্যাট্রিকগুলি 1000+ বাই 1000+ এবং ম্যাট্রিক্সের মান 0 এবং 1000 এর মধ্যে হবে I আমি সূচক ক্যালকুলাস আলগোরিদম সম্পাদন করব তাই আমি ম্যাট্রিক্সের ক্রমিক (স্পার্স) সারি ভেক্টরগুলি ক্রমিক উত্পন্ন করব। আমি প্রতিটি সারিতে বিকাশ …

12 linear-algebra  algorithms  matrix 

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

12 simulation 

আমি 3 ডি-তে 4-নোড এলিমেন্টের মাধ্যমে বহুবচনীয় ভাবটি সংহত করতে চাই। এফআইএর বেশ কয়েকটি বইয়ের ক্ষেত্রে এমন একটি মামলা রয়েছে যেখানে একটি স্বেচ্ছাসেবী সমতল 4-নন উপাদানকে সংহত করা হয়। এই ক্ষেত্রে স্বাভাবিক প্রক্রিয়াটি হ'ল জ্যাকোব্যাটিক ম্যাট্রিক্স সন্ধান করা এবং এটি নির্ধারণকারীকে ইন্টিগ্রেশন ভিত্তিকে সাধারণীকরণের পরিবর্তনে ব্যবহার করতে হয় যার মধ্যে …

12 finite-element  quadrature 

Nocedal & রাইট বই মতে সংখ্যাসূচক অপ্টিমাইজেশান (2006), একটি বেঠিক লাইন অনুসন্ধানের জন্য Wolfe এর অবস্থার একটি বংশদ্ভুত নির্দেশনার জন্য রয়েছে ,ppp পর্যাপ্ত হ্রাস: বক্রতা শর্ত: ∇ f ( x + α p ) টি পি ≥ সি 2 ∇ ফ ( এক্স ) টি পি জন্য 0 < গ …

12 optimization 

আমি আশ্চর্য হই: যেখানে আসল ম্যাট্রিক্স সমাধানের জন্য সেরা অ্যালগরিদম কী । এ স্পষ্টভাবে সময় নির্ভর নয়, সাধারণত বিরল তবে অগত্যা ব্যান্ডেড হয় না। এর ইগেনুয়ালুগুলির অ-ইতিবাচক বাস্তব অংশ রয়েছে। এটি তির্যকীয়ও তবে সম্পূর্ণ তির্যক গণনার জন্য দক্ষ হওয়ার পক্ষে এটি খুব বড় হতে পারে।dudt=Aududt=Au\begin{equation} \frac{du}{dt} = Au \end{equation}AAAn×nn×nn\times n …

12 linear-algebra  ode 

আমার নিম্নোক্ত অপ্টিমাইজেশান সমস্যা রয়েছে যেখানে আমার সীমাবদ্ধতার মধ্যে আমার নিখুঁত মান রয়েছে: x∈Rnx∈Rn\mathbf{x} \in \mathbb{R}^nf0,f1,…,fmf0,f1,…,fm\mathbf{f}_0, \mathbf{f}_1, \ldots, \mathbf{f}_mnnnmins.t.fT0x|fT1x|≤|fT2x|≤…≤|fTmx|minf0Txs.t.|f1Tx|≤|f2Tx|≤…≤|fmTx|\begin{align} \min &\mathbf{f}_0^T \mathbf{x} \notag \\ \text{s.t.} &|\mathbf{f}_1^T \mathbf{x}| \leq |\mathbf{f}_2^T \mathbf{x}| \leq \ldots \leq |\mathbf{f}_m^T \mathbf{x}| \end{align} আমি জানি যে সম্ভাব্য স্থানটি উত্তল হবে না এবং সমস্যাটি সমাধান করার জন্য আমার …

12 optimization  linear-programming