Pith. sign in

REVIEW 2 major objections 5 minor 41 references

A new multi-component short-pulse equation admits Pfaffian N-solitons and an integrable fully discrete self-adaptive mesh scheme that matches exact solutions.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-11 14:00 UTC pith:OLDX54KX

load-bearing objection Clean Pfaffian full discretization of multi-component short-pulse equations that actually yields a usable self-adaptive mesh scheme, with solid algebra and small soliton errors. the 2 major comments →

arxiv 2607.04756 v1 pith:OLDX54KX submitted 2026-07-06 nlin.SI math-phmath.MP

Integrable full discretization of the multi-component short pulse equation

Ayako Hori , Yuta Tanaka , Ken-ichi Maruno , Yasuhiro Ohta , Bao-Feng Feng This is my paper
classification nlin.SI math-phmath.MP MSC 37K1037K4035Q5165M06
keywords integrable full discretizationmulti-component short pulse equationcoupled complex short pulse equationself-adaptive moving mesh schemePfaffianHirota bilinear methodN-soliton solutions
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a multi-component short-pulse (MCSP) equation that contains the coupled complex short pulse equation as a special case. It derives a bilinear form, constructs exact N-soliton solutions as Pfaffians, and then builds integrable semi-discrete and fully discrete versions that keep the same soliton solutions. The fully discrete system is used as a practical moving-mesh numerical scheme whose mesh intervals automatically tighten where the solution changes rapidly. For the one- and two-soliton cases tested, the numerical profiles stay within relative peak errors of order 10^{-3} or smaller of the exact continuous solutions. The work therefore supplies both an analytic multi-component model and a ready-to-use integrable discretization for accurate simulation of short-pulse dynamics.

Core claim

There exists a formulation of the multi-component short-pulse equation, admitting Pfaffian N-soliton solutions and including the coupled complex short-pulse equation by reduction, that possesses integrable semi-discrete and fully discrete analogues; the fully discrete system supplies a self-adaptive moving-mesh scheme whose numerical solutions agree with the exact continuous solitons for the parameter sets examined.

What carries the argument

Hirota bilinearization of the MCSP equation together with Pfaffian solutions of the resulting bilinear system; the same Pfaffians, after discrete shifts of the spectral factors, solve the semi-discrete and fully discrete bilinear equations and generate the self-adaptive mesh via a discrete hodograph transformation.

Load-bearing premise

That continuous exact solutions sampled at two successive time levels, together with the discrete mesh update started from the continuous hodograph mesh, keep every denominator nonzero and that the observed accuracy near soliton peaks continues to hold beyond the specific one- and two-soliton parameter sets that were tested.

What would settle it

Run the fully discrete scheme on a multi-soliton or non-soliton initial condition for which an independent high-accuracy continuous solution is known; if the numerical peak relative error grows systematically above the reported 10^{-3} level or if a denominator vanishes, the claim of a practical high-accuracy integrable scheme fails.

Watch this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 5 minor

Summary. The paper introduces a multi-component short-pulse (MCSP) system (1.2) that reduces to the coupled complex short-pulse (CCSP) equation under conjugacy. Via a hodograph map it is bilinearized to (2.12); N-soliton solutions are given in Pfaffian form (Theorem 2.1, Lemma 2.2). Integrable semi-discrete (3.1)/(3.12) and fully discrete (4.1)/(4.30) analogues are constructed that admit the same Pfaffian solutions (Theorems 4.1–4.3). The fully discrete system is proposed as a self-adaptive moving-mesh scheme; Section 5 reports numerical tests against continuous one- and two-soliton solutions for 2-SP, 1-CSP and 2-CSP reductions, with peak relative errors of order 10^{-3}–10^{-4}.

Significance. An integrable full discretization of multi-component short-pulse equations that remains compact enough for practical numerics has been missing; earlier Wronskian-based full discretizations of the scalar SP equation were too cumbersome to extend. The Pfaffian construction yields a clean bilinear system, explicit N-soliton formulae, and a self-adaptive mesh whose intervals automatically concentrate where the solution varies rapidly. Concrete maxerr values near soliton peaks, together with continuum limits recovering the semi-discrete and continuous equations (Remark 4.4), make the work a solid contribution to discrete integrable systems and to numerical methods for ultrashort-pulse models.

