Skip to content

GitLab

  • Menu
Projects Groups Snippets
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
  • Sign in
  • S src
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
    • Locked Files
  • Issues 59
    • Issues 59
    • List
    • Boards
    • Service Desk
    • Milestones
    • Iterations
  • Merge requests 4
    • Merge requests 4
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Schedules
  • Deployments
    • Deployments
    • Environments
    • Releases
  • Monitor
    • Monitor
    • Incidents
  • Analytics
    • Analytics
    • Value stream
    • CI/CD
    • Code review
    • Issue
    • Repository
  • Wiki
    • Wiki
  • Snippets
    • Snippets
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
Collapse sidebar
  • OPAL
  • src
  • Issues
  • #469

Closed
Open
Created Feb 18, 2020 by snuverink_j@snuverink_jDeveloper

asciih5block not working with recent H5Hut library?

Summary

The asciih5block script returns error messages and an empty (header only) file

Steps to reproduce

module use unstable; module load gcc/7.3.0 mpich/3.3 hdf5/1.10.4 H5hut/2.0.0rc5

I believe this will happen with any input file, an example is given here.

Example input file: NewChimney_R000_D000E.txt

Modify ascii2h5block_asgic.cpp for the right dimensions: ascii2h5block_asgic.cpp

$ export LDLIBS=-lH5hut
$ make ascii2h5block_asgic
$ ./ascii2h5block_asgic NewChimney_R000_D000E.txt out
--------------------------------------------------------
Using NewChimney_R000_D000E.txt to create test_CYC.h5part
Lines in finE: 67370
Hardcoded limits: x(-90/-10) cm
Hardcoded limits: y(-30/30) cm
Hardcoded limits: z(-1/1) cm
Hardcoded spacing: 0.5 cm
Hardcoded nlines: 67367
File nlines: 67365
Grid dimensions: Px = 23 , Py = 29 , Pz = 101
E_max = (29.0302, -5.6267, -23.2274) V/cm at index 20641.
Converting from V/cm and cm to kV/mm and mm before saving h5part
[proc 0] E: H5Block3dSetView: Iteration is invalid! Have you set the time step?
[proc 0] E: H5Block3dWriteVector3dFieldFloat64: No view has been defined!
[proc 0] E: H5Block3dWriteVector3dFieldFloat64: No view has been defined!
Done, bye ...

What is the current bug behavior?

Header only H5-file

What is the expected correct behavior?

H5-file with correct field

Possible fixes

The H5SetStep should be set before the View:

        H5SetStep(file, 0);
        h5_err_t h5err = H5Block3dSetView(file,
                                          0, gridPx - 1,
                                          0, gridPy - 1,
                                          0, gridPz - 1);

This modified ascii2h5block_asgic.cpp runs successfully.

Assignee
Assign to
Time tracking