“…For restricted systems, such as BECs in a harmonic trap, in a double or triple well, or in an optical lattice, some numerical or theoretical methods have been employed to solve the Gross-Pitaevskii equation (GPE), for example, two-or three-mode approximation, the symplectic method, and the particle swarm optimization (PSO) algorithm. [8][9][10][11][12][13][14] For non-restricted systems, a few methods have been proposed to deal with the nonlinear tunneling problem from various aspects, for example, involving the nonlinear effect only within the barrier, [15] or gradually increasing nonlinear interaction as BECs approaching an atomic quantum dot in a waveguide. [16] The dispersive wave shocks, formed in the upstream when BECs collide with a potential, were suggested to be studied by hydromechanics methods (Riemann invariants).…”