99爱在线视频这里只有精品_窝窝午夜看片成人精品_日韩精品久久久毛片一区二区_亚洲一区二区久久

合肥生活安徽新聞合肥交通合肥房產(chǎn)生活服務(wù)合肥教育合肥招聘合肥旅游文化藝術(shù)合肥美食合肥地圖合肥社保合肥醫(yī)院企業(yè)服務(wù)合肥法律

AERO20542代做、代寫Python/Java編程

時(shí)間:2024-03-07  來源:合肥網(wǎng)hfw.cc  作者:hfw.cc 我要糾錯(cuò)



MECH20042/AERO20542 Numerical Methods and Computing
Laboratory exercise 1: Direct methods for the solution of
tridiagonal systems of linear equations
Solution of systems of linear equations is one of the most frequently encountered problems in
numerical modelling and simulation. Efficient numerical methods, both in terms of the execution time
and memory storage are essential to complete this task. Sparse systems of linear equations arise in
many applications, such as finite element or finite volume solution of differential equations. Sparse
linear systems have coefficient matrices that are sparse, i.e., a large proportion of the elements are
equal to zero. Banded matrices are a special class of sparse matrices in which the non-zero coefficients
are concentrated about the main diagonal.
Storing sparse matrices in computer memory as two-dimensional arrays is inefficient, as many zero
elements are kept needlessly in computer memory. Banded matrices can be stored by their diagonals,
where each diagonal is stored as a one-dimensional array (a vector). With this setup a tridiagonal
matrix 𝑇 of size 𝑛 × 𝑛

can be stored using three vectors as follows:
𝐴 = [𝑎11 𝑎22 ⋯ 𝑎𝑛𝑛]
𝑇 ∈ 𝑅
𝑛
,
w**; = [𝑎21 𝑎** ⋯ 𝑎𝑛,𝑛−1]
𝑇 ∈ 𝑅
𝑛−1
,
𝐶 = [𝑎12 𝑎23 ⋯ 𝑎𝑛−1,𝑛]
𝑇 ∈ 𝑅
𝑛−1
.
The Gaussian elimination technique applied to a tridiagonal system 𝑇𝒙 = 𝒇 is particularly simple,
because only the non-zero elements in the sub-diagonal held in vector w**; need to be eliminated. This
algorithm, known as the Thomas algorithm, proceeds as follows:
FORWARD ELIMINATION BACKSUBSTITUTION
𝑎𝑖𝑖 = 𝑎𝑖𝑖 −
𝑎𝑖,𝑖−1
𝑎𝑖−1,𝑖−1
𝑎𝑖−1,𝑖 w**9;𝑛 =
𝑓𝑛
𝑎𝑛𝑛
𝑓𝑖 = 𝑓𝑖 −
𝑎𝑖,𝑖−1
𝑎𝑖−1,𝑖−1
𝑓𝑖−1 w**9;𝑖 =
1
𝑎𝑖𝑖
(𝑓𝑖 − 𝑎𝑖,𝑖+1 w**9;𝑖+1)
𝑖 = 2, … , 𝑛 𝑖 = 𝑛 − 1, … ,1
TASK 1. Calculate the number of arithmetic operations that are required to solve a tridiagonal system
𝑇𝒙 = 𝒇 of size 𝑛 using the Thomas algorithm. Based on this result, determine the asymptotic
complexity of the Thomas algorithm, and compare it to the asymptotic complexity of the standard
Gaussian elimination.
TASK 2. Rewrite the Thomas algorithm in terms of the arrays 𝐴,w**;, and 𝐶 introduced to store the matrix
𝑇 efficiently.
TASK 3. Implement the Thomas algorithm from TASK 2 as a Python function. The input parameters to
the function should be the coefficient matrix 𝑇 (stored as three arrays 𝐴,w**;, and 𝐶) and the right-hand
side vector 𝒇. The output should be the solution vector 𝒙. The coefficient matrix and the right-hand
side should be defined in the main script and passed to the function that solves the system.
TASK 4. Test your code by solving the linear system of size 𝑛 = 10 with the values 𝐴 = 2, and w**; = 𝐶 =
−1. Set the right-hand side to 𝒇 = 𝟏. To verify the correctness of your code, compare the solution
vector obtained from the Thomas algorithm to that obtained by applying the direct solver
numpy.linalg.solve(). For the latter, the coefficient matrix should be assembled.
TASK 5. Solve five linear systems 𝑇𝒙 = 𝒇 with 𝐴 = 2, w**; = 𝐶 = −1 and 𝒇 = 𝟏 varying the problem size
𝑛 between 106
and 108
. Record the execution times in seconds for each case. To accomplish this task,
explore the Python function timer() from the package timeit (refer to the code for matrix
multiplication covered in lectures). Plot a graph where the obtained execution times are represented
as the function of the problem size 𝑛. What are your conclusions about the cost of the Thomas
請加QQ:99515681  郵箱:99515681@qq.com   WX:codehelp

