code: def test_te_3smooths():
      from patsy.highlevel import incr_dbuilder, build_design_matrices
  
      x1 = (-1.5) ** np.arange(20)
      x2 = (1.6) ** np.arange(20)
      x3 = (-1.2) ** np.arange(20)
      # Hard coded R results for smooth:  te(x1, x2, x3, bs=c("cr", "cs", "cc"), k=c(3,3,4))
      design_matrix_R = np.array(
          [
              [
                  7.2077663709837084334e-05,
                  2.0648333344343273131e-03,
                  -4.7934014082310591768e-04,
                  2.3923430783992746568e-04,
                  6.8534265421922660466e-03,
                  -1.5909867344112936776e-03,
                  -6.8057712777151204314e-09,
                  -1.9496724335203412851e-07,
                  4.5260614658693259131e-08,
                  0.0101479754187435277507,
                  0.290712501531622591333,
                  -0.067487370093906928759,
                  0.03368233306025386619709,
                  0.9649092451763204847381,
                  -0.2239985793289433757547,
                  -9.5819975394704535133e-07,
                  -2.7449874082511405643e-05,
                  6.3723431275833230217e-06,
                  -1.5205851762850489204e-04,
                  -0.00435607204539782688624,
                  0.00101123909269346416370,
                  -5.0470024059694933508e-04,
                  -1.4458319360584082416e-02,
                  3.3564223914790921634e-03,
                  1.4357783514933466209e-08,
                  4.1131230514870551983e-07,
                  -9.5483976834512651038e-08,
              ]
          ]
      )
      new_data = {
          "x1": -38.443359375000000000,
          "x2": 68.719476736000032702,
          "x3": -5.1597803519999985156,
      }
      data_chunked = [
          {"x1": x1[:10], "x2": x2[:10], "x3": x3[:10]},
          {"x1": x1[10:], "x2": x2[10:], "x3": x3[10:]},
      ]
      builder = incr_dbuilder(
          "te(cr(x1, df=3), cr(x2, df=3), cc(x3, df=3)) - 1",
          lambda: iter(data_chunked),
      )
      design_matrix = build_design_matrices([builder], new_data)[0]
      assert np.allclose(
          design_matrix, design_matrix_R, rtol=1e-12, atol=0.0
      )
  
imports: import numpy
  from patsy.util import have_pandas, atleast_2d_column_default, no_pickling, assert_no_pickling, safe_string_eq
  from patsy.state import stateful_transform
label_desc: Use image node_id 24 for calling a global function with example usage: test_te_3smooths() without return types