major comments (2)
  1. Section 5 and Theorem 4.3: the fully discrete scheme (4.30) and the bidiagonal solves are well-defined only when all displayed denominators remain nonzero. The manuscript assumes this (and 1-ab eq0) but never verifies it for the initialized data of Section 5, nor supplies a practical safeguard (e.g., a check or a fallback when a denominator vanishes). A short statement or numerical monitor confirming that the denominators stayed safely away from zero throughout the reported runs would remove the remaining load-bearing caveat.
  2. Section 5: all accuracy claims rest on continuous exact one- and two-soliton initial data (and the corresponding continuous hodograph mesh). While the reported maxerr values are convincing for those families, the paper does not test non-soliton or multi-soliton (N>2) data, nor does it examine long-time stability under the adaptive mesh. A brief additional experiment or an explicit caveat that the observed accuracy is demonstrated only for the soliton families of Section 5 would keep the numerical claim proportionate to the evidence.
minor comments (5)
  1. Abstract and Introduction: the phrase “practical self-adaptive moving mesh scheme” is repeated; a single, precise statement would suffice.
  2. Equation (2.17) and subsequent Pfaffian elements: the condition that all denominators are nonzero is stated only once; a short global remark at the beginning of Section 2 would avoid later repetition.
  3. Figures 5.1–5.22: the captions report maxerr but do not indicate the final time T or the number of time steps; adding these data would make the error figures self-contained.
  4. References [33] and [37–38] are earlier works by overlapping authors on semi-discretizations; a one-sentence clarification of what is new relative to those papers would help the reader.
  5. Typographical: “sho rt” appears in the running title; “T wo-soliton” and similar line-break artifacts occur in several places (e.g., Examples 2.3, 4.5).

Circularity Check

0 steps flagged

No significant circularity: discrete equations and Pfaffian solutions are constructed independently; numerics only verify consistency.

full rationale

The paper's central claims rest on an explicit bilinear construction. Continuous MCSP is mapped via hodograph (2.3) and dependent-variable change (2.11) to bilinear system (2.12); Theorem 2.1 and Lemma 2.2 then exhibit Pfaffian N-soliton solutions by direct expansion and Pfaffian identities. Semi-discrete (3.1) and fully discrete (4.1) bilinear systems are postulated by the same method; Theorems 4.1–4.3 and Lemmas 3.1–3.2 verify that the same Pfaffian elements (with discrete exponential factors) satisfy them and that the dependent-variable map (4.25) produces the nonlinear difference equations (4.30). Continuous and semi-discrete limits are recovered by Taylor expansion (Remark 4.4). Numerical Section 5 initializes the scheme from the continuous exact solutions and reports peak relative errors; this is consistency checking, not a fitted prediction. Self-citations to earlier SP/CSP discretizations by overlapping authors supply background technique (hodograph + bilinear) but are not load-bearing for the new multi-component full discretization or its Pfaffian solutions. No parameter is fitted to data and then re-presented as a prediction; no uniqueness theorem is imported solely by self-citation to forbid alternatives. The derivation is therefore self-contained against its own algebraic benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 3 axioms · 1 invented entities

The work rests on standard Hirota bilinear calculus, known Pfaffian identities, and the classical continuous MCSP/CCSP equations. No free parameters are fitted to data; lattice spacings and soliton parameters are free choices of the numerical experiments, not of the analytic claims. The only invented objects are the new bilinear and discrete systems themselves, which are defined by explicit equations and therefore carry independent algebraic content.

free parameters (2)
  • lattice spacings a, b = a=0.005, b=0.0005
    Chosen by hand for the numerical tests (a=0.005, b=0.0005); they do not enter the analytic integrability statements.
  • soliton spectral and phase parameters (p_j, B_j, a_j^{(μ)}, b_k^{(ν)})
    Free constants that label particular exact solutions used for numerical illustration; not fitted to external data.
