核心架构¶

PhysEngine/
├── Solver/              # 求解器模块
│   ├── SPHSolver.cc     # SPH流体求解器
│   ├── FESolver.cc      # FEM结构求解器
│   ├── FSISolver.cc     # 流固耦合求解器
│   └── Octree.cc        # 八叉树邻居搜索
├── Mesh/                # 网格模块
│   ├── MeshReader.cc    # 网格读取（ABAQUS .inp）
│   └── MeshData.cc      # 网格数据结构
├── Method/              # 有限元方法
│   ├── ShapeFunction.cc # 形状函数计算
│   ├── Integration.cc   # 数值积分规则
│   ├── H1R_Quad.cc      # 四边形单元
│   └── H1R_Hex.cc       # 六面体单元
├── LinearAlgebra/       # 线性代数库
│   ├── Vector.cc        # 向量运算
│   ├── DenseMatrix.cc   # 稠密矩阵
│   ├── SparseMat.cc     # 稀疏矩阵（CSR格式）
│   └── LUDecompose.cc   # LU分解求解器
└── General/             # 通用基础设施
    ├── Array.cc         # 数组操作
    ├── Memory.cc        # 内存管理
    └── Backends/        # 计算后端抽象
        ├── CPU.cc       # CPU后端
        ├── CUDA.cc      # CUDA GPU后端
        └── HIP.cc       # HIP GPU后端

graph TD
    A[输入网格文件] --> B[Mesh模块解析]
    B --> C[网格数据初始化]
    C --> D[物理参数设置]
    D --> E{求解类型}
    E -->|纯流体| F[SPHSolver]
    E -->|纯结构| G[FESolver]
    E -->|流固耦合| H[FSISolver]
    F --> I[八叉树搜索]
    G --> J[刚度矩阵组装]
    H --> K[界面数据交换]
    I --> L[时间步迭代]
    J --> L
    K --> L
    L --> M[结果输出]