/ Hex Artifact Content
Login

Artifact 148dd07286b16e50f232bb638a47850085ad37d51f270429905bd865e595d1ca:


0000: 2f 2a 0a 2a 2a 20 32 30 31 33 2d 30 35 2d 32 38  /*.** 2013-05-28
0010: 0a 2a 2a 0a 2a 2a 20 54 68 65 20 61 75 74 68 6f  .**.** The autho
0020: 72 20 64 69 73 63 6c 61 69 6d 73 20 63 6f 70 79  r disclaims copy
0030: 72 69 67 68 74 20 74 6f 20 74 68 69 73 20 73 6f  right to this so
0040: 75 72 63 65 20 63 6f 64 65 2e 20 20 49 6e 20 70  urce code.  In p
0050: 6c 61 63 65 20 6f 66 0a 2a 2a 20 61 20 6c 65 67  lace of.** a leg
0060: 61 6c 20 6e 6f 74 69 63 65 2c 20 68 65 72 65 20  al notice, here 
0070: 69 73 20 61 20 62 6c 65 73 73 69 6e 67 3a 0a 2a  is a blessing:.*
0080: 2a 0a 2a 2a 20 20 20 20 4d 61 79 20 79 6f 75 20  *.**    May you 
0090: 64 6f 20 67 6f 6f 64 20 61 6e 64 20 6e 6f 74 20  do good and not 
00a0: 65 76 69 6c 2e 0a 2a 2a 20 20 20 20 4d 61 79 20  evil..**    May 
00b0: 79 6f 75 20 66 69 6e 64 20 66 6f 72 67 69 76 65  you find forgive
00c0: 6e 65 73 73 20 66 6f 72 20 79 6f 75 72 73 65 6c  ness for yoursel
00d0: 66 20 61 6e 64 20 66 6f 72 67 69 76 65 20 6f 74  f and forgive ot
00e0: 68 65 72 73 2e 0a 2a 2a 20 20 20 20 4d 61 79 20  hers..**    May 
00f0: 79 6f 75 20 73 68 61 72 65 20 66 72 65 65 6c 79  you share freely
0100: 2c 20 6e 65 76 65 72 20 74 61 6b 69 6e 67 20 6d  , never taking m
0110: 6f 72 65 20 74 68 61 6e 20 79 6f 75 20 67 69 76  ore than you giv
0120: 65 2e 0a 2a 2a 0a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a  e..**.**********
0130: 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a  ****************
0140: 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a  ****************
0150: 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a  ****************
0160: 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a  ****************
0170: 2a 2a 2a 2a 0a 2a 2a 0a 2a 2a 20 54 68 69 73 20  ****.**.** This 
0180: 66 69 6c 65 20 63 6f 6e 74 61 69 6e 73 20 63 6f  file contains co
0190: 64 65 20 74 6f 20 69 6d 70 6c 65 6d 65 6e 74 20  de to implement 
01a0: 74 68 65 20 70 65 72 63 65 6e 74 69 6c 65 28 59  the percentile(Y
01b0: 2c 50 29 20 53 51 4c 20 66 75 6e 63 74 69 6f 6e  ,P) SQL function
01c0: 0a 2a 2a 20 61 73 20 64 65 73 63 72 69 62 65 64  .** as described
01d0: 20 62 65 6c 6f 77 3a 0a 2a 2a 0a 2a 2a 20 20 20   below:.**.**   
01e0: 28 31 29 20 20 54 68 65 20 70 65 72 63 65 6e 74  (1)  The percent
01f0: 69 6c 65 28 59 2c 50 29 20 66 75 6e 63 74 69 6f  ile(Y,P) functio
0200: 6e 20 69 73 20 61 6e 20 61 67 67 72 65 67 61 74  n is an aggregat
0210: 65 20 66 75 6e 63 74 69 6f 6e 20 74 61 6b 69 6e  e function takin
0220: 67 0a 2a 2a 20 20 20 20 20 20 20 20 65 78 61 63  g.**        exac
0230: 74 6c 79 20 74 77 6f 20 61 72 67 75 6d 65 6e 74  tly two argument
0240: 73 2e 0a 2a 2a 0a 2a 2a 20 20 20 28 32 29 20 20  s..**.**   (2)  
0250: 49 66 20 74 68 65 20 50 20 61 72 67 75 6d 65 6e  If the P argumen
0260: 74 20 74 6f 20 70 65 72 63 65 6e 74 69 6c 65 28  t to percentile(
0270: 59 2c 50 29 20 69 73 20 6e 6f 74 20 74 68 65 20  Y,P) is not the 
0280: 73 61 6d 65 20 66 6f 72 20 65 76 65 72 79 0a 2a  same for every.*
0290: 2a 20 20 20 20 20 20 20 20 72 6f 77 20 69 6e 20  *        row in 
02a0: 74 68 65 20 61 67 67 72 65 67 61 74 65 20 74 68  the aggregate th
02b0: 65 6e 20 61 6e 20 65 72 72 6f 72 20 69 73 20 74  en an error is t
02c0: 68 72 6f 77 6e 2e 20 20 54 68 65 20 77 6f 72 64  hrown.  The word
02d0: 20 22 73 61 6d 65 22 0a 2a 2a 20 20 20 20 20 20   "same".**      
02e0: 20 20 69 6e 20 74 68 65 20 70 72 65 76 69 6f 75    in the previou
02f0: 73 20 73 65 6e 74 65 6e 63 65 20 6d 65 61 6e 73  s sentence means
0300: 20 74 68 61 74 20 74 68 65 20 76 61 6c 75 65 20   that the value 
0310: 64 69 66 66 65 72 20 62 79 20 6c 65 73 73 0a 2a  differ by less.*
0320: 2a 20 20 20 20 20 20 20 20 74 68 61 6e 20 30 2e  *        than 0.
0330: 30 30 31 2e 0a 2a 2a 0a 2a 2a 20 20 20 28 33 29  001..**.**   (3)
0340: 20 20 49 66 20 74 68 65 20 50 20 61 72 67 75 6d    If the P argum
0350: 65 6e 74 20 74 6f 20 70 65 72 63 65 6e 74 69 6c  ent to percentil
0360: 65 28 59 2c 50 29 20 65 76 61 6c 75 61 74 65 73  e(Y,P) evaluates
0370: 20 74 6f 20 61 6e 79 74 68 69 6e 67 20 6f 74 68   to anything oth
0380: 65 72 0a 2a 2a 20 20 20 20 20 20 20 20 74 68 61  er.**        tha
0390: 6e 20 61 20 6e 75 6d 62 65 72 20 69 6e 20 74 68  n a number in th
03a0: 65 20 72 61 6e 67 65 20 6f 66 20 30 2e 30 20 74  e range of 0.0 t
03b0: 6f 20 31 30 30 2e 30 20 69 6e 63 6c 75 73 69 76  o 100.0 inclusiv
03c0: 65 20 74 68 65 6e 20 61 6e 0a 2a 2a 20 20 20 20  e then an.**    
03d0: 20 20 20 20 65 72 72 6f 72 20 69 73 20 74 68 72      error is thr
03e0: 6f 77 6e 2e 0a 2a 2a 0a 2a 2a 20 20 20 28 34 29  own..**.**   (4)
03f0: 20 20 49 66 20 61 6e 79 20 59 20 61 72 67 75 6d    If any Y argum
0400: 65 6e 74 20 74 6f 20 70 65 72 63 65 6e 74 69 6c  ent to percentil
0410: 65 28 59 2c 50 29 20 65 76 61 6c 75 61 74 65 73  e(Y,P) evaluates
0420: 20 74 6f 20 61 20 76 61 6c 75 65 20 74 68 61 74   to a value that
0430: 0a 2a 2a 20 20 20 20 20 20 20 20 69 73 20 6e 6f  .**        is no
0440: 74 20 4e 55 4c 4c 20 61 6e 64 20 69 73 20 6e 6f  t NULL and is no
0450: 74 20 6e 75 6d 65 72 69 63 20 74 68 65 6e 20 61  t numeric then a
0460: 6e 20 65 72 72 6f 72 20 69 73 20 74 68 72 6f 77  n error is throw
0470: 6e 2e 0a 2a 2a 0a 2a 2a 20 20 20 28 35 29 20 20  n..**.**   (5)  
0480: 49 66 20 61 6e 79 20 59 20 61 72 67 75 6d 65 6e  If any Y argumen
0490: 74 20 74 6f 20 70 65 72 63 65 6e 74 69 6c 65 28  t to percentile(
04a0: 59 2c 50 29 20 65 76 61 6c 75 61 74 65 73 20 74  Y,P) evaluates t
04b0: 6f 20 70 6c 75 73 20 6f 72 20 6d 69 6e 75 73 0a  o plus or minus.
04c0: 2a 2a 20 20 20 20 20 20 20 20 69 6e 66 69 6e 69  **        infini
04d0: 74 79 20 74 68 65 6e 20 61 6e 20 65 72 72 6f 72  ty then an error
04e0: 20 69 73 20 74 68 72 6f 77 6e 2e 20 20 28 53 51   is thrown.  (SQ
04f0: 4c 69 74 65 20 61 6c 77 61 79 73 20 69 6e 74 65  Lite always inte
0500: 72 70 72 65 74 73 20 4e 61 4e 0a 2a 2a 20 20 20  rprets NaN.**   
0510: 20 20 20 20 20 76 61 6c 75 65 73 20 61 73 20 4e       values as N
0520: 55 4c 4c 2e 29 0a 2a 2a 0a 2a 2a 20 20 20 28 36  ULL.).**.**   (6
0530: 29 20 20 42 6f 74 68 20 59 20 61 6e 64 20 50 20  )  Both Y and P 
0540: 69 6e 20 70 65 72 63 65 6e 74 69 6c 65 28 59 2c  in percentile(Y,
0550: 50 29 20 63 61 6e 20 62 65 20 61 72 62 69 74 72  P) can be arbitr
0560: 61 72 79 20 65 78 70 72 65 73 73 69 6f 6e 73 2c  ary expressions,
0570: 0a 2a 2a 20 20 20 20 20 20 20 20 69 6e 63 6c 75  .**        inclu
0580: 64 69 6e 67 20 43 41 53 45 20 57 48 45 4e 20 65  ding CASE WHEN e
0590: 78 70 72 65 73 73 69 6f 6e 73 2e 0a 2a 2a 0a 2a  xpressions..**.*
05a0: 2a 20 20 20 28 37 29 20 20 54 68 65 20 70 65 72  *   (7)  The per
05b0: 63 65 6e 74 69 6c 65 28 59 2c 50 29 20 61 67 67  centile(Y,P) agg
05c0: 72 65 67 61 74 65 20 69 73 20 61 62 6c 65 20 74  regate is able t
05d0: 6f 20 68 61 6e 64 6c 65 20 69 6e 70 75 74 73 20  o handle inputs 
05e0: 6f 66 20 61 74 20 6c 65 61 73 74 0a 2a 2a 20 20  of at least.**  
05f0: 20 20 20 20 20 20 6f 6e 65 20 6d 69 6c 6c 69 6f        one millio
0600: 6e 20 28 31 2c 30 30 30 2c 30 30 30 29 20 72 6f  n (1,000,000) ro
0610: 77 73 2e 0a 2a 2a 0a 2a 2a 20 20 20 28 38 29 20  ws..**.**   (8) 
0620: 20 49 66 20 74 68 65 72 65 20 61 72 65 20 6e 6f   If there are no
0630: 20 6e 6f 6e 2d 4e 55 4c 4c 20 76 61 6c 75 65 73   non-NULL values
0640: 20 66 6f 72 20 59 2c 20 74 68 65 6e 20 70 65 72   for Y, then per
0650: 63 65 6e 74 69 6c 65 28 59 2c 50 29 0a 2a 2a 20  centile(Y,P).** 
0660: 20 20 20 20 20 20 20 72 65 74 75 72 6e 73 20 4e         returns N
0670: 55 4c 4c 2e 0a 2a 2a 0a 2a 2a 20 20 20 28 39 29  ULL..**.**   (9)
0680: 20 20 49 66 20 74 68 65 72 65 20 69 73 20 65 78    If there is ex
0690: 61 63 74 6c 79 20 6f 6e 65 20 6e 6f 6e 2d 4e 55  actly one non-NU
06a0: 4c 4c 20 76 61 6c 75 65 20 66 6f 72 20 59 2c 20  LL value for Y, 
06b0: 74 68 65 20 70 65 72 63 65 6e 74 69 6c 65 28 59  the percentile(Y
06c0: 2c 50 29 0a 2a 2a 20 20 20 20 20 20 20 20 72 65  ,P).**        re
06d0: 74 75 72 6e 73 20 74 68 65 20 6f 6e 65 20 59 20  turns the one Y 
06e0: 76 61 6c 75 65 2e 0a 2a 2a 0a 2a 2a 20 20 28 31  value..**.**  (1
06f0: 30 29 20 20 49 66 20 74 68 65 72 65 20 4e 20 6e  0)  If there N n
0700: 6f 6e 2d 4e 55 4c 4c 20 76 61 6c 75 65 73 20 6f  on-NULL values o
0710: 66 20 59 20 77 68 65 72 65 20 4e 20 69 73 20 74  f Y where N is t
0720: 77 6f 20 6f 72 20 6d 6f 72 65 20 61 6e 64 0a 2a  wo or more and.*
0730: 2a 20 20 20 20 20 20 20 20 74 68 65 20 59 20 76  *        the Y v
0740: 61 6c 75 65 73 20 61 72 65 20 6f 72 64 65 72 65  alues are ordere
0750: 64 20 66 72 6f 6d 20 6c 65 61 73 74 20 74 6f 20  d from least to 
0760: 67 72 65 61 74 65 73 74 20 61 6e 64 20 61 20 67  greatest and a g
0770: 72 61 70 68 20 69 73 0a 2a 2a 20 20 20 20 20 20  raph is.**      
0780: 20 20 64 72 61 77 6e 20 66 72 6f 6d 20 30 20 74    drawn from 0 t
0790: 6f 20 4e 2d 31 20 73 75 63 68 20 74 68 61 74 20  o N-1 such that 
07a0: 74 68 65 20 68 65 69 67 68 74 20 6f 66 20 74 68  the height of th
07b0: 65 20 67 72 61 70 68 20 61 74 20 4a 20 69 73 0a  e graph at J is.
07c0: 2a 2a 20 20 20 20 20 20 20 20 74 68 65 20 4a 2d  **        the J-
07d0: 74 68 20 59 20 76 61 6c 75 65 20 61 6e 64 20 73  th Y value and s
07e0: 75 63 68 20 74 68 61 74 20 73 74 72 61 69 67 68  uch that straigh
07f0: 74 20 6c 69 6e 65 73 20 61 72 65 20 64 72 61 77  t lines are draw
0800: 6e 20 62 65 74 77 65 65 6e 0a 2a 2a 20 20 20 20  n between.**    
0810: 20 20 20 20 61 64 6a 61 63 65 6e 74 20 59 20 76      adjacent Y v
0820: 61 6c 75 65 73 2c 20 74 68 65 6e 20 74 68 65 20  alues, then the 
0830: 70 65 72 63 65 6e 74 69 6c 65 28 59 2c 50 29 20  percentile(Y,P) 
0840: 66 75 6e 63 74 69 6f 6e 20 72 65 74 75 72 6e 73  function returns
0850: 0a 2a 2a 20 20 20 20 20 20 20 20 74 68 65 20 68  .**        the h
0860: 65 69 67 68 74 20 6f 66 20 74 68 65 20 67 72 61  eight of the gra
0870: 70 68 20 61 74 20 50 2a 28 4e 2d 31 29 2f 31 30  ph at P*(N-1)/10
0880: 30 2e 0a 2a 2a 0a 2a 2a 20 20 28 31 31 29 20 20  0..**.**  (11)  
0890: 54 68 65 20 70 65 72 63 65 6e 74 69 6c 65 28 59  The percentile(Y
08a0: 2c 50 29 20 66 75 6e 63 74 69 6f 6e 20 61 6c 77  ,P) function alw
08b0: 61 79 73 20 72 65 74 75 72 6e 73 20 65 69 74 68  ays returns eith
08c0: 65 72 20 61 20 66 6c 6f 61 74 69 6e 67 0a 2a 2a  er a floating.**
08d0: 20 20 20 20 20 20 20 20 70 6f 69 6e 74 20 6e 75          point nu
08e0: 6d 62 65 72 20 6f 72 20 4e 55 4c 4c 2e 0a 2a 2a  mber or NULL..**
08f0: 0a 2a 2a 20 20 28 31 32 29 20 20 54 68 65 20 70  .**  (12)  The p
0900: 65 72 63 65 6e 74 69 6c 65 28 59 2c 50 29 20 69  ercentile(Y,P) i
0910: 73 20 69 6d 70 6c 65 6d 65 6e 74 65 64 20 61 73  s implemented as
0920: 20 61 20 73 69 6e 67 6c 65 20 43 39 39 20 73 6f   a single C99 so
0930: 75 72 63 65 2d 63 6f 64 65 0a 2a 2a 20 20 20 20  urce-code.**    
0940: 20 20 20 20 66 69 6c 65 20 74 68 61 74 20 63 6f      file that co
0950: 6d 70 69 6c 65 73 20 69 6e 74 6f 20 61 20 73 68  mpiles into a sh
0960: 61 72 65 64 2d 6c 69 62 72 61 72 79 20 6f 72 20  ared-library or 
0970: 44 4c 4c 20 74 68 61 74 20 63 61 6e 20 62 65 20  DLL that can be 
0980: 6c 6f 61 64 65 64 0a 2a 2a 20 20 20 20 20 20 20  loaded.**       
0990: 20 69 6e 74 6f 20 53 51 4c 69 74 65 20 75 73 69   into SQLite usi
09a0: 6e 67 20 74 68 65 20 73 71 6c 69 74 65 33 5f 6c  ng the sqlite3_l
09b0: 6f 61 64 5f 65 78 74 65 6e 73 69 6f 6e 28 29 20  oad_extension() 
09c0: 69 6e 74 65 72 66 61 63 65 2e 0a 2a 2f 0a 23 69  interface..*/.#i
09d0: 6e 63 6c 75 64 65 20 22 73 71 6c 69 74 65 33 65  nclude "sqlite3e
09e0: 78 74 2e 68 22 0a 53 51 4c 49 54 45 5f 45 58 54  xt.h".SQLITE_EXT
09f0: 45 4e 53 49 4f 4e 5f 49 4e 49 54 31 0a 23 69 6e  ENSION_INIT1.#in
0a00: 63 6c 75 64 65 20 3c 61 73 73 65 72 74 2e 68 3e  clude <assert.h>
0a10: 0a 23 69 6e 63 6c 75 64 65 20 3c 73 74 72 69 6e  .#include <strin
0a20: 67 2e 68 3e 0a 23 69 6e 63 6c 75 64 65 20 3c 73  g.h>.#include <s
0a30: 74 64 6c 69 62 2e 68 3e 0a 0a 2f 2a 20 54 68 65  tdlib.h>../* The
0a40: 20 66 6f 6c 6c 6f 77 69 6e 67 20 6f 62 6a 65 63   following objec
0a50: 74 20 69 73 20 74 68 65 20 73 65 73 73 69 6f 6e  t is the session
0a60: 20 63 6f 6e 74 65 78 74 20 66 6f 72 20 61 20 73   context for a s
0a70: 69 6e 67 6c 65 20 70 65 72 63 65 6e 74 69 6c 65  ingle percentile
0a80: 28 29 0a 2a 2a 20 66 75 6e 63 74 69 6f 6e 2e 20  ().** function. 
0a90: 20 57 65 20 68 61 76 65 20 74 6f 20 72 65 6d 65   We have to reme
0aa0: 6d 62 65 72 20 61 6c 6c 20 69 6e 70 75 74 20 59  mber all input Y
0ab0: 20 76 61 6c 75 65 73 20 75 6e 74 69 6c 20 74 68   values until th
0ac0: 65 20 76 65 72 79 20 65 6e 64 2e 0a 2a 2a 20 54  e very end..** T
0ad0: 68 6f 73 65 20 76 61 6c 75 65 73 20 61 72 65 20  hose values are 
0ae0: 61 63 63 75 6d 75 6c 61 74 65 64 20 69 6e 20 74  accumulated in t
0af0: 68 65 20 50 65 72 63 65 6e 74 69 6c 65 2e 61 5b  he Percentile.a[
0b00: 5d 20 61 72 72 61 79 2e 0a 2a 2f 0a 74 79 70 65  ] array..*/.type
0b10: 64 65 66 20 73 74 72 75 63 74 20 50 65 72 63 65  def struct Perce
0b20: 6e 74 69 6c 65 20 50 65 72 63 65 6e 74 69 6c 65  ntile Percentile
0b30: 3b 0a 73 74 72 75 63 74 20 50 65 72 63 65 6e 74  ;.struct Percent
0b40: 69 6c 65 20 7b 0a 20 20 75 6e 73 69 67 6e 65 64  ile {.  unsigned
0b50: 20 6e 41 6c 6c 6f 63 3b 20 20 20 20 20 2f 2a 20   nAlloc;     /* 
0b60: 4e 75 6d 62 65 72 20 6f 66 20 73 6c 6f 74 73 20  Number of slots 
0b70: 61 6c 6c 6f 63 61 74 65 64 20 66 6f 72 20 61 5b  allocated for a[
0b80: 5d 20 2a 2f 0a 20 20 75 6e 73 69 67 6e 65 64 20  ] */.  unsigned 
0b90: 6e 55 73 65 64 3b 20 20 20 20 20 20 2f 2a 20 4e  nUsed;      /* N
0ba0: 75 6d 62 65 72 20 6f 66 20 73 6c 6f 74 73 20 61  umber of slots a
0bb0: 63 74 75 61 6c 6c 79 20 75 73 65 64 20 69 6e 20  ctually used in 
0bc0: 61 5b 5d 20 2a 2f 0a 20 20 64 6f 75 62 6c 65 20  a[] */.  double 
0bd0: 72 50 63 74 3b 20 20 20 20 20 20 20 20 20 2f 2a  rPct;         /*
0be0: 20 31 2e 30 20 6d 6f 72 65 20 74 68 61 6e 20 74   1.0 more than t
0bf0: 68 65 20 76 61 6c 75 65 20 66 6f 72 20 50 20 2a  he value for P *
0c00: 2f 0a 20 20 64 6f 75 62 6c 65 20 2a 61 3b 20 20  /.  double *a;  
0c10: 20 20 20 20 20 20 20 20 20 2f 2a 20 41 72 72 61           /* Arra
0c20: 79 20 6f 66 20 59 20 76 61 6c 75 65 73 20 2a 2f  y of Y values */
0c30: 0a 7d 3b 0a 0a 2f 2a 0a 2a 2a 20 52 65 74 75 72  .};../*.** Retur
0c40: 6e 20 54 52 55 45 20 69 66 20 74 68 65 20 69 6e  n TRUE if the in
0c50: 70 75 74 20 66 6c 6f 61 74 69 6e 67 2d 70 6f 69  put floating-poi
0c60: 6e 74 20 6e 75 6d 62 65 72 20 69 73 20 61 6e 20  nt number is an 
0c70: 69 6e 66 69 6e 69 74 79 2e 0a 2a 2f 0a 73 74 61  infinity..*/.sta
0c80: 74 69 63 20 69 6e 74 20 69 73 49 6e 66 69 6e 69  tic int isInfini
0c90: 74 79 28 64 6f 75 62 6c 65 20 72 29 7b 0a 20 20  ty(double r){.  
0ca0: 73 71 6c 69 74 65 33 5f 75 69 6e 74 36 34 20 75  sqlite3_uint64 u
0cb0: 3b 0a 20 20 61 73 73 65 72 74 28 20 73 69 7a 65  ;.  assert( size
0cc0: 6f 66 28 75 29 3d 3d 73 69 7a 65 6f 66 28 72 29  of(u)==sizeof(r)
0cd0: 20 29 3b 0a 20 20 6d 65 6d 63 70 79 28 26 75 2c   );.  memcpy(&u,
0ce0: 20 26 72 2c 20 73 69 7a 65 6f 66 28 75 29 29 3b   &r, sizeof(u));
0cf0: 0a 20 20 72 65 74 75 72 6e 20 28 28 75 3e 3e 35  .  return ((u>>5
0d00: 32 29 26 30 78 37 66 66 29 3d 3d 30 78 37 66 66  2)&0x7ff)==0x7ff
0d10: 3b 0a 7d 0a 0a 2f 2a 0a 2a 2a 20 52 65 74 75 72  ;.}../*.** Retur
0d20: 6e 20 54 52 55 45 20 69 66 20 74 77 6f 20 64 6f  n TRUE if two do
0d30: 75 62 6c 65 73 20 64 69 66 66 65 72 20 62 79 20  ubles differ by 
0d40: 30 2e 30 30 31 20 6f 72 20 6c 65 73 73 0a 2a 2f  0.001 or less.*/
0d50: 0a 73 74 61 74 69 63 20 69 6e 74 20 73 61 6d 65  .static int same
0d60: 56 61 6c 75 65 28 64 6f 75 62 6c 65 20 61 2c 20  Value(double a, 
0d70: 64 6f 75 62 6c 65 20 62 29 7b 0a 20 20 61 20 2d  double b){.  a -
0d80: 3d 20 62 3b 0a 20 20 72 65 74 75 72 6e 20 61 3e  = b;.  return a>
0d90: 3d 2d 30 2e 30 30 31 20 26 26 20 61 3c 3d 30 2e  =-0.001 && a<=0.
0da0: 30 30 31 3b 0a 7d 0a 0a 2f 2a 0a 2a 2a 20 54 68  001;.}../*.** Th
0db0: 65 20 22 73 74 65 70 22 20 66 75 6e 63 74 69 6f  e "step" functio
0dc0: 6e 20 66 6f 72 20 70 65 72 63 65 6e 74 69 6c 65  n for percentile
0dd0: 28 59 2c 50 29 20 69 73 20 63 61 6c 6c 65 64 20  (Y,P) is called 
0de0: 6f 6e 63 65 20 66 6f 72 20 65 61 63 68 0a 2a 2a  once for each.**
0df0: 20 69 6e 70 75 74 20 72 6f 77 2e 0a 2a 2f 0a 73   input row..*/.s
0e00: 74 61 74 69 63 20 76 6f 69 64 20 70 65 72 63 65  tatic void perce
0e10: 6e 74 53 74 65 70 28 73 71 6c 69 74 65 33 5f 63  ntStep(sqlite3_c
0e20: 6f 6e 74 65 78 74 20 2a 70 43 74 78 2c 20 69 6e  ontext *pCtx, in
0e30: 74 20 61 72 67 63 2c 20 73 71 6c 69 74 65 33 5f  t argc, sqlite3_
0e40: 76 61 6c 75 65 20 2a 2a 61 72 67 76 29 7b 0a 20  value **argv){. 
0e50: 20 50 65 72 63 65 6e 74 69 6c 65 20 2a 70 3b 0a   Percentile *p;.
0e60: 20 20 64 6f 75 62 6c 65 20 72 50 63 74 3b 0a 20    double rPct;. 
0e70: 20 69 6e 74 20 65 54 79 70 65 3b 0a 20 20 64 6f   int eType;.  do
0e80: 75 62 6c 65 20 79 3b 0a 20 20 61 73 73 65 72 74  uble y;.  assert
0e90: 28 20 61 72 67 63 3d 3d 32 20 29 3b 0a 0a 20 20  ( argc==2 );..  
0ea0: 2f 2a 20 52 65 71 75 69 72 65 6d 65 6e 74 20 33  /* Requirement 3
0eb0: 3a 20 20 50 20 6d 75 73 74 20 62 65 20 61 20 6e  :  P must be a n
0ec0: 75 6d 62 65 72 20 62 65 74 77 65 65 6e 20 30 20  umber between 0 
0ed0: 61 6e 64 20 31 30 30 20 2a 2f 0a 20 20 65 54 79  and 100 */.  eTy
0ee0: 70 65 20 3d 20 73 71 6c 69 74 65 33 5f 76 61 6c  pe = sqlite3_val
0ef0: 75 65 5f 6e 75 6d 65 72 69 63 5f 74 79 70 65 28  ue_numeric_type(
0f00: 61 72 67 76 5b 31 5d 29 3b 0a 20 20 72 50 63 74  argv[1]);.  rPct
0f10: 20 3d 20 73 71 6c 69 74 65 33 5f 76 61 6c 75 65   = sqlite3_value
0f20: 5f 64 6f 75 62 6c 65 28 61 72 67 76 5b 31 5d 29  _double(argv[1])
0f30: 3b 0a 20 20 69 66 28 20 28 65 54 79 70 65 21 3d  ;.  if( (eType!=
0f40: 53 51 4c 49 54 45 5f 49 4e 54 45 47 45 52 20 26  SQLITE_INTEGER &
0f50: 26 20 65 54 79 70 65 21 3d 53 51 4c 49 54 45 5f  & eType!=SQLITE_
0f60: 46 4c 4f 41 54 29 0a 20 20 20 7c 7c 20 72 50 63  FLOAT).   || rPc
0f70: 74 3c 30 2e 30 20 7c 7c 20 72 50 63 74 3e 31 30  t<0.0 || rPct>10
0f80: 30 2e 30 20 29 7b 0a 20 20 20 20 73 71 6c 69 74  0.0 ){.    sqlit
0f90: 65 33 5f 72 65 73 75 6c 74 5f 65 72 72 6f 72 28  e3_result_error(
0fa0: 70 43 74 78 2c 20 22 32 6e 64 20 61 72 67 75 6d  pCtx, "2nd argum
0fb0: 65 6e 74 20 74 6f 20 70 65 72 63 65 6e 74 69 6c  ent to percentil
0fc0: 65 28 29 20 69 73 20 6e 6f 74 20 22 0a 20 20 20  e() is not ".   
0fd0: 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20                  
0fe0: 20 20 20 20 20 20 22 61 20 6e 75 6d 62 65 72 20        "a number 
0ff0: 62 65 74 77 65 65 6e 20 30 2e 30 20 61 6e 64 20  between 0.0 and 
1000: 31 30 30 2e 30 22 2c 20 2d 31 29 3b 0a 20 20 20  100.0", -1);.   
1010: 20 72 65 74 75 72 6e 3b 0a 20 20 7d 0a 0a 20 20   return;.  }..  
1020: 2f 2a 20 41 6c 6c 6f 63 61 74 65 20 74 68 65 20  /* Allocate the 
1030: 73 65 73 73 69 6f 6e 20 63 6f 6e 74 65 78 74 2e  session context.
1040: 20 2a 2f 0a 20 20 70 20 3d 20 28 50 65 72 63 65   */.  p = (Perce
1050: 6e 74 69 6c 65 2a 29 73 71 6c 69 74 65 33 5f 61  ntile*)sqlite3_a
1060: 67 67 72 65 67 61 74 65 5f 63 6f 6e 74 65 78 74  ggregate_context
1070: 28 70 43 74 78 2c 20 73 69 7a 65 6f 66 28 2a 70  (pCtx, sizeof(*p
1080: 29 29 3b 0a 20 20 69 66 28 20 70 3d 3d 30 20 29  ));.  if( p==0 )
1090: 20 72 65 74 75 72 6e 3b 0a 0a 20 20 2f 2a 20 52   return;..  /* R
10a0: 65 6d 65 6d 62 65 72 20 74 68 65 20 50 20 76 61  emember the P va
10b0: 6c 75 65 2e 20 20 54 68 72 6f 77 20 61 6e 20 65  lue.  Throw an e
10c0: 72 72 6f 72 20 69 66 20 74 68 65 20 50 20 76 61  rror if the P va
10d0: 6c 75 65 20 69 73 20 64 69 66 66 65 72 65 6e 74  lue is different
10e0: 0a 20 20 2a 2a 20 66 72 6f 6d 20 61 6e 79 20 70  .  ** from any p
10f0: 72 69 6f 72 20 72 6f 77 2c 20 70 65 72 20 52 65  rior row, per Re
1100: 71 75 69 72 65 6d 65 6e 74 20 28 32 29 2e 20 2a  quirement (2). *
1110: 2f 0a 20 20 69 66 28 20 70 2d 3e 72 50 63 74 3d  /.  if( p->rPct=
1120: 3d 30 2e 30 20 29 7b 0a 20 20 20 20 70 2d 3e 72  =0.0 ){.    p->r
1130: 50 63 74 20 3d 20 72 50 63 74 2b 31 2e 30 3b 0a  Pct = rPct+1.0;.
1140: 20 20 7d 65 6c 73 65 20 69 66 28 20 21 73 61 6d    }else if( !sam
1150: 65 56 61 6c 75 65 28 70 2d 3e 72 50 63 74 2c 72  eValue(p->rPct,r
1160: 50 63 74 2b 31 2e 30 29 20 29 7b 0a 20 20 20 20  Pct+1.0) ){.    
1170: 73 71 6c 69 74 65 33 5f 72 65 73 75 6c 74 5f 65  sqlite3_result_e
1180: 72 72 6f 72 28 70 43 74 78 2c 20 22 32 6e 64 20  rror(pCtx, "2nd 
1190: 61 72 67 75 6d 65 6e 74 20 74 6f 20 70 65 72 63  argument to perc
11a0: 65 6e 74 69 6c 65 28 29 20 69 73 20 6e 6f 74 20  entile() is not 
11b0: 74 68 65 20 22 0a 20 20 20 20 20 20 20 20 20 20  the ".          
11c0: 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20                  
11d0: 20 20 20 20 20 22 73 61 6d 65 20 66 6f 72 20 61       "same for a
11e0: 6c 6c 20 69 6e 70 75 74 20 72 6f 77 73 22 2c 20  ll input rows", 
11f0: 2d 31 29 3b 0a 20 20 20 20 72 65 74 75 72 6e 3b  -1);.    return;
1200: 0a 20 20 7d 0a 0a 20 20 2f 2a 20 49 67 6e 6f 72  .  }..  /* Ignor
1210: 65 20 72 6f 77 73 20 66 6f 72 20 77 68 69 63 68  e rows for which
1220: 20 59 20 69 73 20 4e 55 4c 4c 20 2a 2f 0a 20 20   Y is NULL */.  
1230: 65 54 79 70 65 20 3d 20 73 71 6c 69 74 65 33 5f  eType = sqlite3_
1240: 76 61 6c 75 65 5f 74 79 70 65 28 61 72 67 76 5b  value_type(argv[
1250: 30 5d 29 3b 0a 20 20 69 66 28 20 65 54 79 70 65  0]);.  if( eType
1260: 3d 3d 53 51 4c 49 54 45 5f 4e 55 4c 4c 20 29 20  ==SQLITE_NULL ) 
1270: 72 65 74 75 72 6e 3b 0a 0a 20 20 2f 2a 20 49 66  return;..  /* If
1280: 20 6e 6f 74 20 4e 55 4c 4c 2c 20 74 68 65 6e 20   not NULL, then 
1290: 59 20 6d 75 73 74 20 62 65 20 6e 75 6d 65 72 69  Y must be numeri
12a0: 63 2e 20 20 4f 74 68 65 72 77 69 73 65 20 74 68  c.  Otherwise th
12b0: 72 6f 77 20 61 6e 20 65 72 72 6f 72 2e 0a 20 20  row an error..  
12c0: 2a 2a 20 52 65 71 75 69 72 65 6d 65 6e 74 20 34  ** Requirement 4
12d0: 20 2a 2f 0a 20 20 69 66 28 20 65 54 79 70 65 21   */.  if( eType!
12e0: 3d 53 51 4c 49 54 45 5f 49 4e 54 45 47 45 52 20  =SQLITE_INTEGER 
12f0: 26 26 20 65 54 79 70 65 21 3d 53 51 4c 49 54 45  && eType!=SQLITE
1300: 5f 46 4c 4f 41 54 20 29 7b 0a 20 20 20 20 73 71  _FLOAT ){.    sq
1310: 6c 69 74 65 33 5f 72 65 73 75 6c 74 5f 65 72 72  lite3_result_err
1320: 6f 72 28 70 43 74 78 2c 20 22 31 73 74 20 61 72  or(pCtx, "1st ar
1330: 67 75 6d 65 6e 74 20 74 6f 20 70 65 72 63 65 6e  gument to percen
1340: 74 69 6c 65 28 29 20 69 73 20 6e 6f 74 20 22 0a  tile() is not ".
1350: 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20                  
1360: 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 22                 "
1370: 6e 75 6d 65 72 69 63 22 2c 20 2d 31 29 3b 0a 20  numeric", -1);. 
1380: 20 20 20 72 65 74 75 72 6e 3b 0a 20 20 7d 0a 0a     return;.  }..
1390: 20 20 2f 2a 20 54 68 72 6f 77 20 61 6e 20 65 72    /* Throw an er
13a0: 72 6f 72 20 69 66 20 74 68 65 20 59 20 76 61 6c  ror if the Y val
13b0: 75 65 20 69 73 20 69 6e 66 69 6e 69 74 79 20 6f  ue is infinity o
13c0: 72 20 4e 61 4e 20 2a 2f 0a 20 20 79 20 3d 20 73  r NaN */.  y = s
13d0: 71 6c 69 74 65 33 5f 76 61 6c 75 65 5f 64 6f 75  qlite3_value_dou
13e0: 62 6c 65 28 61 72 67 76 5b 30 5d 29 3b 0a 20 20  ble(argv[0]);.  
13f0: 69 66 28 20 69 73 49 6e 66 69 6e 69 74 79 28 79  if( isInfinity(y
1400: 29 20 29 7b 0a 20 20 20 20 73 71 6c 69 74 65 33  ) ){.    sqlite3
1410: 5f 72 65 73 75 6c 74 5f 65 72 72 6f 72 28 70 43  _result_error(pC
1420: 74 78 2c 20 22 49 6e 66 20 69 6e 70 75 74 20 74  tx, "Inf input t
1430: 6f 20 70 65 72 63 65 6e 74 69 6c 65 28 29 22 2c  o percentile()",
1440: 20 2d 31 29 3b 0a 20 20 20 20 72 65 74 75 72 6e   -1);.    return
1450: 3b 0a 20 20 7d 0a 0a 20 20 2f 2a 20 41 6c 6c 6f  ;.  }..  /* Allo
1460: 63 61 74 65 20 61 6e 64 20 73 74 6f 72 65 20 74  cate and store t
1470: 68 65 20 59 20 2a 2f 0a 20 20 69 66 28 20 70 2d  he Y */.  if( p-
1480: 3e 6e 55 73 65 64 3e 3d 70 2d 3e 6e 41 6c 6c 6f  >nUsed>=p->nAllo
1490: 63 20 29 7b 0a 20 20 20 20 75 6e 73 69 67 6e 65  c ){.    unsigne
14a0: 64 20 6e 20 3d 20 70 2d 3e 6e 41 6c 6c 6f 63 2a  d n = p->nAlloc*
14b0: 32 20 2b 20 32 35 30 3b 0a 20 20 20 20 64 6f 75  2 + 250;.    dou
14c0: 62 6c 65 20 2a 61 20 3d 20 73 71 6c 69 74 65 33  ble *a = sqlite3
14d0: 5f 72 65 61 6c 6c 6f 63 36 34 28 70 2d 3e 61 2c  _realloc64(p->a,
14e0: 20 73 69 7a 65 6f 66 28 64 6f 75 62 6c 65 29 2a   sizeof(double)*
14f0: 6e 29 3b 0a 20 20 20 20 69 66 28 20 61 3d 3d 30  n);.    if( a==0
1500: 20 29 7b 0a 20 20 20 20 20 20 73 71 6c 69 74 65   ){.      sqlite
1510: 33 5f 66 72 65 65 28 70 2d 3e 61 29 3b 0a 20 20  3_free(p->a);.  
1520: 20 20 20 20 6d 65 6d 73 65 74 28 70 2c 20 30 2c      memset(p, 0,
1530: 20 73 69 7a 65 6f 66 28 2a 70 29 29 3b 0a 20 20   sizeof(*p));.  
1540: 20 20 20 20 73 71 6c 69 74 65 33 5f 72 65 73 75      sqlite3_resu
1550: 6c 74 5f 65 72 72 6f 72 5f 6e 6f 6d 65 6d 28 70  lt_error_nomem(p
1560: 43 74 78 29 3b 0a 20 20 20 20 20 20 72 65 74 75  Ctx);.      retu
1570: 72 6e 3b 0a 20 20 20 20 7d 0a 20 20 20 20 70 2d  rn;.    }.    p-
1580: 3e 6e 41 6c 6c 6f 63 20 3d 20 6e 3b 0a 20 20 20  >nAlloc = n;.   
1590: 20 70 2d 3e 61 20 3d 20 61 3b 0a 20 20 7d 0a 20   p->a = a;.  }. 
15a0: 20 70 2d 3e 61 5b 70 2d 3e 6e 55 73 65 64 2b 2b   p->a[p->nUsed++
15b0: 5d 20 3d 20 79 3b 0a 7d 0a 0a 2f 2a 0a 2a 2a 20  ] = y;.}../*.** 
15c0: 43 6f 6d 70 61 72 65 20 74 6f 20 64 6f 75 62 6c  Compare to doubl
15d0: 65 73 20 66 6f 72 20 73 6f 72 74 69 6e 67 20 75  es for sorting u
15e0: 73 69 6e 67 20 71 73 6f 72 74 28 29 0a 2a 2f 0a  sing qsort().*/.
15f0: 73 74 61 74 69 63 20 69 6e 74 20 53 51 4c 49 54  static int SQLIT
1600: 45 5f 43 44 45 43 4c 20 64 6f 75 62 6c 65 43 6d  E_CDECL doubleCm
1610: 70 28 63 6f 6e 73 74 20 76 6f 69 64 20 2a 70 41  p(const void *pA
1620: 2c 20 63 6f 6e 73 74 20 76 6f 69 64 20 2a 70 42  , const void *pB
1630: 29 7b 0a 20 20 64 6f 75 62 6c 65 20 61 20 3d 20  ){.  double a = 
1640: 2a 28 64 6f 75 62 6c 65 2a 29 70 41 3b 0a 20 20  *(double*)pA;.  
1650: 64 6f 75 62 6c 65 20 62 20 3d 20 2a 28 64 6f 75  double b = *(dou
1660: 62 6c 65 2a 29 70 42 3b 0a 20 20 69 66 28 20 61  ble*)pB;.  if( a
1670: 3d 3d 62 20 29 20 72 65 74 75 72 6e 20 30 3b 0a  ==b ) return 0;.
1680: 20 20 69 66 28 20 61 3c 62 20 29 20 72 65 74 75    if( a<b ) retu
1690: 72 6e 20 2d 31 3b 0a 20 20 72 65 74 75 72 6e 20  rn -1;.  return 
16a0: 2b 31 3b 0a 7d 0a 0a 2f 2a 0a 2a 2a 20 43 61 6c  +1;.}../*.** Cal
16b0: 6c 65 64 20 74 6f 20 63 6f 6d 70 75 74 65 20 74  led to compute t
16c0: 68 65 20 66 69 6e 61 6c 20 6f 75 74 70 75 74 20  he final output 
16d0: 6f 66 20 70 65 72 63 65 6e 74 69 6c 65 28 29 20  of percentile() 
16e0: 61 6e 64 20 74 6f 20 63 6c 65 61 6e 0a 2a 2a 20  and to clean.** 
16f0: 75 70 20 61 6c 6c 20 61 6c 6c 6f 63 61 74 65 64  up all allocated
1700: 20 6d 65 6d 6f 72 79 2e 0a 2a 2f 0a 73 74 61 74   memory..*/.stat
1710: 69 63 20 76 6f 69 64 20 70 65 72 63 65 6e 74 46  ic void percentF
1720: 69 6e 61 6c 28 73 71 6c 69 74 65 33 5f 63 6f 6e  inal(sqlite3_con
1730: 74 65 78 74 20 2a 70 43 74 78 29 7b 0a 20 20 50  text *pCtx){.  P
1740: 65 72 63 65 6e 74 69 6c 65 20 2a 70 3b 0a 20 20  ercentile *p;.  
1750: 75 6e 73 69 67 6e 65 64 20 69 31 2c 20 69 32 3b  unsigned i1, i2;
1760: 0a 20 20 64 6f 75 62 6c 65 20 76 31 2c 20 76 32  .  double v1, v2
1770: 3b 0a 20 20 64 6f 75 62 6c 65 20 69 78 2c 20 76  ;.  double ix, v
1780: 78 3b 0a 20 20 70 20 3d 20 28 50 65 72 63 65 6e  x;.  p = (Percen
1790: 74 69 6c 65 2a 29 73 71 6c 69 74 65 33 5f 61 67  tile*)sqlite3_ag
17a0: 67 72 65 67 61 74 65 5f 63 6f 6e 74 65 78 74 28  gregate_context(
17b0: 70 43 74 78 2c 20 30 29 3b 0a 20 20 69 66 28 20  pCtx, 0);.  if( 
17c0: 70 3d 3d 30 20 29 20 72 65 74 75 72 6e 3b 0a 20  p==0 ) return;. 
17d0: 20 69 66 28 20 70 2d 3e 61 3d 3d 30 20 29 20 72   if( p->a==0 ) r
17e0: 65 74 75 72 6e 3b 0a 20 20 69 66 28 20 70 2d 3e  eturn;.  if( p->
17f0: 6e 55 73 65 64 20 29 7b 0a 20 20 20 20 71 73 6f  nUsed ){.    qso
1800: 72 74 28 70 2d 3e 61 2c 20 70 2d 3e 6e 55 73 65  rt(p->a, p->nUse
1810: 64 2c 20 73 69 7a 65 6f 66 28 64 6f 75 62 6c 65  d, sizeof(double
1820: 29 2c 20 64 6f 75 62 6c 65 43 6d 70 29 3b 0a 20  ), doubleCmp);. 
1830: 20 20 20 69 78 20 3d 20 28 70 2d 3e 72 50 63 74     ix = (p->rPct
1840: 2d 31 2e 30 29 2a 28 70 2d 3e 6e 55 73 65 64 2d  -1.0)*(p->nUsed-
1850: 31 29 2a 30 2e 30 31 3b 0a 20 20 20 20 69 31 20  1)*0.01;.    i1 
1860: 3d 20 28 75 6e 73 69 67 6e 65 64 29 69 78 3b 0a  = (unsigned)ix;.
1870: 20 20 20 20 69 32 20 3d 20 69 78 3d 3d 28 64 6f      i2 = ix==(do
1880: 75 62 6c 65 29 69 31 20 7c 7c 20 69 31 3d 3d 70  uble)i1 || i1==p
1890: 2d 3e 6e 55 73 65 64 2d 31 20 3f 20 69 31 20 3a  ->nUsed-1 ? i1 :
18a0: 20 69 31 2b 31 3b 0a 20 20 20 20 76 31 20 3d 20   i1+1;.    v1 = 
18b0: 70 2d 3e 61 5b 69 31 5d 3b 0a 20 20 20 20 76 32  p->a[i1];.    v2
18c0: 20 3d 20 70 2d 3e 61 5b 69 32 5d 3b 0a 20 20 20   = p->a[i2];.   
18d0: 20 76 78 20 3d 20 76 31 20 2b 20 28 76 32 2d 76   vx = v1 + (v2-v
18e0: 31 29 2a 28 69 78 2d 69 31 29 3b 0a 20 20 20 20  1)*(ix-i1);.    
18f0: 73 71 6c 69 74 65 33 5f 72 65 73 75 6c 74 5f 64  sqlite3_result_d
1900: 6f 75 62 6c 65 28 70 43 74 78 2c 20 76 78 29 3b  ouble(pCtx, vx);
1910: 0a 20 20 7d 0a 20 20 73 71 6c 69 74 65 33 5f 66  .  }.  sqlite3_f
1920: 72 65 65 28 70 2d 3e 61 29 3b 0a 20 20 6d 65 6d  ree(p->a);.  mem
1930: 73 65 74 28 70 2c 20 30 2c 20 73 69 7a 65 6f 66  set(p, 0, sizeof
1940: 28 2a 70 29 29 3b 0a 7d 0a 0a 0a 23 69 66 64 65  (*p));.}...#ifde
1950: 66 20 5f 57 49 4e 33 32 0a 5f 5f 64 65 63 6c 73  f _WIN32.__decls
1960: 70 65 63 28 64 6c 6c 65 78 70 6f 72 74 29 0a 23  pec(dllexport).#
1970: 65 6e 64 69 66 0a 69 6e 74 20 73 71 6c 69 74 65  endif.int sqlite
1980: 33 5f 70 65 72 63 65 6e 74 69 6c 65 5f 69 6e 69  3_percentile_ini
1990: 74 28 0a 20 20 73 71 6c 69 74 65 33 20 2a 64 62  t(.  sqlite3 *db
19a0: 2c 20 0a 20 20 63 68 61 72 20 2a 2a 70 7a 45 72  , .  char **pzEr
19b0: 72 4d 73 67 2c 20 0a 20 20 63 6f 6e 73 74 20 73  rMsg, .  const s
19c0: 71 6c 69 74 65 33 5f 61 70 69 5f 72 6f 75 74 69  qlite3_api_routi
19d0: 6e 65 73 20 2a 70 41 70 69 0a 29 7b 0a 20 20 69  nes *pApi.){.  i
19e0: 6e 74 20 72 63 20 3d 20 53 51 4c 49 54 45 5f 4f  nt rc = SQLITE_O
19f0: 4b 3b 0a 20 20 53 51 4c 49 54 45 5f 45 58 54 45  K;.  SQLITE_EXTE
1a00: 4e 53 49 4f 4e 5f 49 4e 49 54 32 28 70 41 70 69  NSION_INIT2(pApi
1a10: 29 3b 0a 20 20 28 76 6f 69 64 29 70 7a 45 72 72  );.  (void)pzErr
1a20: 4d 73 67 3b 20 20 2f 2a 20 55 6e 75 73 65 64 20  Msg;  /* Unused 
1a30: 70 61 72 61 6d 65 74 65 72 20 2a 2f 0a 20 20 72  parameter */.  r
1a40: 63 20 3d 20 73 71 6c 69 74 65 33 5f 63 72 65 61  c = sqlite3_crea
1a50: 74 65 5f 66 75 6e 63 74 69 6f 6e 28 64 62 2c 20  te_function(db, 
1a60: 22 70 65 72 63 65 6e 74 69 6c 65 22 2c 20 32 2c  "percentile", 2,
1a70: 20 53 51 4c 49 54 45 5f 55 54 46 38 2c 20 30 2c   SQLITE_UTF8, 0,
1a80: 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20  .               
1a90: 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20                  
1aa0: 30 2c 20 70 65 72 63 65 6e 74 53 74 65 70 2c 20  0, percentStep, 
1ab0: 70 65 72 63 65 6e 74 46 69 6e 61 6c 29 3b 0a 20  percentFinal);. 
1ac0: 20 72 65 74 75 72 6e 20 72 63 3b 0a 7d 0a         return rc;.}.