zhangyiss/gctl_mesh
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

bug fixed

Browse Source
This commit is contained in:
张壹 2025-01-16 09:51:10 +08:00
parent cc2f1d49b5
commit bde2b5fa88
2 changed files with 5 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
lib/mesh/regular_grid.cpp
View File

@ -192,6 +192,7 @@ int gctl::regular_grid::view(std::string datname)
{ {
meshdata &data = get_data(datname); meshdata &data = get_data(datname);
mathgl_dens plt; ///< mathgl绘图对象
plt.range(rg_xmin, rg_xmin + (rg_xnum - 1)*rg_dx, rg_ymin, rg_ymin + (rg_ynum - 1)*rg_dy); plt.range(rg_xmin, rg_xmin + (rg_xnum - 1)*rg_dx, rg_ymin, rg_ymin + (rg_ynum - 1)*rg_dy);
plt.demension(rg_xnum, rg_ynum); plt.demension(rg_xnum, rg_ynum);
plt.add_dens(data.datval_, datname); plt.add_dens(data.datval_, datname);
@ -211,8 +212,9 @@ int gctl::regular_grid::view(std::string datname)
void gctl::regular_grid::plot(std::string datname) void gctl::regular_grid::plot(std::string datname)
{ {
meshdata data = get_data(datname); meshdata &data = get_data(datname);
gmt_JX_single pic;
pic.set_command("psconvert", "-A -TG -E300"); pic.set_command("psconvert", "-A -TG -E300");
pic.plot(datname, data.datval_, pic.plot(datname, data.datval_,
rg_xmin, rg_xmin + (rg_xnum - 1)*rg_dx, rg_xmin, rg_xmin + (rg_xnum - 1)*rg_dx,
@ -329,7 +331,7 @@ void gctl::regular_grid::load_netcdf_grid(std::string filename, mesh_data_type_e
new_data.create(NodeData, Scalar, node_num_, in_name[i], true, GCTL_BDL_MAX); new_data.create(NodeData, Scalar, node_num_, in_name[i], true, GCTL_BDL_MAX);
for (size_t j = 0; j < node_num_; j++) for (size_t j = 0; j < node_num_; j++)
{ {
new_data.datval_[i] = in_arr[i][j]; new_data.datval_[j] = in_arr[i][j];
} }
datalist_.push_back(new_data); datalist_.push_back(new_data);
@ -348,7 +350,7 @@ void gctl::regular_grid::load_netcdf_grid(std::string filename, mesh_data_type_e
new_data.create(ElemData, Scalar, ele_num_, in_name[i], true, GCTL_BDL_MAX); new_data.create(ElemData, Scalar, ele_num_, in_name[i], true, GCTL_BDL_MAX);
for (size_t j = 0; j < ele_num_; j++) for (size_t j = 0; j < ele_num_; j++)
{ {
new_data.datval_[i] = in_arr[i][j]; new_data.datval_[j] = in_arr[i][j];
} }
datalist_.push_back(new_data); datalist_.push_back(new_data);

8
lib/mesh/regular_grid.h
View File

@ -342,14 +342,6 @@ namespace gctl
double rg_xmin, rg_ymin; ///< 规则网格的x方向最小值和y方向最小值 double rg_xmin, rg_ymin; ///< 规则网格的x方向最小值和y方向最小值
double rg_dx, rg_dy; ///< 规则网格的x方向单元长度和y方向单元长度 double rg_dx, rg_dy; ///< 规则网格的x方向单元长度和y方向单元长度
#ifdef GCTL_GRAPHIC_MATHGL
mathgl_dens plt; ///< mathgl绘图对象
#endif // GCTL_GRAPHIC_MATHGL
#ifdef GCTL_GRAPHIC_GMT
gmt_JX_single pic; ///< gmt绘图对象
#endif // GCTL_GRAPHIC_GMT
array<vertex2dc> nodes; ///< 规则网格的节点数组 array<vertex2dc> nodes; ///< 规则网格的节点数组
array<rectangle2d> elements; ///< 规则网格的单元数组 array<rectangle2d> elements; ///< 规则网格的单元数组
}; };