YES
We show the termination of the TRS R:

  le(|0|(),y) -> true()
  le(s(x),|0|()) -> false()
  le(s(x),s(y)) -> le(x,y)
  minus(|0|(),y) -> |0|()
  minus(s(x),y) -> if_minus(le(s(x),y),s(x),y)
  if_minus(true(),s(x),y) -> |0|()
  if_minus(false(),s(x),y) -> s(minus(x,y))
  mod(|0|(),y) -> |0|()
  mod(s(x),|0|()) -> |0|()
  mod(s(x),s(y)) -> if_mod(le(y,x),s(x),s(y))
  if_mod(true(),s(x),s(y)) -> mod(minus(x,y),s(y))
  if_mod(false(),s(x),s(y)) -> s(x)


-- SCC decomposition.

Consider the dependency pair problem (P, R), where P consists of

p1: le#(s(x),s(y)) -> le#(x,y)
p2: minus#(s(x),y) -> if_minus#(le(s(x),y),s(x),y)
p3: minus#(s(x),y) -> le#(s(x),y)
p4: if_minus#(false(),s(x),y) -> minus#(x,y)
p5: mod#(s(x),s(y)) -> if_mod#(le(y,x),s(x),s(y))
p6: mod#(s(x),s(y)) -> le#(y,x)
p7: if_mod#(true(),s(x),s(y)) -> mod#(minus(x,y),s(y))
p8: if_mod#(true(),s(x),s(y)) -> minus#(x,y)

and R consists of:

r1: le(|0|(),y) -> true()
r2: le(s(x),|0|()) -> false()
r3: le(s(x),s(y)) -> le(x,y)
r4: minus(|0|(),y) -> |0|()
r5: minus(s(x),y) -> if_minus(le(s(x),y),s(x),y)
r6: if_minus(true(),s(x),y) -> |0|()
r7: if_minus(false(),s(x),y) -> s(minus(x,y))
r8: mod(|0|(),y) -> |0|()
r9: mod(s(x),|0|()) -> |0|()
r10: mod(s(x),s(y)) -> if_mod(le(y,x),s(x),s(y))
r11: if_mod(true(),s(x),s(y)) -> mod(minus(x,y),s(y))
r12: if_mod(false(),s(x),s(y)) -> s(x)

The estimated dependency graph contains the following SCCs:

  {p5, p7}
  {p2, p4}
  {p1}


-- Reduction pair.

Consider the dependency pair problem (P, R), where P consists of

p1: if_mod#(true(),s(x),s(y)) -> mod#(minus(x,y),s(y))
p2: mod#(s(x),s(y)) -> if_mod#(le(y,x),s(x),s(y))

and R consists of:

r1: le(|0|(),y) -> true()
r2: le(s(x),|0|()) -> false()
r3: le(s(x),s(y)) -> le(x,y)
r4: minus(|0|(),y) -> |0|()
r5: minus(s(x),y) -> if_minus(le(s(x),y),s(x),y)
r6: if_minus(true(),s(x),y) -> |0|()
r7: if_minus(false(),s(x),y) -> s(minus(x,y))
r8: mod(|0|(),y) -> |0|()
r9: mod(s(x),|0|()) -> |0|()
r10: mod(s(x),s(y)) -> if_mod(le(y,x),s(x),s(y))
r11: if_mod(true(),s(x),s(y)) -> mod(minus(x,y),s(y))
r12: if_mod(false(),s(x),s(y)) -> s(x)

The set of usable rules consists of

  r1, r4, r5, r6, r7, r8, r9

