jaabell / shakermaker Goto Github PK

Hi José,
I was trying to install shakermaker on my Mac (Ventura osx ) with Intel processors, but after many tries, I still get errors as below when I run 'python setup.py install'. I am wondering if you happen to come across this error in the setup stage. This is most likely mac-specific/clang/gcc related issue.

error.docx

Hello Dr. @jaabell ,

Recently, I installed the latest version of ShakerMaker, and during running examples script, the program shows the following error:

>> python3 example0_readme_example.py 

ShakerMaker Run begin. dt=0.05 nfft=4096 dk=0.3 tb=1000 tmin=0.0 tmax=100
-------------------------------------------------------------------------
 *** 0 psource.tt=0 t0[0]=-49.166668 dh=4.0 dz=4
Traceback (most recent call last):
  File "/usr/local/lib/python3.10/dist-packages/shakermaker/shakermaker.py", line 254, in run
    station.add_to_response(z_stf, e_stf, n_stf, t, tmin, tmax)
  File "/usr/local/lib/python3.10/dist-packages/shakermaker/station.py", line 75, in add_to_response
    self._z[nskip:(nskip+len(z))] = z
ValueError: could not broadcast input array from shape (8192,) into shape (983,)
0 of 1 done t[0]=-49.1667 t[-1]=360.3833 (tmin=0.0000 tmax=100.0000)

ShakerMaker Run done. Total time: 0.13463470100032282 s
------------------------------------------------

However the example scripts work when I go back to previous commit of station.py in 1a37a02 .

Could you please check this commit.

Regards,
Hossein Bazeghi

Hello
I would like to check whether the H5DRM is functioning properly in OpenSees. Specifically, I wanted to compare the response at the center point of the DRM box example provided in the OpenSees documentation between OpenSees and Shakermaker. However, while the acceleration results match well, the velocity and displacement results are different. It's surprising because the numerical integration of the accelerations yields the same results as OpenSees, but Shakermaker produces a completely different answer.

I also have another question regarding the output of station.getresponse() in Shakermaker. From what I have seen in the Shakermaker code, it appears that the output represents velocities. Is that correct?