axioms (3)
  • standard math Hirota D-operator bilinear calculus and the standard Pfaffian expansion identities (cited as [41]) hold for the introduced elements.
    Used throughout Sections 2–4 to verify the bilinear equations.
  • domain assumption The continuous multi-component short-pulse equation (1.2) and its complex reduction (1.3) are the correct continuum models.
    Taken as given from the literature; the paper’s contribution is their discretization.
  • ad hoc to paper All denominators appearing in the discrete equations and Pfaffian elements remain nonzero for the solutions under consideration.
    Stated as a standing assumption in Section 4; required for the dependent-variable transformations and the numerical scheme to be well-defined.
invented entities (1)
  • Pfaffian multi-component short-pulse bilinear system (2.12) and its fully discrete analogue (4.1)/(4.30) independent evidence
    purpose: Provide a compact integrable formulation that admits Pfaffian N-solitons and yields a practical self-adaptive mesh scheme.
    Defined by explicit equations; independent algebraic content is the existence of the Pfaffian solutions and the continuum limits.

pith-pipeline@v1.1.0-grok45 · 43529 in / 2610 out tokens · 20974 ms · 2026-07-11T14:00:39.931741+00:00 · methodology

0 comments
read the original abstract

We propose a new formulation of the multi-component short pulse (MCSP) equation that includes the coupled complex short pulse (CCSP) equation as a reduction. Using Hirota's bilinear method, we construct its $N$-soliton solutions in Pfaffian form. We then derive integrable semi-discrete and fully discrete analogues of the MCSP equation admitting Pfaffian $N$-soliton solutions. The resulting fully discrete system provides a practical self-adaptive moving mesh scheme for numerical simulations. For the parameter sets considered, numerical simulations demonstrate excellent agreement between the numerical and exact solutions, confirming the robustness and high accuracy of the proposed scheme.

Figures

Figures reproduced from arXiv: 2607.04756 by Ayako Hori, Bao-feng Feng, Ken-ichi Maruno, Yasuhiro Ohta, Yuta Tanaka.

