| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 2668 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | ||
pmsl666至尊木蟲(chóng) (著名寫(xiě)手)
|
[求助]
FreeFem++ 怎么將所有單元的所有節(jié)點(diǎn)的坐標(biāo)及位移值輸出 已有1人參與
|
|
|
我使用FreeFem++ 求解一個(gè)偏微分方程,采用的是P2單元,P2單元按理說(shuō)應(yīng)該是6節(jié)點(diǎn)三角形單元,我想按單元順序,將所有單元的所有節(jié)點(diǎn)的坐標(biāo)及位移值輸出, 我看了幫助文檔,貌似里面只有三角形3個(gè)頂點(diǎn)坐標(biāo)及位移輸出,怎么樣才能將三角形頂點(diǎn)以及中點(diǎn)坐標(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 ] |

新蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
{ 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]點(diǎn)的位移,邊界中點(diǎn)坐標(biāo)可以根據(jù)三角形頂點(diǎn)求得 |

新蟲(chóng) (初入文壇)
銀蟲(chóng) (正式寫(xiě)手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[基金申請(qǐng)] 面上5B能上會(huì)嗎? +4 | redcom 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
[考研] 求收留 +5 | 1943443204 2026-03-28 | 5/250 |
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚(yú) 2026-03-26 | 9/450 |
|
|
[考研] 一志愿北化085600材料專(zhuān)碩275|有文章專(zhuān)利|求調(diào)劑 +7 | Micky11223 2026-03-25 | 7/350 |
|
|
[考研] 085602 化工專(zhuān)碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
[考研] 311求調(diào)劑 +3 | 希望上岸阿小楊 2026-03-23 | 3/150 |
|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
10+4
|
chibby 2026-03-25 | 4/200 |
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
[考研] 085600,材料與化工321分,求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
[考研]
材料調(diào)劑
5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 材料專(zhuān)碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
[考研] 282求調(diào)劑 +3 | wcq131415 2026-03-24 | 3/150 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專(zhuān)碩) +4 | allen-yin 2026-03-23 | 6/300 |
|