tecio-python 0.3.2__tar.gz → 0.4.0__tar.gz

{tecio_python-0.3.2/tecio_python.egg-info → tecio_python-0.4.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: tecio-python
-Version: 0.3.2
+Version: 0.4.0
 Summary: Python interface for reading and writing Tecplot data files
 Project-URL: Homepage, https://github.com/meersman/tecio
 Project-URL: Documentation, https://meersman.github.io/tecio/

tecio_python-0.4.0/plugins/matlab/contourfe.m

@@ -0,0 +1,111 @@
+function f = contourfe(dim, x, y, z, c, node_map, varargin)
+% Filled patch contour for unstructured 2D/3D surface zones
+%
+% Call:
+%     h = contourfe(dim, x, y, z, c, node_map, ...)
+%
+% Args:
+%     dim: 2 or 3
+%     x,y[,z]: nodal coordinate column vectors (Nx1)
+%     c: contour values; either nodal (Nx1) or cell-centered (Mcells x 1)
+%     node_map: M x K integer matrix of node indices per cell (K=3 or 4)
+%     varargin: name-value pairs forwarded to patch (e.g. 'EdgeColor','none')
+%
+% Returns:
+%     f: handle to the patch object
+%
+% Notes:
+%  - If c is nodal, uses 'FaceVertexCData' with 'FaceColor'='interp'
+%  - If c is cell-centered, uses 'CData' per-face with 'FaceColor'='flat'
+%  - For dim==3, z must be provided and plotting uses 3D patch (viewable with lighting)
+%  - node_map may contain NaNs for unequal-sided cells; those faces are handled
+
+% Validate dim
+if ~ismember(dim, [2,3])
+    error('dim must be 2 or 3');
+end
+
+% Ensure column vectors
+x = x(:); y = y(:);
+if dim == 3
+    if nargin < 4 || isempty(z)
+        error('z must be provided for dim==3');
+    end
+    z = z(:);
+end
+
+% Determine number of nodes and faces
+nNodes = numel(x);
+[nFaces, ~] = size(node_map);
+
+% Validate node_map indices
+if any(node_map(:) > nNodes) || any(node_map(:) < 0 & ~isnan(node_map(:)))
+    error('node_map contains invalid node indices.');
+end
+
+% Prepare vertices matrix
+if dim == 2
+    verts = [x, y];
+else
+    verts = [x, y, z];
+end
+
+% Determine if c is nodal or cell-centered
+c = c(:);
+isNodal = (numel(c) == nNodes);
+isCell   = (numel(c) == nFaces);
+if ~(isNodal || isCell)
+    error('Length of c must match number of nodes (nodal) or number of faces (cell-centered).');
+end
+
+% Clean node_map: replace zeros -> NaN (if any), ensure double
+node_map = double(node_map);
+node_map(node_map==0) = NaN;
+
+% If node_map contains NaNs for some cells, patch accepts NaN-separated faces.
+% Create faces in appropriate format: patch accepts face vertex indices matrix.
+faces = node_map;
+
+% Build patch arguments
+pArgs = varargin;
+
+if isNodal
+    % Nodal data: use FaceVertexCData with interpolation
+    % Create patch with vertices and faces, supply FaceVertexCData
+    % Set FaceColor to 'interp' and EdgeColor as provided or 'none' by default
+    if ~any(strcmpi('FaceColor',pArgs))
+        pArgs = [{'FaceColor','interp'}, pArgs];
+    end
+    if ~any(strcmpi('EdgeColor',pArgs))
+        pArgs = [{'EdgeColor','none'}, pArgs];
+    end
+    % Create patch
+    f = patch('Vertices', verts, 'Faces', faces, 'FaceVertexCData', c, pArgs{:});
+
+else
+    % Cell-centered: supply per-face CData and use flat coloring
+    if ~any(strcmpi('FaceColor',pArgs))
+        pArgs = [{'FaceColor','flat'}, pArgs];
+    end
+    if ~any(strcmpi('EdgeColor',pArgs))
+        pArgs = [{'EdgeColor','none'}, pArgs];
+    end
+    % patch expects CData as Mx1 or Mx3 color; provide as Mx1
+    faceCData = c;
+    % Create patch
+    f = patch('Vertices', verts, 'Faces', faces, 'CData', faceCData, pArgs{:});
+
+end
+
+% Set colormap and colorbar behavior consistent with contourf
+if ~any(strcmpi('EdgeColor',pArgs))
+    set(f,'EdgeColor','none');
+end
+axis equal
+if dim == 3
+    view(3)
+    camlight headlight
+    lighting gouraud
+end
+colorbar
+end

tecio_python-0.4.0/plugins/matlab/tecmat2struct.m

@@ -0,0 +1,137 @@
+function zones = tecmat2struct(d)
+% Reshape a loaded tec2mat .mat struct into a struct array of zones with
+% dataset variables accessible by name.
+%
+% tec2mat outputs (info, zone_1, zone_2, ...) and variables by position
+% (var_1, var_2, ...) because valid Tecplot variable names can include
+% characters that are not valid MATLAB identifiers (e.g. "X [ft]"). When the
+% variable names are known to be valid identifiers, this function lets you
+% reformat the structures such that variables can be accessed by name, e.g.
+% plot(zones(1).x, zones(1).pressure) instead of plot(zone_1.var_1,
+% zone_1.var_3).
+%
+% Call:
+%     zones = tecmat2struct(d)
+%
+% Args:
+%     d: struct returned by load() on a tec2mat .mat file, i.e. the struct
+%        has fields info, zone_1, zone_2, ...
+%
+% Returns:
+%     zones: 1 x info.num_zones struct array. Every element has the same
+%            fields, so this is a normal MATLAB struct array (zones(i),
+%            [zones.pressure], etc. all work):
+%       .title, .zone_type, .solution_time, .strand_id
+%       .I, .J, .K                     ordered zones; [] for FE zones
+%       .num_nodes, .num_elements      FE zones; [] for ordered zones
+%       .node_map                      FE zones (resolved through
+%                                       node_map_shared_from if the zone
+%                                       shares connectivity); [] otherwise
+%       .face_neighbor_mode, .num_face_connections, .face_connections
+%                                       FE zones with face-neighbor data;
+%                                       '' / [] / [] otherwise
+%       .<variable name>                one field per dataset variable
+%                                       whose real name is a valid MATLAB
+%                                       identifier (see isvarname). Passive
+%                                       in a zone -> []. Shared -> resolved
+%                                       to the source zone's data, so every
+%                                       zones(i) is self-contained.
+%
+% Notes:
+%  - A variable name that is not a valid MATLAB identifier is skipped with a
+%    warning; it remains accessible the original way, via zone_<n>.var_<k> and
+%    info.var_names{k}.
+%  - Shared variables/connectivity are branched (copied from referenced zone's
+%    data into every zone that shares it). This can create large memory
+%    usage for files whose zones heavily share a static grid across many time
+%    steps.
+%
+% Example:
+%     d = load('flow.mat');
+%     zones = tecmat2struct(d);
+%     plot(zones(1).x, zones(1).pressure);
+
+if ~isfield(d, 'info')
+    error('tecmat2struct:badInput', ...
+        'Input does not look like a tec2mat output struct (no ''info'' field).');
+end
+
+numZones = double(d.info.num_zones);
+numVars = double(d.info.num_vars);
+varNames = d.info.var_names;
+
+% Reserved field names already used for zone metadata (variable names that
+% collides with a reserved metadat name is skipped).
+reserved = {'title', 'zone_type', 'solution_time', 'strand_id', ...
+    'I', 'J', 'K', 'num_nodes', 'num_elements', 'node_map', ...
+    'face_neighbor_mode', 'num_face_connections', 'face_connections'};
+
+% Validate variable names up once so every zone gets the exact same set of
+% variable fields in the same order.
+validName = false(1, numVars);
+for vi = 1:numVars
+    name = strtrim(char(varNames{vi}));
+    if ~isvarname(name)
+        warning('tecmat2struct:invalidName', ...
+            ['Variable %d (''%s'') is not a valid MATLAB identifier; skipping. ' ...
+             'Access it via zone_<n>.var_%d instead.'], vi, name, vi);
+    elseif any(strcmp(name, reserved))
+        warning('tecmat2struct:reservedName', ...
+            ['Variable %d (''%s'') collides with a zone metadata field; skipping. ' ...
+             'Access it via zone_<n>.var_%d instead.'], vi, name, vi);
+    else
+        validName(vi) = true;
+    end
+end
+varIdx = find(validName);
+
+for zi = numZones:-1:1  % descending so the first iteration preallocates zones(1)
+    zsrc = d.(sprintf('zone_%d', zi));
+
+    out = struct( ...
+        'title', zsrc.title, ...
+        'zone_type', zsrc.zone_type, ...
+        'solution_time', zsrc.solution_time, ...
+        'strand_id', zsrc.strand_id, ...
+        'I', [], 'J', [], 'K', [], ...
+        'num_nodes', [], 'num_elements', [], 'node_map', [], ...
+        'face_neighbor_mode', '', 'num_face_connections', [], ...
+        'face_connections', []);
+
+    if isfield(zsrc, 'I')
+        out.I = zsrc.I;
+        out.J = zsrc.J;
+        out.K = zsrc.K;
+    end
+    if isfield(zsrc, 'num_nodes')
+        out.num_nodes = zsrc.num_nodes;
+        out.num_elements = zsrc.num_elements;
+    end
+    if isfield(zsrc, 'node_map')
+        out.node_map = zsrc.node_map;
+    elseif isfield(zsrc, 'node_map_shared_from')
+        srcZone = d.(sprintf('zone_%d', zsrc.node_map_shared_from));
+        out.node_map = srcZone.node_map;
+    end
+    if isfield(zsrc, 'face_neighbor_mode')
+        out.face_neighbor_mode = zsrc.face_neighbor_mode;
+        out.num_face_connections = zsrc.num_face_connections;
+        out.face_connections = zsrc.face_connections;
+    end
+
+    for vi = varIdx
+        name = strtrim(char(varNames{vi}));
+        status = zsrc.var_status{vi};
+        if strcmp(status, 'shared')
+            srcIdx = zsrc.var_shared_from(vi);
+            srcZone = d.(sprintf('zone_%d', srcIdx));
+            out.(name) = srcZone.(sprintf('var_%d', vi));
+        else
+            % 'active' carries real data; 'passive' is already [].
+            out.(name) = zsrc.(sprintf('var_%d', vi));
+        end
+    end
+
+    zones(zi) = out;
+end
+end

{tecio_python-0.3.2 → tecio_python-0.4.0}/plugins/paraview/TecplotTecioReader.py

@@ -170,6 +170,9 @@ _AXIS_SYNONYMS: dict[str, frozenset[str]] = {
         "coordx",
         "xcoord",
         "x_coord",
+        "xgrid",
+        "x-grid",
+        "x_grid",
     }),
     "y": frozenset({
         "y",
@@ -180,6 +183,9 @@ _AXIS_SYNONYMS: dict[str, frozenset[str]] = {
         "coordy",
         "ycoord",
         "y_coord",
+        "ygrid",
+        "y-grid",
+        "y_grid",
     }),
     "z": frozenset({
         "z",
@@ -190,6 +196,9 @@ _AXIS_SYNONYMS: dict[str, frozenset[str]] = {
         "coordz",
         "zcoord",
         "z_coord",
+        "zgrid",
+        "z-grid",
+        "z_grid",
     }),
 }