Figure 5.1
Figure 5.1. Figure 5.1: Numerical simulation of the u (1)-profile of the one-soliton solution of the 2-SP equation. maxerr(u (1)) = 9.33 × 10−4 [PITH_FULL_IMAGE:figures/full_fig_p033_5_1.png] view at source ↗
Figure 5.2
Figure 5.2. Figure 5.2: Numerical simulation of the v (1)-profile of the one-soliton solution of the 2-SP equation. maxerr(v (1)) = 9.80 × 10−4 [PITH_FULL_IMAGE:figures/full_fig_p033_5_2.png] view at source ↗
Figure 5.3
Figure 5.3. Figure 5.3: Numerical simulation of the u (1)-profile of the two-soliton solution of the 2-SP equation. maxerr(u (1)) = 1.00 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p034_5_3.png] view at source ↗
Figure 5.4
Figure 5.4. Figure 5.4: Numerical simulation of the v (1)-profile of the two-soliton solution of the 2-SP equation. maxerr(v (1)) = 3.70 × 10−4 [PITH_FULL_IMAGE:figures/full_fig_p034_5_4.png] view at source ↗
Figure 5.5
Figure 5.5. Figure 5.5: Numerical simulation of the |u (1)|-profile of the one-soliton solution of the 1-CSP equation. maxerr(|u (1)|) = 1.93 × 10−4 . One-soliton: Applying the constraints p1 = p ∗ 2 , B1 = B∗ 2 , and a (1) 1 = (b (1) 2 ) ∗ to (5.13) gives u (1) = g (1) f , x = X − 2(log f)T , t = T, f = 1 + 1 4  p1p ∗ 1 p1 + p ∗ 1 2 ϕ1ϕ ∗ 1a (1) 1 (a (1) 1 ) ∗ , g(1) = a (1) 1 ϕ1, (5.17) where ϕ1 = B1e p1X+p −1 1 T [PITH_F… view at source ↗
Figure 5.6
Figure 5.6. Figure 5.6: Numerical simulation of the Re(u (1))-profile of the one-soliton solution of the 1-CSP equation. maxerr(Re(u (1))) = 2.99 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p036_5_6.png] view at source ↗
Figure 5.7
Figure 5.7. Figure 5.7: Numerical simulation of the Im(u (1))-profile of the one-soliton solution of the 1-CSP equation. maxerr(Im(u (1))) = 4.32 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p036_5_7.png] view at source ↗
Figure 5.8
Figure 5.8. Figure 5.8: Numerical simulation of the |u (1)|-profile of the two-soliton solution of the 1-CSP equation. maxerr(|u (1)|) = 4.03 × 10−4 . Figures 5.5–5.7 show the numerical |u (1)|, Re(u (1)), and Im(u (1)) profiles for the one￾soliton solution with p1 = 0.5 + 0.1i, a (1) 1 = 1 + 0.5i, and B1 = exp(−10). Figures 5.8–5.10 show the corresponding profiles for the two-soliton solution with p1 = 0.5+0.1i, p2 = 2−0.5i, a… view at source ↗
Figure 5.9
Figure 5.9. Figure 5.9: Numerical simulation of the Re(u (1))-profile of the two-soliton solution of the 1-CSP equation. maxerr(Re(u (1))) = 2.07 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p038_5_9.png] view at source ↗
Figure 5.10
Figure 5.10. Figure 5.10: Numerical simulation of the Im(u (1))-profile of the two-soliton solution of the 1-CSP equation. maxerr(Im(u (1))) = 3.87 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p038_5_10.png] view at source ↗
Figure 5.11
Figure 5.11. Figure 5.11: Numerical simulation of the |u (1)|-profile of the one-soliton solution of the 2-CSP equation. maxerr(|u (1)|) = 1.67 × 10−4 [PITH_FULL_IMAGE:figures/full_fig_p040_5_11.png] view at source ↗
Figure 5.12
Figure 5.12. Figure 5.12: Numerical simulation of the Re(u (1))-profile of the one-soliton solution of the 2-CSP equation. maxerr(Re(u (1))) = 2.87 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p040_5_12.png] view at source ↗
Figure 5.13
Figure 5.13. Figure 5.13: Numerical simulation of the Im(u (1))-profile of the one-soliton solution of the 2-CSP equation. maxerr(Im(u (1))) = 4.22 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p041_5_13.png] view at source ↗
Figure 5.14
Figure 5.14. Figure 5.14: Numerical simulation of the |u (2)|-profile of the one-soliton solution of the 2-CSP equation. maxerr(|u (2)|) = 1.67 × 10−4 [PITH_FULL_IMAGE:figures/full_fig_p041_5_14.png] view at source ↗
Figure 5.15
Figure 5.15. Figure 5.15: Numerical simulation of the Re(u (2))-profile of the one-soliton solution of the 2-CSP equation. maxerr(Re(u (2))) = 2.85 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p042_5_15.png] view at source ↗
Figure 5.16
Figure 5.16. Figure 5.16: Numerical simulation of the Im(u (2))-profile of the one-soliton solution of the 2-CSP equation. maxerr(Im(u (2))) = 4.21 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p042_5_16.png] view at source ↗
Figure 5.17
Figure 5.17. Figure 5.17: Numerical simulation of the |u (1)|-profile of the two-soliton solution of the 2-CSP equation. maxerr(|u (1)|) = 1.63 × 10−4 . These numerical experiments show that the mesh automatically refines in regions where the solutions vary rapidly, confirming that the fully discrete MCSP system (4.30) provides a self-adaptive moving mesh scheme. For the parameter sets considered, the numerical and exact solutio… view at source ↗
Figure 5.18
Figure 5.18. Figure 5.18: Numerical simulation of the Re(u (1))-profile of the two-soliton solution of the 2-CSP equation. maxerr(Re(u (1))) = 2.27 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p044_5_18.png] view at source ↗
Figure 5.19
Figure 5.19. Figure 5.19: Numerical simulation of the Im(u (1))-profile of the two-soliton solution of the 2-CSP equation. maxerr(Im(u (1))) = 2.44 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p044_5_19.png] view at source ↗
Figure 5.20
Figure 5.20. Figure 5.20: Numerical simulation of the |u (2)|-profile of the two-soliton solution of the 2-CSP equation. maxerr(|u (2)|) = 1.02 × 10−4 [PITH_FULL_IMAGE:figures/full_fig_p045_5_20.png] view at source ↗
Figure 5.21
Figure 5.21. Figure 5.21: Numerical simulation of the Re(u (2))-profile of the two-soliton solution of the 2-CSP equation. maxerr(Re(u (2))) = 2.27 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p045_5_21.png] view at source ↗
Figure 5.22
Figure 5.22. Figure 5.22: Numerical simulation of the Im(u (2))-profile of the two-soliton solution of the 2-CSP equation. maxerr(Im(u (2))) = 2.27 × 10−3 . error relative to the one-soliton cases. 6. Conclusion We have established a new formulation of the MCSP equation that naturally encapsulates the CCSP equation as a reduction. Using Hirota’s bilinear method, we derived a bilinear form of the MCSP equation and constructed its… view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

