Introduce signal_shutdown().
[paraslash.git] / mood.c
1 /*
2  * Copyright (C) 2007 Andre Noll <maan@systemlinux.org>
3  *
4  * Licensed under the GPL v2. For licencing details see COPYING.
5  */
6
7 /** \file mood.c Paraslash's mood handling functions. */
8
9 #include <fnmatch.h>
10 #include "para.h"
11 #include "error.h"
12 #include "string.h"
13 #include "afh.h"
14 #include "afs.h"
15 #include "list.h"
16
17 /**
18  * Contains statistical data of the currently admissible audio files.
19  *
20  * It is used to assign normalized score values to each admissible audio file.
21  */
22 struct afs_statistics {
23         /** Sum of num played over all admissible files. */
24         int64_t num_played_sum;
25         /** Sum of last played times over all admissible files. */
26         int64_t last_played_sum;
27         /** Quadratic deviation of num played time. */
28         int64_t num_played_qd;
29         /** Quadratic deviation of last played time. */
30         int64_t last_played_qd;
31         /** Number of admissible files */
32         unsigned num;
33 };
34 struct afs_statistics statistics;
35
36 /**
37  * Assign scores according to a mood_method.
38  *
39  * Each mood_method has its own mood_score_function. The first three parameters
40  * passed to that function are informations about the audio file whose score is
41  * to be computed. The data argument depends on the mood method this function
42  * is used for. It usually is the argument given at the end of a mood line.
43  *
44  * Mood score functions must return values between -100 and +100 inclusively.
45  * Boolean score functions should always return either -100 or +100.
46  *
47  * \sa struct mood_method, mood_parser.
48  */
49 typedef int mood_score_function(const char *path, const struct afs_info *afsi,
50                 const struct afh_info *afhi, const void *data);
51
52 /**
53  * Pre-process a mood line.
54  *
55  * The mood_parser of a mood_method is called once at mood open time for each
56  * line of the current mood definition that contains the mood_method's name as
57  * a keyword. The line is passed to the mood_parser as the first argument. The
58  * mood_parser must determine whether the line is syntactically correct and
59  * return a positive value if so and a negative value otherwise.
60  *
61  * Some mood parsers pre-process the data given in the mood line to compute a
62  * structure which depends of the particular mood_method and which is used
63  * later in the mood_score_function of the mood_method. The mood_parser may
64  * store a pointer to its structure via the second argument.
65  *
66  * \sa mood_open(), mood_cleanup_function, mood_score_function.
67  */
68 typedef int mood_parser(const char *, void **);
69
70 /**
71  * Deallocate resources which were allocated by the mood_parser.
72  *
73  * This optional function of a mood_method is used to free any resources
74  * allocated in mood_open() by the mood_parser. The argument passed is a
75  * pointer to the mood_method specific data structure that was returned by the
76  * mood_parser.
77  *
78  * \sa mood_parser.
79  */
80 typedef void mood_cleanup_function(void *);
81
82 /**
83  * Used for scoring and to determine whether a file is admissible.
84  */
85 struct mood_method {
86         /* The name of the method. */
87         const char *name;
88         /** Pointer to the mood parser. */
89         mood_parser *parser;
90         /** Pointer to the score function */
91         mood_score_function *score_function;
92         /** Optional cleanup function. */
93         mood_cleanup_function *cleanup;
94 };
95
96 /**
97  * Each line of the current mood corresponds to a mood_item.
98  */
99 struct mood_item {
100         /** The method this line is referring to. */
101         const struct mood_method *method;
102         /** The data structure computed by the mood parser. */
103         void *parser_data;
104         /** The given score value, or zero if none was given. */
105         int32_t score_arg;
106         /** Non-zero if random scoring was requested. */
107         int random_score;
108         /** Whether the "not" keyword was given in the mood line. */
109         int logical_not;
110         /** The position in the list of items. */
111         struct list_head mood_item_node;
112 };
113
114 /**
115  * Created from the mood definition by mood_open().
116  *
117  * When a mood is opened, each line of its definition is investigated, and a
118  * corresponding mood item is produced. Each mood line starts with \p accept,
119  * \p deny, or \p score which determines the type of the mood line.  For each
120  * such type a linked list is maintained whose entries are the mood items.
121  *
122  * \sa mood_item, mood_open().
123  */
124 struct mood {
125         /** The name of this mood. */
126         char *name;
127         /** The list of mood items of type \p accept. */
128         struct list_head accept_list;
129         /** The list of mood items of type \p deny. */
130         struct list_head deny_list;
131         /** The list of mood items of type \p score. */
132         struct list_head score_list;
133 };
134
135 static struct mood *current_mood;
136
137 /**
138  *  Rough approximation to sqrt.
139  *
140  *  \param x Integer of which to calculate the sqrt.
141  *
142  *  \return An integer res with res * res <= x.
143  */
144 static uint64_t int_sqrt(uint64_t x)
145 {
146         uint64_t op, res, one = 1;
147         op = x;
148         res = 0;
149
150         one = one << 62;
151         while (one > op)
152                 one >>= 2;
153
154         while (one != 0) {
155                 if (op >= res + one) {
156                         op = op - (res + one);
157                         res = res +  2 * one;
158                 }
159                 res /= 2;
160                 one /= 4;
161         }
162 //      PARA_NOTICE_LOG("sqrt(%llu) = %llu\n", x, res);
163         return res;
164 }
165
166 static int mm_no_attributes_set_parser(const char *arg, __a_unused void **ignored)
167 {
168         if (arg && *arg)
169                 PARA_WARNING_LOG("ignored junk at eol: %s\n", arg);
170         return 1;
171 }
172
173 static int mm_no_attributes_set_score_function(__a_unused const char *path,
174                 const struct afs_info *afsi,
175                 __a_unused const struct afh_info *afhi,
176                 __a_unused const void *data)
177 {
178         if (!afsi->attributes)
179                 return 100;
180         return -100;
181 }
182
183 static int mm_played_rarely_score_function(__a_unused const char *path,
184                 const struct afs_info *afsi,
185                 __a_unused const struct afh_info *afhi,
186                 __a_unused const void *data)
187 {
188         unsigned num;
189         int ret = get_num_admissible_files(&num);
190
191         if (ret < 0)
192                 return 0;
193         if (statistics.num_played_sum - num * afsi->num_played
194                         > int_sqrt(statistics.num_played_qd * num))
195                 return 100;
196         return -100;
197 }
198
199 static int mm_played_rarely_parser(const char *arg, __a_unused void **ignored)
200 {
201         if (arg && *arg)
202                 PARA_WARNING_LOG("ignored junk at eol: %s\n", arg);
203         return 1;
204 }
205
206 static int mm_name_like_score_function(const char *path,
207                 __a_unused const struct afs_info *afsi,
208                 __a_unused const struct afh_info *afhi,
209                 const void *data)
210 {
211         if (fnmatch(data, path, 0))
212                 return -100;
213         return 100;
214 }
215
216 static int mm_name_like_parser(const char *arg, void **data)
217 {
218         *data = para_strdup(arg);
219         return 1;
220 }
221
222 static void mm_name_like_cleanup(void *data)
223 {
224         free(data);
225 }
226
227 static int mm_is_set_parser(const char *arg, void **bitnum)
228 {
229         unsigned char *c = para_malloc(1);
230         int ret = get_attribute_bitnum_by_name(arg, c);
231
232         if (ret >= 0)
233                 *bitnum = c;
234         else
235                 free(c);
236         return ret;
237 }
238
239 static int mm_is_set_score_function(__a_unused const char *path,
240                 __a_unused const struct afs_info *afsi,
241                 __a_unused const struct afh_info *afhi,
242                 const void *data)
243 {
244         const unsigned char *bn = data;
245         if (afsi->attributes & (1ULL << *bn))
246                 return 100;
247         return -100;
248 }
249
250 /* returns 1 if row matches score item, negative otherwise */
251 static int add_item_score(const struct osl_row *row, struct mood_item *item, long *score,
252                 long *score_arg_sum)
253 {
254         struct afs_info afsi;
255         struct afh_info afhi;
256         char *path;
257         int ret;
258
259         *score_arg_sum += item->random_score? 100 : PARA_ABS(item->score_arg);
260         ret = 100;
261         if (item->method) {
262                 ret = get_afsi_of_row(row, &afsi);
263                 if (ret< 0)
264                         return ret;
265                 ret = get_afhi_of_row(row, &afhi);
266                 if (ret< 0)
267                         return ret;
268                 ret = get_audio_file_path_of_row(row, &path);
269                 if (ret< 0)
270                         return ret;
271                 ret = item->method->score_function(path, &afsi, &afhi,
272                         item->parser_data);
273                 if ((ret < 0 && !item->logical_not) || (ret >= 0 && item->logical_not))
274                         return -1; /* no match */
275         }
276         if (item->random_score)
277                 *score += PARA_ABS(ret) * para_random(100);
278         else
279                 *score += PARA_ABS(ret) * item->score_arg;
280         return 1;
281 }
282
283 static int compute_mood_score(const struct osl_row *aft_row, struct mood *m,
284                 long *result)
285 {
286         struct mood_item *item;
287         int match = 0;
288         long score_arg_sum = 0, score = 0;
289
290         if (!m)
291                 return -E_NO_MOOD;
292         /* reject audio file if it matches any entry in the deny list */
293         list_for_each_entry(item, &m->deny_list, mood_item_node)
294                 if (add_item_score(aft_row, item, &score, &score_arg_sum) > 0)
295                         return -E_NOT_ADMISSIBLE;
296         list_for_each_entry(item, &m->accept_list, mood_item_node)
297                 if (add_item_score(aft_row, item, &score, &score_arg_sum) > 0)
298                         match = 1;
299         /* reject if there is no matching entry in the accept list */
300         if (!match && !list_empty(&m->accept_list))
301                 return -E_NOT_ADMISSIBLE;
302         list_for_each_entry(item, &m->score_list, mood_item_node)
303                 add_item_score(aft_row, item, &score, &score_arg_sum);
304         if (score_arg_sum)
305                 score /= score_arg_sum;
306         *result = score;
307         return 1;
308 }
309
310 #define DEFINE_MOOD_METHOD(_name) \
311 .parser = mm_ ## _name ## _parser, \
312 .score_function = mm_ ## _name ## _score_function, \
313 .name = #_name
314
315 #define DEFINE_MOOD_METHOD_WITH_CLEANUP(_name) \
316         DEFINE_MOOD_METHOD(_name), \
317         .cleanup = mm_ ## _name ## _cleanup
318
319 static const struct mood_method mood_methods[] = {
320         {DEFINE_MOOD_METHOD(no_attributes_set)},
321         {DEFINE_MOOD_METHOD(played_rarely)},
322         {DEFINE_MOOD_METHOD(is_set)},
323         {DEFINE_MOOD_METHOD_WITH_CLEANUP(name_like)},
324         {.parser = NULL}
325 };
326
327 static void cleanup_list_entry(struct mood_item *item)
328 {
329         if (item->method && item->method->cleanup)
330                 item->method->cleanup(item->parser_data);
331         else
332                 free(item->parser_data);
333         list_del(&item->mood_item_node);
334         free(item);
335 }
336
337 static void destroy_mood(struct mood *m)
338 {
339         struct mood_item *tmp, *item;
340
341         if (!m)
342                 return;
343         list_for_each_entry_safe(item, tmp, &m->accept_list, mood_item_node)
344                 cleanup_list_entry(item);
345         list_for_each_entry_safe(item, tmp, &m->deny_list, mood_item_node)
346                 cleanup_list_entry(item);
347         list_for_each_entry_safe(item, tmp, &m->score_list, mood_item_node)
348                 cleanup_list_entry(item);
349         free(m->name);
350         free(m);
351 }
352
353 static struct mood *alloc_new_mood(const char *name)
354 {
355         struct mood *m = para_calloc(sizeof(struct mood));
356         m->name = para_strdup(name);
357         INIT_LIST_HEAD(&m->accept_list);
358         INIT_LIST_HEAD(&m->deny_list);
359         INIT_LIST_HEAD(&m->score_list);
360         return m;
361 }
362
363 /** The different types of a mood line. */
364 enum mood_line_type {
365         /** Invalid. */
366         ML_INVALID,
367         /** Accept line. */
368         ML_ACCEPT,
369         /** Deny line. */
370         ML_DENY,
371         /** Score line. */
372         ML_SCORE
373 };
374
375 struct mood_line_parser_data {
376         struct mood *m;
377         unsigned line_num;
378 };
379
380 /*
381  * <accept [with score <score>] | deny [with score <score>]  | score <score>>
382  *      [if] [not] <mood_method> [options]
383  * <score> is either an integer or "random" which assigns a random score to
384  * all matching files
385  */
386
387 static int parse_mood_line(char *mood_line, void *data)
388 {
389         struct mood_line_parser_data *mlpd = data;
390         char **argv;
391         char *delim = " \t";
392         unsigned num_words;
393         char **w;
394         int i, ret;
395         enum mood_line_type mlt = ML_INVALID;
396         struct mood_item *mi = NULL;
397         char *buf = para_strdup(mood_line);
398
399         mlpd->line_num++;
400         num_words = split_args(buf, &argv, delim);
401         ret = 1;
402         if (!num_words) /* empty line */
403                 goto out;
404         w = argv;
405         if (**w == '#') /* comment */
406                 goto out;
407         if (!strcmp(*w, "accept"))
408                 mlt = ML_ACCEPT;
409         else if (!strcmp(*w, "deny"))
410                 mlt = ML_DENY;
411         else if (!strcmp(*w, "score"))
412                 mlt = ML_SCORE;
413         ret = -E_MOOD_SYNTAX;
414         if (mlt == ML_INVALID)
415                 goto out;
416         mi = para_calloc(sizeof(struct mood_item));
417         if (mlt != ML_SCORE) {
418                 ret = -E_MOOD_SYNTAX;
419                 w++;
420                 if (!*w)
421                         goto out;
422                 if (strcmp(*w, "with"))
423                         goto check_for_if;
424                 w++;
425                 if (!*w)
426                         goto out;
427                 if (strcmp(*w, "score"))
428                         goto out;
429         }
430         if (mlt == ML_SCORE || !strcmp(*w, "score")) {
431                 ret = -E_MOOD_SYNTAX;
432                 w++;
433                 if (!*w)
434                         goto out;
435                 if (strcmp(*w, "random")) {
436                         mi->random_score = 0;
437                         ret = para_atoi32(*w, &mi->score_arg);
438                         if (ret < 0)
439                                 goto out;
440                 } else {
441                         mi->random_score = 1;
442                         if (!*(w + 1))
443                         goto success; /* the line "score random" is valid */
444                 }
445         } else
446                 mi->score_arg = 0;
447         ret = -E_MOOD_SYNTAX;
448         w++;
449         if (!*w)
450                 goto out;
451 check_for_if:
452         if (!strcmp(*w, "if")) {
453                 ret = -E_MOOD_SYNTAX;
454                 w++;
455                 if (!*w)
456                         goto out;
457         }
458         if (!strcmp(*w, "not")) {
459                 ret = -E_MOOD_SYNTAX;
460                 w++;
461                 if (!*w)
462                         goto out;
463                 mi->logical_not = 1;
464         } else
465                 mi->logical_not = 0;
466         for (i = 0; mood_methods[i].parser; i++) {
467                 if (strcmp(*w, mood_methods[i].name))
468                         continue;
469                 break;
470         }
471         ret = -E_MOOD_SYNTAX;
472         if (!mood_methods[i].parser)
473                 goto out;
474         w++;
475         ret = mood_methods[i].parser(*w, &mi->parser_data);
476         if (ret < 0)
477                 goto out;
478         mi->method = &mood_methods[i];
479 success:
480         if (mlpd->m) {
481                 if (mlt == ML_ACCEPT)
482                         para_list_add(&mi->mood_item_node, &mlpd->m->accept_list);
483                 else if (mlt == ML_DENY)
484                         para_list_add(&mi->mood_item_node, &mlpd->m->deny_list);
485                 else
486                         para_list_add(&mi->mood_item_node, &mlpd->m->score_list);
487         }
488         PARA_DEBUG_LOG("%s entry added, method: %p\n", mlt == ML_ACCEPT? "accept" :
489                 (mlt == ML_DENY? "deny" : "score"), mi->method);
490         ret = 1;
491 out:
492         free(argv);
493         free(buf);
494         if (ret >= 0)
495                 return ret;
496         if (mi) {
497                 free(mi->parser_data);
498                 free(mi);
499         }
500         return ret;
501 }
502
503 static int load_mood(const struct osl_row *mood_row, struct mood **m)
504 {
505         char *mood_name;
506         struct osl_object mood_def;
507         struct mood_line_parser_data mlpd = {.line_num = 0};
508         int ret = mood_get_name_and_def_by_row(mood_row, &mood_name, &mood_def);
509
510         if (ret < 0)
511                 return ret;
512         if (!*mood_name)
513                 return -E_DUMMY_ROW;
514         mlpd.m = alloc_new_mood(mood_name);
515         ret = for_each_line_ro(mood_def.data, mood_def.size,
516                 parse_mood_line, &mlpd);
517         osl_close_disk_object(&mood_def);
518         if (ret < 0) {
519                 PARA_ERROR_LOG("unable to load mood %s: %s\n", mlpd.m->name,
520                         PARA_STRERROR(-ret));
521                 destroy_mood(mlpd.m);
522                 return ret;
523         }
524         *m = mlpd.m;
525         return 1;
526 }
527
528 static int check_mood(struct osl_row *mood_row, void *data)
529 {
530         struct para_buffer *pb = data;
531         char *mood_name;
532         struct osl_object mood_def;
533         struct mood_line_parser_data mlpd = {.line_num = 0};
534
535         int ret = mood_get_name_and_def_by_row(mood_row, &mood_name, &mood_def);
536
537         if (ret < 0) {
538                 para_printf(pb, "failed to get mood definition\n");
539                 return ret;
540         }
541         if (!*mood_name) /* ignore dummy row */
542                 goto out;
543         para_printf(pb, "checking mood %s...\n", mood_name);
544         ret = for_each_line_ro(mood_def.data, mood_def.size,
545                 parse_mood_line, &mlpd);
546         if (ret < 0)
547                 para_printf(pb, "%s line %u: %s\n", mood_name, mlpd.line_num,
548                         PARA_STRERROR(-ret));
549 out:
550         osl_close_disk_object(&mood_def);
551         return 1;
552 }
553
554 /**
555  * Check all moods for syntax errors.
556  *
557  * \param query Unused.
558  * \param result: Contains check messages.
559  */
560 int mood_check_callback(__a_unused const struct osl_object *query,
561         struct osl_object *result)
562 {
563         struct para_buffer pb = {.buf = NULL};
564
565         para_printf(&pb, "checking moods...\n");
566         osl_rbtree_loop(moods_table, BLOBCOL_ID, &pb,
567                 check_mood);
568         result->data = pb.buf;
569         result->size = pb.size;
570         return 1;
571 }
572
573 #if 0
574 static unsigned int_log2(uint64_t x)
575 {
576         unsigned res = 0;
577
578         while (x) {
579                 x /= 2;
580                 res++;
581         }
582         return res;
583 }
584 #endif
585
586 static int64_t normalized_value(int64_t x, int64_t n, int64_t sum, int64_t qd)
587 {
588         if (!n || !qd)
589                 return 0;
590         return 100 * (n * x - sum) / (int64_t)int_sqrt(n * qd);
591 }
592
593 static long compute_num_played_score(struct afs_info *afsi)
594 {
595         return -normalized_value(afsi->num_played, statistics.num,
596                 statistics.num_played_sum, statistics.num_played_qd);
597 }
598
599 static long compute_last_played_score(struct afs_info *afsi)
600 {
601         return -normalized_value(afsi->last_played, statistics.num,
602                 statistics.last_played_sum, statistics.last_played_qd);
603 }
604
605 static long compute_dynamic_score(const struct osl_row *aft_row)
606 {
607         struct afs_info afsi;
608         int64_t score, nscore = 0, lscore = 0;
609         int ret;
610
611         ret = get_afsi_of_row(aft_row, &afsi);
612         if (ret < 0)
613                 return -100;
614         nscore = compute_num_played_score(&afsi);
615         lscore = compute_last_played_score(&afsi);
616         score = nscore + lscore;
617         return score;
618 }
619
620 static int add_afs_statistics(const struct osl_row *row)
621 {
622         uint64_t n, x, s;
623         struct afs_info afsi;
624         int ret;
625
626         ret = get_afsi_of_row(row, &afsi);
627         if (ret < 0)
628                 return ret;
629         n = statistics.num;
630         x = afsi.last_played;
631         s = statistics.last_played_sum;
632         if (n > 0)
633                 statistics.last_played_qd += (x - s / n) * (x - s / n) * n / (n + 1);
634         statistics.last_played_sum += x;
635
636         x = afsi.num_played;
637         s = statistics.num_played_sum;
638         if (n > 0)
639                 statistics.num_played_qd += (x - s / n) * (x - s / n) * n / (n + 1);
640         statistics.num_played_sum += x;
641         statistics.num++;
642         return 1;
643 }
644
645 static int del_afs_statistics(const struct osl_row *row)
646 {
647         uint64_t n, s, q, a, new_s;
648         struct afs_info afsi;
649         int ret;
650         ret = get_afsi_of_row(row, &afsi);
651         if (ret < 0)
652                 return ret;
653         n = statistics.num;
654         assert(n);
655         if (n == 1) {
656                 memset(&statistics, 0, sizeof(statistics));
657                 return 1;
658         }
659
660         s = statistics.last_played_sum;
661         q = statistics.last_played_qd;
662         a = afsi.last_played;
663         new_s = s - a;
664         statistics.last_played_sum = new_s;
665         statistics.last_played_qd = q + s * s / n - a * a
666                 - new_s * new_s / (n - 1);
667
668         s = statistics.num_played_sum;
669         q = statistics.num_played_qd;
670         a = afsi.num_played;
671         new_s = s - a;
672         statistics.num_played_sum = new_s;
673         statistics.num_played_qd = q + s * s / n - a * a
674                 - new_s * new_s / (n - 1);
675
676         statistics.num--;
677         return 1;
678 }
679
680 /**
681  * Structure used during mood_open().
682  *
683  * At mood open time, we look at each file in the audio file table in order to
684  * determine whether it is admissible. If a file happens to be admissible, its
685  * mood score is computed by calling each relevant mood_score_function. Next,
686  * we update the afs_statistics and add a struct admissible_file_info to a
687  * temporary array.
688  *
689  * If all files have been processed that way, the final score of each
690  * admissible file is computed by adding the dynamic score (which depends on
691  * the afs_statistics) to the mood score.  Finally, all audio files in the
692  * array are added to the score table and the admissible array is freed.
693  *
694  * \sa mood_method, admissible_array.
695  */
696 struct admissible_file_info
697 {
698         /** The admissible audio file. */
699         struct osl_row *aft_row;
700         /** Its score. */
701         long score;
702 };
703
704 /** The temporary array of admissible files. */
705 struct admissible_array {
706         /** Files are admissible wrt. this mood. */
707         struct mood *m;
708         /** The size of the array */
709         unsigned size;
710         /** Pointer to the array of admissible files. */
711         struct admissible_file_info *array;
712 };
713
714 /**
715  * Add an entry to the array of admissible files.
716  *
717  * \param aft_row The audio file to be added.
718  * \param private_data Pointer to a struct admissible_file_info.
719  *
720  * \return Negative on errors, positive on success.
721  */
722 static int add_if_admissible(struct osl_row *aft_row, void *data)
723 {
724         struct admissible_array *aa = data;
725         int ret;
726         long score = 0;
727
728         ret = compute_mood_score(aft_row, aa->m, &score);
729         if (ret < 0)
730                 return (ret == -E_NOT_ADMISSIBLE)? 1 : ret;
731         if (statistics.num >= aa->size) {
732                 aa->size *= 2;
733                 aa->size += 100;
734                 aa->array = para_realloc(aa->array,
735                         aa->size * sizeof(struct admissible_file_info));
736         }
737         aa->array[statistics.num].aft_row = aft_row;
738         aa->array[statistics.num].score = score;
739         ret = add_afs_statistics(aft_row);
740         if (ret < 0)
741                 return ret;
742         return 1;
743 }
744
745 /**
746  * Compute the new quadratic deviation in case one element changes.
747  *
748  * \param n Number of elements.
749  * \param old_qd The quadratic deviation before the change.
750  * \param old_val The value that was replaced.
751  * \param new_val The replacement value.
752  * \param old_sum The sum of all elements before the update.
753  *
754  * \return The new quadratic deviation resulting from replacing old_val
755  * by new_val.
756  *
757  * Given n real numbers a_1, ..., a_n, their sum S = a_1 + ... + a_n,
758  * their quadratic deviation
759  *
760  * q = (a_1 - S/n)^2 + ... + (a_n - S/n)^2,
761  *
762  * and a real number b, the quadratic deviation q' of a_1,...a_{n-1}, b (ie.
763  * the last number a_n was replaced by b) may be computed in O(1) time in terms
764  * of n, q, a_n, b, and S as
765  *
766  *      q' = q + d * s - (2 * S + d) * d / n,
767  *
768  * where d = b - a_n, and s = b + a_n.
769  *
770  * Example: n = 3, a_1 = 3, a_2 = 5, a_3 = 7, b = 10. Then S = 15, q = 8, d = 3,
771  * s = 17, so
772  *
773  *      q + d * s - (2 * S + d) * d / n = 8 + 51 - 33 = 26,
774  *
775  * which equals q' = (3 - 6)^2 + (5 - 6)^2 + (10 - 6)^2.
776  *
777  */
778 _static_inline_ int64_t update_quadratic_deviation(int64_t n, int64_t old_qd,
779                 int64_t old_val, int64_t new_val, int64_t old_sum)
780 {
781         int64_t delta = new_val - old_val;
782         int64_t sigma = new_val + old_val;
783         return old_qd + delta * sigma - (2 * old_sum + delta) * delta / n;
784 }
785
786 static int update_afs_statistics(struct afs_info *old_afsi, struct afs_info *new_afsi)
787 {
788         unsigned n;
789         int ret = get_num_admissible_files(&n);
790
791         if (ret < 0)
792                 return ret;
793         assert(n);
794
795         statistics.last_played_qd = update_quadratic_deviation(n,
796                 statistics.last_played_qd, old_afsi->last_played,
797                 new_afsi->last_played, statistics.last_played_sum);
798         statistics.last_played_sum += new_afsi->last_played - old_afsi->last_played;
799
800         statistics.num_played_qd = update_quadratic_deviation(n,
801                 statistics.num_played_qd, old_afsi->num_played,
802                 new_afsi->num_played, statistics.num_played_sum);
803         statistics.num_played_sum += new_afsi->num_played - old_afsi->num_played;
804         return 1;
805 }
806
807 static int add_to_score_table(const struct osl_row *aft_row, long mood_score)
808 {
809         long score = (compute_dynamic_score(aft_row) + mood_score) / 3;
810         return score_add(aft_row, score);
811 }
812
813 static int delete_from_statistics_and_score_table(const struct osl_row *aft_row)
814 {
815         int ret = del_afs_statistics(aft_row);
816         if (ret < 0)
817                 return ret;
818         return score_delete(aft_row);
819 }
820
821 /**
822  * Delete one entry from the statistics and from the score table.
823  *
824  * \param aft_row The audio file which is no longer admissible.
825  *
826  * \return Positive on success, negative on errors.
827  *
828  * \sa score_delete().
829  */
830 static int mood_delete_audio_file(const struct osl_row *aft_row)
831 {
832         int ret;
833
834         ret = row_belongs_to_score_table(aft_row, NULL);
835         if (ret < 0)
836                 return ret;
837         if (!ret) /* not admissible, nothing to do */
838                 return 1;
839         return delete_from_statistics_and_score_table(aft_row);
840 }
841
842 /**
843  * Compute the new score of an audio file wrt. the current mood.
844  *
845  * \param aft_row Determines the audio file.
846  * \param old_afsi The audio file selector info before updating.
847  *
848  * The \a old_afsi argument may be \p NULL which indicates that no changes to
849  * the audio file info were made.
850  *
851  * \return Positive on success, negative on errors.
852  */
853 static int mood_update_audio_file(const struct osl_row *aft_row,
854                 struct afs_info *old_afsi)
855 {
856         long score, percent;
857         int ret, is_admissible, was_admissible = 0;
858         struct afs_info afsi;
859         unsigned rank;
860
861         if (!current_mood)
862                 return 1; /* nothing to do */
863         ret = row_belongs_to_score_table(aft_row, &rank);
864         if (ret < 0)
865                 return ret;
866         was_admissible = ret;
867         ret = compute_mood_score(aft_row, current_mood, &score);
868         is_admissible = (ret > 0);
869         if (!was_admissible && !is_admissible)
870                 return 1;
871         if (was_admissible && !is_admissible)
872                 return delete_from_statistics_and_score_table(aft_row);
873         if (!was_admissible && is_admissible) {
874                 ret = add_afs_statistics(aft_row);
875                 if (ret < 0)
876                         return ret;
877                 return add_to_score_table(aft_row, score);
878         }
879         /* update score */
880         ret = get_afsi_of_row(aft_row, &afsi);
881         if (ret < 0)
882                 return ret;
883         if (old_afsi) {
884                 ret = update_afs_statistics(old_afsi, &afsi);
885                 if (ret < 0)
886                         return ret;
887         }
888         score += compute_num_played_score(&afsi);
889         score += compute_last_played_score(&afsi);
890         score /= 3;
891         PARA_DEBUG_LOG("score: %li\n", score);
892         percent = (score + 100) / 3;
893         if (percent > 100)
894                 percent = 100;
895         else if (percent < 0)
896                 percent = 0;
897         PARA_DEBUG_LOG("moving from rank %u to %lu%%\n", rank, percent);
898         return score_update(aft_row, percent);
899 }
900
901 static void log_statistics(void)
902 {
903         unsigned n = statistics.num;
904
905         if (!n) {
906                 PARA_NOTICE_LOG("no admissible files\n");
907                 return;
908         }
909         PARA_INFO_LOG("last_played mean: %lli, last_played sigma: %llu\n",
910                 (long long int)(statistics.last_played_sum / n),
911                 (long long unsigned)int_sqrt(statistics.last_played_qd / n));
912         PARA_INFO_LOG("num_played mean: %lli, num_played sigma: %llu\n",
913                 (long long int)statistics.num_played_sum / n,
914                 (long long unsigned)int_sqrt(statistics.num_played_qd / n));
915 }
916
917 /**
918  * Close the current mood.
919  *
920  * Free all resources of the current mood which were allocated during
921  * mood_open().
922  */
923 void close_current_mood(void)
924 {
925         destroy_mood(current_mood);
926         current_mood = NULL;
927         memset(&statistics, 0, sizeof(statistics));
928 }
929
930
931 /**
932  * Change the current mood.
933  *
934  * \param mood_name The name of the mood to open.
935  *
936  * If \a mood_name is \a NULL, load the dummy mood that accepts every audio file
937  * and uses a scoring method based only on the \a last_played information.
938  *
939  * If there is already an open mood, it will be closed first.
940  *
941  * \return Positive on success, negative on errors. Loading the dummy mood
942  * always succeeds.
943  *
944  * \sa struct admissible_file_info, struct admissible_array, struct
945  * afs_info::last_played, mood_close().
946  */
947 int change_current_mood(char *mood_name)
948 {
949         int i, ret;
950         struct admissible_array aa = {
951                 .size = 0,
952                 .array = NULL
953         };
954
955         if (mood_name) {
956                 struct mood *m;
957                 struct osl_row *row;
958                 struct osl_object obj = {
959                         .data = mood_name,
960                         .size = strlen(mood_name) + 1
961                 };
962                 ret = osl_get_row(moods_table, BLOBCOL_NAME, &obj, &row);
963                 if (ret < 0) {
964                         PARA_NOTICE_LOG("no such mood: %s\n", mood_name);
965                         return ret;
966                 }
967                 ret = load_mood(row, &m);
968                 if (ret < 0)
969                         return ret;
970                 close_current_mood();
971                 current_mood = m;
972         } else {
973                 close_current_mood();
974                 current_mood = alloc_new_mood("dummy");
975         }
976         aa.m = current_mood;
977         PARA_NOTICE_LOG("computing statistics of admissible files\n");
978         ret = audio_file_loop(&aa, add_if_admissible);
979         if (ret < 0)
980                 return ret;
981         log_statistics();
982         PARA_INFO_LOG("%d admissible files \n", statistics.num);
983         for (i = 0; i < statistics.num; i++) {
984                 struct admissible_file_info *a = aa.array + i;
985                 ret = add_to_score_table(a->aft_row, a->score);
986                 if (ret < 0)
987                         goto out;
988         }
989         PARA_NOTICE_LOG("loaded mood %s\n", current_mood->name);
990         ret = statistics.num;
991 out:
992         free(aa.array);
993         return ret;
994 }
995 /**
996  * Close and re-open the current mood.
997  *
998  * This function is used if changes to the audio file table or the
999  * attribute table were made that render the current list of admissible
1000  * files useless. For example, if an attribute is removed from the
1001  * attribute table, this function is called.
1002  *
1003  * \return Positive on success, negative on errors. If no mood is currently
1004  * open, the function returns success.
1005  *
1006  * \sa mood_open(), mood_close().
1007  */
1008 int reload_current_mood(void)
1009 {
1010         int ret;
1011         char *mood_name = NULL;
1012
1013         PARA_NOTICE_LOG("reloading current mood\n");
1014         if (!current_mood)
1015                 return 1;
1016         if (current_mood->name)
1017                 mood_name = para_strdup(current_mood->name);
1018         close_current_mood();
1019         ret = change_current_mood(mood_name);
1020         free(mood_name);
1021         return ret;
1022 }
1023
1024 int moods_event_handler(enum afs_events event, __a_unused struct para_buffer *pb,
1025                 void *data)
1026 {
1027         switch(event) {
1028         /*
1029          * The three blob events might change the set of admissible files,
1030          * so we must reload the score list.
1031          */
1032         case BLOB_RENAME:
1033         case BLOB_REMOVE:
1034         case BLOB_ADD:
1035                 if (data == moods_table || data == playlists_table)
1036                         return 1; /* no reload necessary for these */
1037                 return reload_current_mood();
1038         /* these also require reload of the score table */
1039         case ATTRIBUTE_ADD:
1040         case ATTRIBUTE_REMOVE:
1041         case ATTRIBUTE_RENAME:
1042                 return reload_current_mood();
1043         /* changes to the aft only require to re-examine the audio file */
1044         case AFSI_CHANGE: {
1045                 struct afsi_change_event_data *aced = data;
1046                 return mood_update_audio_file(aced->aft_row, aced->old_afsi);
1047                 }
1048         case AFHI_CHANGE:
1049         case AUDIO_FILE_RENAME:
1050         case AUDIO_FILE_ADD:
1051                 return mood_update_audio_file(data, NULL);
1052         case AUDIO_FILE_REMOVE:
1053                 return mood_delete_audio_file(data);
1054         default:
1055                 return 1;
1056         }
1057         return 1;
1058 }
1059