pygpc-polynomial-chaos / pygpc Goto Github PK

Hi there,

I'm new to this python-package for uncertainty and sensitivity tools and I'd love to use it for some research as UA&SA packages in python are quite scarce.
I'm wondering why almost all of the examples in the repo don't run. Either the script can't find a certain model or function (that it requires to run properly) or I get a back to back to back to back ... error of this kind (see below). I have python 3.6 and installed the package as stated here on the home page. At first sight, due to these issues, i'm not sure if the package is that user-friendly at the moment..

Traceback (most recent call last):
File "", line 1, in
File "C:\Program Files\Python 3.6\lib\multiprocessing\spawn.py", line 105, in spawn_main
exitcode = _main(fd)
File "C:\Program Files\Python 3.6\lib\multiprocessing\spawn.py", line 114, in _main
prepare(preparation_data)
File "C:\Program Files\Python 3.6\lib\multiprocessing\spawn.py", line 225, in prepare
_fixup_main_from_path(data['init_main_from_path'])
File "C:\Program Files\Python 3.6\lib\multiprocessing\spawn.py", line 277, in _fixup_main_from_path
run_name="mp_main")
File "C:\Program Files\Python 3.6\lib\runpy.py", line 263, in run_path
pkg_name=pkg_name, script_name=fname)
File "C:\Program Files\Python 3.6\lib\runpy.py", line 96, in _run_module_code
mod_name, mod_spec, pkg_name, script_name)
File "C:\Program Files\Python 3.6\lib\runpy.py", line 85, in _run_code
exec(code, run_globals)
File "C:\Users\myvhove\Documents\GitHub\pygpc\examples\gpc\plot_basis.py", line 45, in
interaction_order=3)
File "C:\Program Files\Python 3.6\lib\site-packages\pygpc\Basis.py", line 168, in init_basis_sgpc
with multiprocessing.Pool(multiprocessing.cpu_count()) as pool:
File "C:\Program Files\Python 3.6\lib\multiprocessing\context.py", line 119, in Pool
context=self.get_context())
File "C:\Program Files\Python 3.6\lib\multiprocessing\pool.py", line 174, in init
self._repopulate_pool()
File "C:\Program Files\Python 3.6\lib\multiprocessing\pool.py", line 239, in _repopulate_pool
w.start()
File "C:\Program Files\Python 3.6\lib\multiprocessing\process.py", line 105, in start
self._popen = self._Popen(self)
File "C:\Program Files\Python 3.6\lib\multiprocessing\context.py", line 322, in _Popen
return Popen(process_obj)
File "C:\Program Files\Python 3.6\lib\multiprocessing\popen_spawn_win32.py", line 33, in init
prep_data = spawn.get_preparation_data(process_obj._name)
File "C:\Program Files\Python 3.6\lib\multiprocessing\spawn.py", line 143, in get_preparation_data
_check_not_importing_main()
File "C:\Program Files\Python 3.6\lib\multiprocessing\spawn.py", line 136, in _check_not_importing_main
is not going to be frozen to produce an executable.''')
RuntimeError:
An attempt has been made to start a new process before the
current process has finished its bootstrapping phase.

    This probably means that you are not using fork to start your
    child processes and you have forgotten to use the proper idiom
    in the main module:

        if __name__ == '__main__':
            freeze_support()
            ...

    The "freeze_support()" line can be omitted if the program
    is not going to be frozen to produce an executable.

I have installed the module using 'pip install pygpc'. However when importing pygpc I have the following error:

File "C:\Users\Alez_\Anaconda3\lib\site-packages\spyder_kernels\customize\spydercustomize.py", line 108, in execfile
exec(compile(f.read(), filename, 'exec'), namespace)

File "C:/.../pygpcTest.py", line 3, in
import pygpc

File "C:..._init_.py", line 2, in
import pygpc.GPC

File "C:...\pygpc\GPC.py", line 20, in
from pygpc.calc_gpc_matrix_cpu import calc_gpc_matrix_cpu

File "init.pxd", line 918, in init pygpc.calc_gpc_matrix_cpu

ValueError: numpy.ufunc size changed, may indicate binary incompatibility. Expected 216 from C header, got 192 from PyObject

I am using python 3.7.1 in Windows 10 Home x64 with Spyder.

plot_gpc giving below error in case of multiple QoI's with 2 random variables

File C:\ProgramData\Anaconda3\envs\pygpcnew\lib\site-packages\pygpc\validation.py:553, in plot_gpc(session, coeffs, random_vars, coords, results, n_grid, output_idx, fn_out, camera_pos, zlim)
550 ax1 = fig.add_subplot(1, 2, 1, projection='3d')
551 im1 = ax1.plot_surface(x1_2d, x2_2d, np.reshape(y_gpc[:, i], (x[1].size, x[0].size), order='f'),
552 cmap="jet", alpha=0.75, linewidth=0, edgecolors=None)
--> 553 ax1.scatter(coords[:, idx_global[0]], coords[:, idx_global[1]], results,
554 'k', alpha=1, edgecolors='k', depthshade=False)
555 ax1.set_title(r'gPC approximation', fontsize=fs)
556 ax1.set_xlabel(r"%s" % random_vars[0], fontsize=fs)

ValueError: shape mismatch: objects cannot be broadcast to a single shape. Mismatch is between arg 0 with shape (151,) and arg 2 with shape (604,).

Thanks

Seshagiri

I am getting a runtime error when I run the plot_gpc_validation.py file in the example gallery . The error log is similar to the one shown in #3. I am using Python 3.7.6 and numpy 1.19.0