Take the reduction pair:

  generalized weighted path order (>=MG, >MG)
  
    base order 1:
  
      max/plus interpretation on natural numbers:
      
        if_mod#_A(x1,x2,x3) = max{x2 + 3, x3 + 17}
        true_A = 2
        s_A(x1) = x1 + 21
        mod#_A(x1,x2) = max{22, x1 + 4, x2 + 17}
        minus_A(x1,x2) = x1 + 13
        le_A(x1,x2) = max{4, x2}
        |0|_A = 1
        if_minus_A(x1,x2,x3) = max{20, x2 + 13}
        false_A = 5
        mod_A(x1,x2) = max{4, x2 - 38}
  
    base order 2:
  
      precedence:
      
        if_mod# = mod# = false > minus > if_minus = mod > s > le > true > |0|
      
  
    partial status:
  
      pi(if_mod#) = [2]
      pi(true) = []
      pi(s) = [1]
      pi(mod#) = [1]
      pi(minus) = []
      pi(le) = [2]
      pi(|0|) = []
      pi(if_minus) = []
      pi(false) = []
      pi(mod) = []

The following rules are removed:

  p1



-- SCC decomposition.

Consider the dependency pair problem (P, R), where P consists of

p1: mod#(s(x),s(y)) -> if_mod#(le(y,x),s(x),s(y))

and R consists of:

r1: le(|0|(),y) -> true()
r2: le(s(x),|0|()) -> false()
r3: le(s(x),s(y)) -> le(x,y)
r4: minus(|0|(),y) -> |0|()
r5: minus(s(x),y) -> if_minus(le(s(x),y),s(x),y)
r6: if_minus(true(),s(x),y) -> |0|()
r7: if_minus(false(),s(x),y) -> s(minus(x,y))
r8: mod(|0|(),y) -> |0|()
r9: mod(s(x),|0|()) -> |0|()
r10: mod(s(x),s(y)) -> if_mod(le(y,x),s(x),s(y))
r11: if_mod(true(),s(x),s(y)) -> mod(minus(x,y),s(y))
r12: if_mod(false(),s(x),s(y)) -> s(x)

The estimated dependency graph contains the following SCCs:

  (no SCCs)

-- Reduction pair.

Consider the dependency pair problem (P, R), where P consists of

p1: if_minus#(false(),s(x),y) -> minus#(x,y)
p2: minus#(s(x),y) -> if_minus#(le(s(x),y),s(x),y)

and R consists of:

r1: le(|0|(),y) -> true()
r2: le(s(x),|0|()) -> false()
r3: le(s(x),s(y)) -> le(x,y)
r4: minus(|0|(),y) -> |0|()
r5: minus(s(x),y) -> if_minus(le(s(x),y),s(x),y)
r6: if_minus(true(),s(x),y) -> |0|()
r7: if_minus(false(),s(x),y) -> s(minus(x,y))
r8: mod(|0|(),y) -> |0|()
r9: mod(s(x),|0|()) -> |0|()
r10: mod(s(x),s(y)) -> if_mod(le(y,x),s(x),s(y))
r11: if_mod(true(),s(x),s(y)) -> mod(minus(x,y),s(y))
r12: if_mod(false(),s(x),s(y)) -> s(x)

The set of usable rules consists of

  r1, r2, r3, r4, r6, r8, r9

Take the reduction pair:

  generalized weighted path order (>=MG, >MG)
  
    base order 1:
  
      max/plus interpretation on natural numbers:
      
        if_minus#_A(x1,x2,x3) = max{38, x2 + 1}
        false_A = 15
        s_A(x1) = max{37, x1 + 20}
        minus#_A(x1,x2) = x1 + 19
        le_A(x1,x2) = max{14, x1 - 21, x2 - 38}
        |0|_A = 2
        true_A = 1
        minus_A(x1,x2) = max{36, x1 - 19}
        if_minus_A(x1,x2,x3) = max{36, x1 - 17, x2 - 38, x3 - 1}
        mod_A(x1,x2) = max{3, x1 - 59, x2 - 49}
  
    base order 2:
  
      precedence:
      
        if_minus# > s = minus# = le = minus = if_minus = mod > true > |0| > false
      
  
    partial status:
  
      pi(if_minus#) = [2]
      pi(false) = []
      pi(s) = [1]
      pi(minus#) = [1]
      pi(le) = []
      pi(|0|) = []
      pi(true) = []
      pi(minus) = []
      pi(if_minus) = []
      pi(mod) = []

The following rules are removed:

  p2



-- SCC decomposition.

Consider the dependency pair problem (P, R), where P consists of

p1: if_minus#(false(),s(x),y) -> minus#(x,y)

and R consists of:

r1: le(|0|(),y) -> true()
r2: le(s(x),|0|()) -> false()
r3: le(s(x),s(y)) -> le(x,y)
r4: minus(|0|(),y) -> |0|()
r5: minus(s(x),y) -> if_minus(le(s(x),y),s(x),y)
r6: if_minus(true(),s(x),y) -> |0|()
r7: if_minus(false(),s(x),y) -> s(minus(x,y))
r8: mod(|0|(),y) -> |0|()
r9: mod(s(x),|0|()) -> |0|()
r10: mod(s(x),s(y)) -> if_mod(le(y,x),s(x),s(y))
r11: if_mod(true(),s(x),s(y)) -> mod(minus(x,y),s(y))
r12: if_mod(false(),s(x),s(y)) -> s(x)

The estimated dependency graph contains the following SCCs:

  (no SCCs)

-- Reduction pair.

Consider the dependency pair problem (P, R), where P consists of

p1: le#(s(x),s(y)) -> le#(x,y)

and R consists of:

r1: le(|0|(),y) -> true()
r2: le(s(x),|0|()) -> false()
r3: le(s(x),s(y)) -> le(x,y)
r4: minus(|0|(),y) -> |0|()
r5: minus(s(x),y) -> if_minus(le(s(x),y),s(x),y)
r6: if_minus(true(),s(x),y) -> |0|()
r7: if_minus(false(),s(x),y) -> s(minus(x,y))
r8: mod(|0|(),y) -> |0|()
r9: mod(s(x),|0|()) -> |0|()
r10: mod(s(x),s(y)) -> if_mod(le(y,x),s(x),s(y))
r11: if_mod(true(),s(x),s(y)) -> mod(minus(x,y),s(y))
r12: if_mod(false(),s(x),s(y)) -> s(x)

The set of usable rules consists of

  r1, r2, r4, r6, r8, r9

Take the reduction pair:

  generalized weighted path order (>=MG, >MG)
  
    base order 1:
  
      max/plus interpretation on natural numbers:
      
        le#_A(x1,x2) = max{1, x1 - 11, x2 - 22}
        s_A(x1) = max{21, x1 + 10}
        le_A(x1,x2) = max{8, x1 - 25, x2 - 23}
        |0|_A = 2
        true_A = 3
        false_A = 1
        minus_A(x1,x2) = max{20, x1 - 22, x2 - 22}
        if_minus_A(x1,x2,x3) = max{20, x1 - 10, x2 - 32, x3 - 22}
        mod_A(x1,x2) = max{9, x1 - 31, x2 - 34}
  
    base order 2:
  
      precedence:
      
        if_minus > minus = mod > le# = le > s > true > |0| > false
      
  
    partial status:
  
      pi(le#) = []
      pi(s) = [1]
      pi(le) = []
      pi(|0|) = []
      pi(true) = []
      pi(false) = []
      pi(minus) = []
      pi(if_minus) = []
      pi(mod) = []

The following rules are removed:

  p1

Then no dependency pair remains.