| 5 | 1/1 | 返回列表 |
| 查看: 2659 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
pmsl666至尊木蟲 (著名寫手)
|
[求助]
FreeFem++ 怎么將所有單元的所有節(jié)點的坐標(biāo)及位移值輸出 已有1人參與
|
|
|
我使用FreeFem++ 求解一個偏微分方程,采用的是P2單元,P2單元按理說應(yīng)該是6節(jié)點三角形單元,我想按單元順序,將所有單元的所有節(jié)點的坐標(biāo)及位移值輸出, 我看了幫助文檔,貌似里面只有三角形3個頂點坐標(biāo)及位移輸出,怎么樣才能將三角形頂點以及中點坐標(biāo)位移都輸出到txt文件中。特此求助。 { ofstream ff("grapht.txt" ); for (int i=0;i<Th.nt;i++) { for (int j=0; j <3; j++) ff<<Th[j].x << " "<< Th[j].y<< " "<<u[][Vh(i,j)]<<endl; ff<<Th[0].x << " "<< Th[0].y<< " "<<u[][Vh(i,0)]<<endl <<endl<<endl; } } [ Last edited by feixiaolin on 2015-6-12 at 16:00 ] |

銀蟲 (正式寫手)
新蟲 (初入文壇)
{ ofstream ff("graph1.txt" ;for (int i=0;i<Th.nt;i++) {ff<<Th[0].x << " "<< Th[0].y<< " "<<u(Th[0].x,Th[0].y)<<endl; ff<<(Th[0].x+Th[1].x)/2 << " "<< (Th[0].y+Th[1].y)/2<< " "<<u((Th[0].x+Th[1].x)/2,(Th[0].y+Th[1].y)/2)<<endl; ff<<Th[1].x << " "<< Th[1].y<< " "<<u(Th[1].x,Th[1].y)<<endl; ff<<(Th[1].x+Th[2].x)/2 << " "<< (Th[1].y+Th[2].y)/2<< " "<<u((Th[1].x+Th[2].x)/2,(Th[1].y+Th[2].y)/2)<<endl; ff<<Th[2].x << " "<< Th[2].y<< " "<<u(Th[2].x,Th[2].y)<<endl; ff<<(Th[0].x+Th[2].x)/2 << " "<< (Th[0].y+Th[2].y)/2<< " "<<u((Th[0].x+Th[2].x)/2,(Th[0].y+Th[2].y)/2)<<endl; ff<<Th[0].x << " "<< Th[0].y<< " "<<u(Th[0].x,Th[0].y)<<endl <<endl<<endl; } } u(x,y)就能直接輸出[x,y]點的位移,邊界中點坐標(biāo)可以根據(jù)三角形頂點求得 |

新蟲 (初入文壇)
新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿哈工大,085400,320,求調(diào)劑 +3 | gdlf9999 2026-03-24 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|
[考研] 291求調(diào)劑 +8 | hhhhxn.. 2026-03-23 | 8/400 |
|
|
[考研] 265求調(diào)劑 +10 | 梁梁校校 2026-03-17 | 10/500 |
|
|
[考研] 招08考數(shù)學(xué) +6 | laoshidan 2026-03-20 | 14/700 |
|
|
[考研] 306求調(diào)劑 +9 | chuanzhu川燭 2026-03-18 | 9/450 |
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|