Unable to run demo files

[chethanram7009]

Dear developers,

Thank you very much for developing this excellent toolbox and making it open source.

I learned about PyHawk during EGU 2026 and became very interested in exploring it further.

I carefully followed the user manual and README instructions from GitHub, and the installation completed successfully. The toolbox is running properly overall, and I was able to execute demo_2.py without any errors.

However, when I tried to run demo_1.py following the provided instructions, I encountered an error. I suspected the issue might be related to the accelerometer configuration parameters. When I changed "isRequired" in AccelerometerConfig to false and set "Parameter_Number" to 0, the script ran successfully. But as per my understanding, enabling accelerometer parameters is very much necessary.

I encountered a similar issue with demo_3.py as well.

I would be very grateful if you could help me understand and resolve this issue.

Best regards,
Chethan V A

(py-hawk) F:\F_Drive_Old_Sys\1_PHD_WORK\PyHawk\PyHawk-main\PyHawk-main\demo>python demo_1.py
GRACE level 1b reading for SatID.A: 100%|███████████████████████████████████████████████████████████| 1/1 [00:00<00:00, 5.64it/s]
GRACE level 1b reading for SatID.B: 100%|███████████████████████████████████████████████████████████| 1/1 [00:00<00:00, 5.78it/s]
Orbit conversion from ITRS into GCRS: 100%|███████████████████████████████████████████████| 17280/17280 [00:04<00:00, 4171.60it/s]

===============================================
Fix Payload.SCA for SatID.A
0 gap!!

===============================================
Fix Payload.SCA for SatID.B
0 gap!!

===============================================
Fix Payload.ACC for SatID.A
0 gap!!

===============================================
Fix Payload.ACC for SatID.B
0 gap!!
Arc unpacking for ACC_A: 100%|█████████████████████████████████████████████████████████████████████| 4/4 [00:00<00:00, 166.82it/s]
Arc unpacking for ACC_B: 100%|████████████████████████████████████████████████████████████████████| 4/4 [00:00<00:00, 1000.43it/s]
Arc unpacking for SCA_A: 100%|█████████████████████████████████████████████████████████████████████| 4/4 [00:00<00:00, 666.64it/s]
Arc unpacking for SCA_B: 100%|█████████████████████████████████████████████████████████████████████| 4/4 [00:00<00:00, 800.06it/s]
Arc unpacking for GNV_A: 100%|████████████████████████████████████████████████████████████████████| 4/4 [00:00<00:00, 1795.51it/s]
Arc unpacking for GNV_B: 100%|████████████████████████████████████████████████████████████████████| 4/4 [00:00<00:00, 1658.98it/s]
Arc unpacking for KBR_A: 100%|████████████████████████████████████████████████████████████████████| 4/4 [00:00<00:00, 3986.98it/s]
multiprocessing.pool.RemoteTraceback:
"""
Traceback (most recent call last):
File "C:\Users\Satishkumar.conda\envs\py-hawk\lib\multiprocessing\pool.py", line 125, in worker
result = (True, func(*args, **kwds))
File "C:\Users\Satishkumar.conda\envs\py-hawk\lib\multiprocessing\pool.py", line 48, in mapstar
return list(map(*args))
File "F:\F_Drive_Old_Sys\1_PHD_WORK\PyHawk\PyHawk-main\PyHawk-main\demo\demo_1.py", line 175, in _adjustOrbit
AdjustOrbit(accConfig=[self.__acc_config_A, self.__acc_config_B], date_span=self.__date_span, sat=self.__sat).
File "F:\F_Drive_Old_Sys\1_PHD_WORK\PyHawk\PyHawk-main\PyHawk-main\demo..\src\Solver\AdjustOrbit.py", line 426, in calibrate
self._ode_res = ode.propagate()
File "F:\F_Drive_Old_Sys\1_PHD_WORK\PyHawk\PyHawk-main\PyHawk-main\demo..\src\ODE\Instance\GravimetryODE.py", line 103, in propagate
propagator = ComplexODE(config).setSecondDerivative(self.__secDer, ArcLen=self._ArcLen - 4).
File "F:\F_Drive_Old_Sys\1_PHD_WORK\PyHawk\PyHawk-main\PyHawk-main\demo..\src\ODE\Common\ComplexODE.py", line 114, in setInitial
x1, r1, v1 = self.__SingleStepOde.propagate(int(self.__order / 2), -self.__stepSize)
File "F:\F_Drive_Old_Sys\1_PHD_WORK\PyHawk\PyHawk-main\PyHawk-main\demo..\src\ODE\Common\SingleStepODE.py", line 116, in propagate
x_0, r_0, v_0 = callfunc(self._func, x_0, r_0, v_0, stepsize)
File "F:\F_Drive_Old_Sys\1_PHD_WORK\PyHawk\PyHawk-main\PyHawk-main\demo..\src\ODE\Common\SingleStepODE.py", line 236, in _callFunc
k_1 = func(x_0, r_0, v_0)
File "F:\F_Drive_Old_Sys\1_PHD_WORK\PyHawk\PyHawk-main\PyHawk-main\demo..\src\SecDerivative\Instance\DualOrbitAndVar.py", line 159, in secDerivative
diff2ndA = MatrixVarEqn.VarEq2ndOrder(PhirSr[0:3], PhivSv[0:3], arA, apA, accA,
File "F:\F_Drive_Old_Sys\1_PHD_WORK\PyHawk\PyHawk-main\PyHawk-main\demo..\src\SecDerivative\MatrixVarEqn.py", line 163, in VarEq2ndOrder
Srp = np.matmul(ar, Sr) + ap
ValueError: operands could not be broadcast together with shapes (3,21) (3,6)
"""

The above exception was the direct cause of the following exception:

Traceback (most recent call last):
File "demo_1.py", line 197, in
demo1.run()
File "demo_1.py", line 129, in run
pool.map(self._adjustOrbit, self.__arclist)
File "C:\Users\Satishkumar.conda\envs\py-hawk\lib\multiprocessing\pool.py", line 364, in map
return self._map_async(func, iterable, mapstar, chunksize).get()
File "C:\Users\Satishkumar.conda\envs\py-hawk\lib\multiprocessing\pool.py", line 771, in get
raise self._value
ValueError: operands could not be broadcast together with shapes (3,21) (3,6)

[HawkFan]

Dear Chethan,

Thanks very much for your interest in our toolbox! Apologies for the late response regarding your issue. Indeed we do not regularly check the information here, but instead we have a forum to gather all reported issues together for the open discussion. Any post in that discussion forum would notify us immediately so that we can see and respond as soon as possible. Could you please kindly copy-paste your post to our discussion forum (https://github.com/orgs/NCSGgroup/discussions), so that others who may have the mutual issues can share.

Regarding your issue, the main developer, @PyHawkWu, is already notified and will come back to you soon. Thanks for your attention!

Fan Yang