掃一掃在手機(jī)打開當(dāng)前頁
  • 上一篇:PROG2007代寫、Python/c++程序語言代做
  • 下一篇:代寫CMSC 323、代做Java/Python編程
  • 無相關(guān)信息
    合肥生活資訊

    合肥圖文信息
    急尋熱仿真分析?代做熱仿真服務(wù)+熱設(shè)計(jì)優(yōu)化
    急尋熱仿真分析?代做熱仿真服務(wù)+熱設(shè)計(jì)優(yōu)化
    出評 開團(tuán)工具
    出評 開團(tuán)工具
    挖掘機(jī)濾芯提升發(fā)動機(jī)性能
    挖掘機(jī)濾芯提升發(fā)動機(jī)性能
    海信羅馬假日洗衣機(jī)亮相AWE  復(fù)古美學(xué)與現(xiàn)代科技完美結(jié)合
    海信羅馬假日洗衣機(jī)亮相AWE 復(fù)古美學(xué)與現(xiàn)代
    合肥機(jī)場巴士4號線
    合肥機(jī)場巴士4號線
    合肥機(jī)場巴士3號線
    合肥機(jī)場巴士3號線
    合肥機(jī)場巴士2號線
    合肥機(jī)場巴士2號線
    合肥機(jī)場巴士1號線
    合肥機(jī)場巴士1號線
  • 短信驗(yàn)證碼 豆包 幣安下載 AI生圖 目錄網(wǎng)

    關(guān)于我們 | 打賞支持 | 廣告服務(wù) | 聯(lián)系我們 | 網(wǎng)站地圖 | 免責(zé)聲明 | 幫助中心 | 友情鏈接 |

    Copyright © 2025 hfw.cc Inc. All Rights Reserved. 合肥網(wǎng) 版權(quán)所有
    ICP備06013414號-3 公安備 42010502001045

    99爱在线视频这里只有精品_窝窝午夜看片成人精品_日韩精品久久久毛片一区二区_亚洲一区二区久久

          亚洲免费在线精品一区| 久久精品国产亚洲精品| 一区二区三区成人| 99国产精品国产精品毛片| 亚洲精品视频中文字幕| 亚洲免费激情| 99热这里只有精品8| aa国产精品| 亚洲一区二区三区在线看| 亚洲视频一二区| 欧美中文字幕在线观看| 久久免费精品视频| 欧美经典一区二区三区| 欧美日韩日本网| 国产精品一区久久| 国内精品模特av私拍在线观看| 国产一区二区三区在线观看精品| 好吊妞这里只有精品| 亚洲片在线资源| 亚洲尤物视频网| 久久国产精品久久精品国产| 美女91精品| 欧美性猛交xxxx乱大交蜜桃| 国产人成精品一区二区三| 在线观看成人小视频| 亚洲免费观看视频| 香蕉久久夜色精品| 欧美成人精品在线播放| 国产精品久久久久一区二区三区| 国产亚洲高清视频| 99精品视频免费在线观看| 欧美在线一区二区| 欧美另类一区二区三区| 国产视频一区在线| 日韩一级视频免费观看在线| 久久黄金**| 欧美日韩黄色大片| 在线观看日韩av先锋影音电影院| 在线视频日韩| 牛人盗摄一区二区三区视频| 国产精品视频yy9099| 亚洲狠狠婷婷| 久久av一区| 国产精品九色蝌蚪自拍| 亚洲激情小视频| 久久精品夜色噜噜亚洲aⅴ| 欧美日韩国产综合久久| 一区二区在线视频| 欧美一级在线播放| 欧美视频网址| 亚洲精品美女在线观看| 美女精品网站| 国内自拍亚洲| 欧美专区福利在线| 国产精品免费观看视频| 在线亚洲观看| 欧美日韩免费观看一区三区| 亚洲激情精品| 免费不卡欧美自拍视频| 好看的日韩视频| 久久久久久9999| 国产一区白浆| 欧美在线视频不卡| 国产一区二区av| 欧美一区中文字幕| 国产欧美二区| 欧美一区二区三区四区在线观看地址| 欧美网站在线观看| 亚洲一区在线观看免费观看电影高清| 欧美日韩不卡一区| 99成人在线| 国产精品成人一区二区三区吃奶| 99热精品在线观看| 欧美日韩一区二区欧美激情| 亚洲精品网站在线播放gif| 欧美精品在线一区二区| 99这里只有久久精品视频| 欧美日韩一区在线观看视频| aa国产精品| 欧美午夜精品理论片a级大开眼界| 99成人精品| 国产精品日本一区二区| 久久se精品一区精品二区| 韩国av一区二区| 欧美承认网站| 日韩一级视频免费观看在线| 亚洲欧美成人一区二区在线电影| 国产精品hd| 亚洲性感美女99在线| 国产精品成人播放| 性色一区二区| 国产私拍一区| 久久国产福利国产秒拍| 国产综合av| 麻豆久久久9性大片| 一区在线免费| 欧美经典一区二区| 在线一区二区三区做爰视频网站 | 男女精品网站| 亚洲国产精品va在线观看黑人| 久久精品视频播放| 亚洲激情电影在线| 欧美日韩成人综合| 亚洲一级网站| 国产午夜精品理论片a级探花| 久久九九国产精品怡红院| 在线成人国产| 国产精品jizz在线观看美国 | 亚洲老板91色精品久久| 欧美日韩免费网站| 亚洲欧美色一区| 在线观看日韩专区| 欧美视频一区| 久久久久久精| 亚洲国产精品ⅴa在线观看| 国产精品毛片va一区二区三区| 亚洲欧美美女| 亚洲国产欧美日韩另类综合| 欧美日韩国产首页| 久久精品国产免费| 亚洲精品国产精品国自产在线| 欧美色视频在线| 美女被久久久| 亚洲欧美制服另类日韩| 亚洲成在线观看| 国产精品xxxav免费视频| 久久精品一二三区| 亚洲一区欧美激情| 91久久视频| 国产日韩精品在线| 性做久久久久久久久| 中文亚洲欧美| 亚洲国产99精品国自产| 国产视频亚洲精品| 欧美调教视频| 欧美激情在线狂野欧美精品| 亚洲免费中文字幕| 亚洲一区二区三区高清 | 亚洲国产精品免费| 国产美女精品视频免费观看| 欧美日韩精品一区二区三区| 久久亚洲国产精品一区二区| 亚洲欧美第一页| 一区二区三区免费观看| 在线观看亚洲一区| 1024国产精品| 樱桃视频在线观看一区| 国产亚洲一区二区三区在线观看| 欧美午夜精品久久久久久超碰| 欧美成人免费观看| 欧美 日韩 国产在线| 久久高清免费观看| 亚洲欧美在线观看| 先锋影院在线亚洲| 午夜久久99| 欧美一级视频| 久久精品导航| 美女任你摸久久| 久久人人爽爽爽人久久久| 久久久久9999亚洲精品| 欧美一区二区三区日韩| 亚洲午夜久久久久久久久电影网| 亚洲精品欧美日韩专区| 亚洲精品一区二区在线观看| 亚洲欧洲三级| 一区二区三区精品国产| 在线综合+亚洲+欧美中文字幕| 亚洲人成网站精品片在线观看| 这里只有精品视频| 亚洲欧美怡红院| 久久九九99视频| 麻豆成人在线观看| 欧美国产三区| 久久都是精品| 欧美精品在线观看| 国产精品porn| 国产一区二区三区久久精品| 激情综合激情| 一区二区三区在线视频免费观看 | 国产精品久久影院| 国产亚洲aⅴaaaaaa毛片| 精品999在线观看| 日韩午夜在线观看视频| 中文日韩在线| 久久精品亚洲一区二区| 欧美大片在线观看一区二区| 欧美视频一区二| 国内精品视频666| 一区二区在线免费观看| 亚洲一区二区三区免费观看| 久久国产日本精品| 免费不卡中文字幕视频| 欧美成年人网站| 欧美极品在线观看| 国产亚洲观看| 日韩亚洲精品电影| 久久久久久久91| 国产精品日韩在线播放| 亚洲国产精品久久久久久女王| 中日韩午夜理伦电影免费|