Hình dung của lưới tứ giác & octree

Cái gọi là lưới quadtree và octree khá hấp dẫn cho các ứng dụng đòi hỏi phải tinh chỉnh lưới thích ứng. Chúng là ví dụ được sử dụng trong Gerris và Paramesh . Có ai biết một định dạng tệp tốt cho các lưới như vậy và hỗ trợ phần mềm trực quan không? Xem thêm câu trả lời này trong Câu hỏi thường gặp về Gerris. Ứng cử viên tiềm năng duy nhất mà tôi hiện đang biết đến là lớp VTK HyperOctree , nhưng dường như gần như không có tài liệu nào.

Như một ví dụ cụ thể về lưới, hãy xem xét lưới octree được tinh chỉnh thưa thớt, bao gồm hộp, mỗi hộp chứa ô. Các chiến lược hiện tại của tôi trong Visit / Paraview là:$10^4$$8 \times 8 \times 8=512$

1. Viết lưới dưới dạng tệp VTK không cấu trúc. Điều này rất tốn kém cho các bộ dữ liệu lớn hơn và không khai thác 'cấu trúc' của quad / octree.
2. Thu thập các 'hộp' lân cận trong cây thành các khối lớn hơn và viết lưới có cấu trúc khối. Điều này liên quan đến khá nhiều mã bổ sung và mất thông tin kết nối được nhúng trong quad / octree.

Cập nhật Nếu hiện tại không có định dạng khả thi để viết và hiển thị trực tiếp lưới octree, tôi sẽ đánh giá cao các đề xuất cho định dạng tệp có thể được sử dụng để viết octree dưới dạng lưới có cấu trúc khối. Lý tưởng nhất là định dạng này:

1. Có một định nghĩa / triển khai tương đối đơn giản, nguồn mở lý tưởng và có khả năng tương thích C / Fortran.
2. Có thể xử lý hiệu quả một số lượng lớn các khối lưới (ví dụ 1000 trở lên).
3. Cho phép dễ dàng chỉ định kết nối của các khối lưới, đặc biệt là tại các ranh giới sàng lọc.

Tôi hiện đang sử dụng Silo cho mục đích này, nhưng tôi đang tìm kiếm thứ gì đó đạt điểm cao hơn ở ba điểm trên (nhỏ hơn, nhanh hơn, đơn giản hơn).

Trong trường hợp của tôi khi tôi đang thực hiện sàng lọc đa cấp, tôi thường viết các tệp lưới ở các cấp khác nhau. Thông thường mức 0 là một lưới thô cơ sở vẫn cố định. Và đối với cấp 1 trở lên, tôi đã sử dụng các tệp khác nhau cho các bản vá khác nhau. Giả sử nếu bạn có 5 bản vá tốt ở cấp 1, chiến lược của tôi là viết năm tệp khác nhau vào cuối chương trình của tôi, đại loại như Level_1_Patch_1, Level_1_Patch_2 ... có thể được thực hiện một cách linh hoạt trong C (không bao giờ hoạt động ở Fortran).

        char filename1[64];
        sprintf(filename1, "Patch%d.dat", patch_number);

        FILE *file;
        file= fopen(filename1, "wb");

Tôi nên đề cập ở đây rằng tôi đã sử dụng một sự khác biệt hữu hạn, hệ thống lưới có cấu trúc.

Vì vậy, một khi tôi đã có tất cả các cấp độ và các tệp .PLT hoặc .DAT của chúng. Tôi nhập tất cả chúng trong TecPLOT. Ví dụ, khi tôi đang viết mã hai cấp của mình với một vài bản vá tốt ở cấp 1, lưới của bạn sẽ trông như thế này trong Tecplot.

Tôi không quen lắm với việc xử lý các lưới không có cấu trúc nhưng các lưới có cấu trúc được xử lý tốt trong Tecplot