41 extracted references · 1 linked inside Pith

  1. [1]

    Hirota R 1977 Nonlinear partial difference equations. I. A difference analogue of the Korteweg-de Vries equation J. Phys. Soc. Jpn. 43 1424–1433

    1977

  2. [2]

    Hirota R 1977 Nonlinear partial difference equations. I I . Discrete–time Toda equation J. Phys. Soc. Jpn. 43 2074–2078

    1977

  3. [3]

    Hirota R 1977 Nonlinear partial difference equations. I I I. Discrete sine–Gordon equation J. Phys. Soc. Jpn. 43 2079–2086

    1977

  4. [4]

    Hirota R 1978 Nonlinear partial difference equations. I V . B ¨acklund transformation for the discrete-time Toda equation J. Phys. Soc. Jpn. 45 321–332

    1978

  5. [5]

    Hirota R 1979 Nonlinear partial difference equations. V . Nonlinear equations reducible to linear equations J. Phys. Soc. Jpn. 46 312–319 Integrable full discretization of the multi-component sho rt pulse equation 48

    1979

  6. [6]

    Hirota R 1981 Discrete analogue of a generalized Toda equ ation J. Phys. Soc. Jpn. 50 3785–3791

    1981

  7. [7]

    Ablowitz M J and Ladik J F 1975 Nonlinear differential-diff erence equations J. Math. Phys. 16 598–603

    1975

  8. [8]

    Ablowitz M J and Ladik J F 1976 Nonlinear differential-diff erence equations and Fourier analysis J. Math. Phys. 17 1011–1018

    1976

  9. [9]

    Ablowitz M J and Ladik J F 1976 A nonlinear difference schem e and inverse scattering Stud. Appl. Math. 55 213–229

    1976

  10. [10]

    Ablowitz M J and Ladik J F 1977 On the solution of a class of nonlinear partial difference equations Stud. Appl. Math. 57 1–12

    1977

  11. [11]

    Ablowitz M J and Segur H 1981 Solitons and the Inverse Scattering Transform (Philadelphia: Society for Industrial and Applied Mathematics)

    1981

  12. [12]

    Suris Y B 2003 The Problem of Integrable Discretization: Hamiltonian App roach (Basel: Springer)

    2003

  13. [13]

    Hietarinta J, Joshi N and Nijhoff F W 2016 Discrete Systems and Integrability (Cambridge: Cambridge University Press)

    2016

  14. [14]

    Wadati M, Konno K and Ichikawa Y 1979 New integrable nonl inear evolution equations J. Phys. Soc. Jpn. 47 1698–1700

    1979

  15. [15]

    Konno K, Ichikawa Y and Wadati M 1981 A loop soliton propa gating along a stretched rope J. Phys. Soc. Jpn. 50 1025–1026

    1981

  16. [16]

    Ishimori Y 1981 On the modified Korteweg-de Vries solito n and the loop soliton J. Phys. Soc. Jpn. 50 2471–2472

    1981

  17. [17]

    Ishimori Y 1982 A relationship between the Ablowitz-Ka up-Newell-Segur and Wadati-Konno-Ichikawa schemes of the inverse scattering method J. Phys. Soc. Jpn. 51 3036–3041

    1982

  18. [18]

    Wadati M and Sogo K 1983 Gauge transformations in solito n theory J. Phys. Soc. Jpn. 52 394–398

    1983

  19. [19]

    Rogers C and Schief W K 2002 B¨acklund and Darboux Transformations: Geometry and Modern Applications in Soliton Theory (Cambridge: Cambridge University Press)

    2002

  20. [20]

    Feng B F, Maruno K and Ohta Y 2010 Integrable discretizat ions of the short pulse equation J. Phys. A: Math. Theor. 43 085203

    2010

  21. [21]

    Feng B F, Inoguchi J, Kajiwara K, Maruno K and Ohta Y 2013 I ntegrable discretizations of the Dym equation Front. Math. China 8 1017–1029

    2013

  22. [22]

    Feng B F, Maruno K and Ohta Y 2014 Self-adaptive moving me sh schemes for short pulse type equations and their Lax pairs Pac. J. Math. Ind. 6 1–14

    2014

  23. [23]

    Ohta Y , Maruno K and Feng B F 2008 An integrable semi-disc retization of the Camassa-Holm equation and its determinant solution J. Phys. A: Math. Theor. 41 355205

    2008

  24. [24]

    Feng B F, Maruno K and Ohta Y 2010 A self-adaptive moving mesh method for the Camassa-Holm equation J. Comput. Appl. Math. 235 229–243

    2010

  25. [25]

    Feng B F, Maruno K and Ohta Y 2010 Integrable discretizat ions for the short-wave model of the Camassa- Holm equation J. Phys. A: Math. Theor. 43 265202

    2010

  26. [26]

    Feng B F, Maruno K and Ohta Y 2015 Integrable semi-discre tizations of the reduced Ostrovsky equation J. Phys. A: Math. Theor. 48 135203

    2015

  27. [27]

    Feng B F, Maruno K and Ohta Y 2017 An integrable semi-discrete Degasperis-Procesi equationNonlinearity 30 2246–2267

    2017

  28. [28]

    Sch ¨afer T and Wayne C E 2004 Propagation of ultra-short optical p ulses in cubic nonlinear media Physica D 196 90–105

    2004

  29. [29]

    Chung Y , Jones C K R T, Sch ¨afer T and Wayne C E 2005 Ultra-short pulses in linear and nonl inear media Nonlinearity 18 1351–1374

    2005

  30. [30]

    Sakovich A and Sakovich S 2005 The short pulse equation i s integrable J. Phys. Soc. Jpn. 74 239–241

    2005

  31. [31]

    Sakovich A and Sakovich S 2006 Solitary wave solutions o f the short pulse equation J. Phys. A: Math. Gen. 39 L361–L367

    2006

  32. [32]

    Matsuno Y 2007 Multiloop soliton and multibreather sol utions of the short pulse model equation J. Phys. Soc. Jpn. 76 084003

    2007

  33. [33]

    Hori A, Maruno K, Ohta Y and Feng B F 2025 Integrable self- adaptive moving mesh schemes for multi- component short pulse type equations with nonzero boundary values arXiv:2507.03880 Integrable full discretization of the multi-component sho rt pulse equation 49

    Pith/arXiv arXiv 2025

  34. [34]

    Feng B F, Inoguchi J, Kajiwara K, Maruno K and Ohta Y 2011 D iscrete integrable systems and hodograph transformations arising from motions of discrete plane cur ves J. Phys. A: Math. Theor. 44 395201

    2011

  35. [35]

    Matsuno Y 2011 A novel multi-component generalization of the short pulse equation and its multisoliton solutions J. Math. Phys. 52 123702

    2011

  36. [36]

    Feng B F 2012 An integrable coupled short pulse equation J. Phys. A: Math. Theor. 45 085202

    2012

  37. [37]

    Feng B F, Maruno K and Ohta Y 2015 Integrable semi-discre tization of a multi-component short pulse equation J. Math. Phys. 56 043502

    2015

  38. [38]

    Feng B F, Chen J, Chen Y , Maruno K and Ohta Y 2015 Integrabl e discretizations and self-adaptive moving mesh method for a coupled short pulse equation J. Phys. A: Math. Theor. 48 385202

    2015

  39. [39]

    Feng B F 2015 Complex short pulse and coupled complex sho rt pulse equations Physica D 297 62–75

    2015

  40. [40]

    Hirota R 2000 Discretization of coupled modified KdV equ ations Chaos, Solitons & Fractals 11 77–84

    2000

  41. [41]

    Hirota R 2004 The Direct Method in Soliton Theory (Cambridge: Cambridge University Press)

    2004