máy bơm nước hỏa tiễn công ty in ấn https://duan-knparadise.com/ Bơm nước hỏa tiễn bóp da nam Túi nilon PP cong ty in bieu mau
Tính toán song song cho bài toán FEA ? - Page 2 - Diễn đàn của các kỹ sư kết cấu Việt Nam










































































Features
Go Back   Diễn đàn của các kỹ sư kết cấu Việt Nam > ỨNG DỤNG TIN HỌC TRONG THIẾT KẾ XÂY DỰNG > Lập trình tính toán, phân tích kết cấu
Tên thành viên
Mật mã
Tin HOT Thông tin BQT diễn đàn

Result  Reload Thống Kê - Diễn đàn của các kỹ sư kết cấu Việt Nam
All Forums | Thông báo | Văn Bản Pháp Luật | Văn Hóa Văn Nghệ Giao Lưu
Please wait...
Please wait...
Loading...
Trả lời
 
Ðiều Chỉnh Xếp Bài
Old 17-04-2006   #21
tnlinh
Thành viên
 
tnlinh's Avatar
 
Tham gia ngày: Nov 2004
Bài gởi: 958
Thanks: 0
Thanked 873 Times in 366 Posts
Send a message via Yahoo to tnlinh
Default Ðề: Tính toán song song cho bài toán FEA ?

Đúng là trong phương pháp lặp đó, rất cần cái ma trận Precondition. Người ta chứng minh được là trong bài toán kết cấu, luôn hội tụ, (hay có lời giải). Một ông bạn nói với tôi là thường lặp đến lần thứ 5 là hội tụ (nếu dùng cái ma trận precondition, tất nhiên cách chọn cũng là một việc cần bàn), trường hợp bất đắc dĩ thì có thể đến 1000 lần. Cái phương pháp lặp (Conjugate Gradients thuần túy) tôi có làm thử vào mùa hè vừa rồi, kết quả tệ lắm.

Tôi có tính thử một số phương pháp lặp trong module SparseMtx v1.5.1 thấy rằng, nếu dùng tính lặp, có thể tăng tốc độ lên 3 lần so với phương pháp trực tiếp. Trong module đó có 9 phương pháp lặp tất cả, một số phương pháp hội tụ, còn một số thì cho kết quả khác lạ. Mô hình khối tôi test có số bậc tự do là 30000. (vì RAM trong máy tôi chỉ đủ dùng có thế, không sử dụng ổ cứng để so tốc độ tính toán).

Nếu so với tốc độ tính của cùng một số mô hình với phần tử khối trong chương trình tôi viết (FEModeler), sử dụng thử MtxVec 1.5.1 so với SAP2000 ver.9 và MIDAS/Civil ver. 6.3.0 thì thấy SAP2000 chạy nhanh gần gấp đôi so với FEModeler, MIDAS/Civil chạy nhanh gấp 3 so với FEModeler. Mô hình càng lớn thì sự chênh lệch về tốc độ giữa SAP2000 và FEModeler càng rút ngắn lại.

Cái MtxVec có nhược điểm: giải trực tiếp thì lại chỉ tính cho ma trận thưa bất kỳ, không tranh thủ được tính chất đối xứng, xác định dương của ma trận độ cứng kết cấu. Do vậy không hiệu quả trong lưu trữ và tính toán và chắc chắn không thể là số một ???!!!

Tôi sẽ tìm hiểu thêm về phương pháp giải lặp trong MtxVec, nghe tên đã thấy phức tạp rồi.

To Hien Nghiem: em cần tính với mô hình càng lớn càng tốt, vì định tập trung tính toán ứng suất cục bộ trong cầu. Thuật giải càng nhanh càng tốt. Có lẽ em sẽ viết lại Cholesky theo Sparse matrix xem. Bác hỗ trợ em phần này được không?
tnlinh vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Old 27-04-2006   #22
phu_ho
Thành viên nhiệt huyết
 
phu_ho's Avatar
 
Tham gia ngày: Aug 2004
Bài gởi: 716
Thanks: 3
Thanked 455 Times in 224 Posts
Default Ðề: Tính toán song song cho bài toán FEA ?

Trích:
Nguyên văn bởi hien nghiem
Bác Phu_ho:
Nếu giải 100000eqs mất <9s thì kiểu gì cũng là rất nhanh, tôi chưa được chứng kiến cái nào nhanh như thế.
Cái quan trọng là số nonzeros chứ số eqs có quan trọng mấy đâu bác ? (kiểu như mà bác mà có 1 cái ma trận chéo 100000 eqs thì có khi 1s đã xong rồi ấy chứ ).

À hôm trước tôi nhớ nhầm 1 số không nên mới chê cái việc giải ma trận 241mB (30 triệu nonzeros) của bác trong 9 phút là chậm. Nếu đúng thế thật thì cái solver của bác nhanh dã man rồi Mấy hôm vừa rồi tôi vừa tí toáy nghịch cái TAUCS đấy, chạy cũng tít phết, nhưng mà nếu NNZ tận 30 triệu thì chắc chắn phải lâu hơn 9 phút nhiều lần.
__________________
Does engineering need science?
phu_ho vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Old 27-04-2006   #23
hien nghiem
Thành viên nhiệt huyết
 
hien nghiem's Avatar
 
Tham gia ngày: Apr 2005
Nơi Cư Ngụ: Từ Liêm - Hà Nội
Bài gởi: 515
Thanks: 6
Thanked 369 Times in 123 Posts
Default Ðề: Tính toán song song cho bài toán FEA ?

Trích:
Nguyên văn bởi phu_ho
Cái quan trọng là số nonzeros chứ số eqs có quan trọng mấy đâu bác ? (kiểu như mà bác mà có 1 cái ma trận chéo 100000 eqs thì có khi 1s đã xong rồi ấy chứ ).

À hôm trước tôi nhớ nhầm 1 số không nên mới chê cái việc giải ma trận 241mB (30 triệu nonzeros) của bác trong 9 phút là chậm. Nếu đúng thế thật thì cái solver của bác nhanh dã man rồi Mấy hôm vừa rồi tôi vừa tí toáy nghịch cái TAUCS đấy, chạy cũng tít phết, nhưng mà nếu NNZ tận 30 triệu thì chắc chắn phải lâu hơn 9 phút nhiều lần.
Bác nói đúng, nó phụ thuộc vào nonzero cho nên khi đưa ra số liệu để so sánh, tôi đưa luôn cả kích thước ma trận để các bác tiện so sánh và đánh giá.
Cái 9 phút đấy là chỉ kể thời gian cholesky thôi (tôi đã ghi rõ ở trên, cái này là do tôi cho nó in-core và thừa RAM, nếu thiếu RAM (256MB) thì mất tới 40 phút). Thằng SAP2000 cũng giải bài này nhưng nó lưu trữ Skyline (912MB) mà Cholesky cũng chỉ trong 10 phút (RAM 1.2G). Tôi phục nó ở chỗ thuật toán đánh số nút của nó (không biết là Cullthin-McKee, hay Sloan's hay cái gì đó mà hiệu quả thế), nếu tôi dùng lưu trữ Skyline, ma trận độ cứng chiếm 2.3G (bản SAP7 mất 1.7G)
Cái nonzero 30tr này là đã kể đến fill-in rồi. Tôi không biết là cái món AMD này đã hiệu quả chưa, vì khi chưa fill-in, nonzero của nó cũng chỉ có hơn 3tr thôi (25MB). Tôi chưa thấy ứng ý lắm về cái này, không biết cái bác làm thế nào. Tôi cũng đang muốn post cái tôi làm lên đây mà chưa có thời gian để tạo thành một module giải riêng như kiểu UMFPACK hoặc TAUCS.
Bác phu_ho đã chạy thử cái Iterative solver của TAUCS chưa? Nếu chạy rồi mong bác cho anh em ít thông tin về kết quả với.
To Linh: Dạo này anh bận quá nên chưa giúp em được

thay đổi nội dung bởi: hien nghiem, 28-04-2006 lúc 08:18 AM
hien nghiem vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Old 01-05-2006   #24
phu_ho
Thành viên nhiệt huyết
 
phu_ho's Avatar
 
Tham gia ngày: Aug 2004
Bài gởi: 716
Thanks: 3
Thanked 455 Times in 224 Posts
Default Ðề: Tính toán song song cho bài toán FEA ?

Bác hien_nghiem: hóa ra 30 triệu là là fill-in, thế mà tôi cứ tưởng cái ma trận độ cứng Nếu bác có mấy cái ma trận to đấy thì bác dùng thử cái thư viện của Metis xem sao bác ạ. Tôi cũng chỉ mới thử dùng trong một vài trường hợp thôi, nhưng thấy rõ ràng với hệ to thì Metis hiệu quả hơn hẳn so với AMD. Ví dụ với 1 cái khối 32x32x32 phần tử khối 8 nút, số eqs khoảng hơn 100K, số nonzeros khoảng > 4 mils, available memory là 768 MB thì nếu ordering bằng AMD thì giải mất 1h trong khi dùng metis chỉ mất có 10 phút (size của fill-in matrix của AMD là 200 triệu và của Metis là 100 triệu). Tôi cũng không hiểu sao Metis nó làm khiếp thế, đúng là "mê tít" thật

Iterative của TAUCS thì tôi chưa thử bác ạ. Thực ra là cũng thử chạy vài lần những thấy vẫn lẩm cẩm lắm. Để lúc nào có thời gian tôi sẽ test kỹ hơn xem sao.
__________________
Does engineering need science?
phu_ho vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Old 01-05-2006   #25
hien nghiem
Thành viên nhiệt huyết
 
hien nghiem's Avatar
 
Tham gia ngày: Apr 2005
Nơi Cư Ngụ: Từ Liêm - Hà Nội
Bài gởi: 515
Thanks: 6
Thanked 369 Times in 123 Posts
Default Ðề: Tính toán song song cho bài toán FEA ?

Trích:
Nguyên văn bởi phu_ho
Bác hien_nghiem: hóa ra 30 triệu là là fill-in, thế mà tôi cứ tưởng cái ma trận độ cứng Nếu bác có mấy cái ma trận to đấy thì bác dùng thử cái thư viện của Metis xem sao bác ạ. Tôi cũng chỉ mới thử dùng trong một vài trường hợp thôi, nhưng thấy rõ ràng với hệ to thì Metis hiệu quả hơn hẳn so với AMD. Ví dụ với 1 cái khối 32x32x32 phần tử khối 8 nút, số eqs khoảng hơn 100K, số nonzeros khoảng > 4 mils, available memory là 768 MB thì nếu ordering bằng AMD thì giải mất 1h trong khi dùng metis chỉ mất có 10 phút (size của fill-in matrix của AMD là 200 triệu và của Metis là 100 triệu). Tôi cũng không hiểu sao Metis nó làm khiếp thế, đúng là "mê tít" thật

Iterative của TAUCS thì tôi chưa thử bác ạ. Thực ra là cũng thử chạy vài lần những thấy vẫn lẩm cẩm lắm. Để lúc nào có thời gian tôi sẽ test kỹ hơn xem sao.
Cảm ơn bác phu_ho. Tôi sẽ thử xem.
hien nghiem vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Old 19-05-2006   #26
phu_ho
Thành viên nhiệt huyết
 
phu_ho's Avatar
 
Tham gia ngày: Aug 2004
Bài gởi: 716
Thanks: 3
Thanked 455 Times in 224 Posts
Default Ðề: Tính toán song song cho bài toán FEA ?

Tôi vừa nhặt được cái bài này trong đó nó so sánh performance của vài cái iterative scheme với lại vài cái preconditioner khác nhau đấy bác hien_nghiem ạ. Bài toán 3D massif (linear elastic and elastoplastic). Không biết bác xem chưa thì ngó thử xem sao.

Mroueh, H. and I. Shahrour (1999). "Use of sparse iterative methods for the resolution of three-dimensional soil/structure interaction problems." International Journal for Numerical and Analytical Methods in Geomechanics 23(15): 1961 - 1975.

http://www3.interscience.wiley.com/c...01711/ABSTRACT
__________________
Does engineering need science?
phu_ho vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Old 20-05-2006   #27
hien nghiem
Thành viên nhiệt huyết
 
hien nghiem's Avatar
 
Tham gia ngày: Apr 2005
Nơi Cư Ngụ: Từ Liêm - Hà Nội
Bài gởi: 515
Thanks: 6
Thanked 369 Times in 123 Posts
Default Ðề: Tính toán song song cho bài toán FEA ?

Trích:
Nguyên văn bởi phu_ho
Tôi vừa nhặt được cái bài này trong đó nó so sánh performance của vài cái iterative scheme với lại vài cái preconditioner khác nhau đấy bác hien_nghiem ạ. Bài toán 3D massif (linear elastic and elastoplastic). Không biết bác xem chưa thì ngó thử xem sao.

Mroueh, H. and I. Shahrour (1999). "Use of sparse iterative methods for the resolution of three-dimensional soil/structure interaction problems." International Journal for Numerical and Analytical Methods in Geomechanics 23(15): 1961 - 1975.

http://www3.interscience.wiley.com/c...01711/ABSTRACT
Cảm ơn bác Phu_Ho. Tôi download rồi. Hy vọng sẽ giải quyết được nốt cái món Iterative solve này
hien nghiem vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Loading...
Trả lời

Ðiều Chỉnh
Xếp Bài

Quyền sử dụng ở Diễn Ðàn
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is Mở
Smilies đang Mở
[IMG] đang Mở
HTML đang Mở

Chuyển đến



Múi giờ GMT. Hiện tại là 12:18 AM.


DIỄN ĐÀN http://ketcau.com/forum NƠI HỘI TỤ CỦA CÁC KỸ SƯ KẾT CÂU VIỆT NAM
WWW.KETCAU.COM - CẦU NỐI CỦA CÁC KỸ SƯ KẾT CẤU CÔNG TRÌNH, ĐỊA KỸ THUẬT VIỆT NAM. DIỄN ĐÀN ĐƯỢC MUA BẢN QUYỀN CỦA JELSOFT ENTERPRISES Ltd. Thiết kế website