]> git.sesse.net Git - plocate/blob - plocate.cpp
Unbreak non-trigram matches after we changed to asynchronous access().
[plocate] / plocate.cpp
1 #include "db.h"
2 #include "dprintf.h"
3 #include "io_uring_engine.h"
4 #include "parse_trigrams.h"
5 #include "turbopfor.h"
6 #include "unique_sort.h"
7
8 #include <algorithm>
9 #include <assert.h>
10 #include <chrono>
11 #include <fcntl.h>
12 #include <fnmatch.h>
13 #include <functional>
14 #include <getopt.h>
15 #include <inttypes.h>
16 #include <iosfwd>
17 #include <iterator>
18 #include <limits>
19 #include <map>
20 #include <memory>
21 #include <queue>
22 #include <regex.h>
23 #include <stdint.h>
24 #include <stdio.h>
25 #include <stdlib.h>
26 #include <string.h>
27 #include <string>
28 #include <string_view>
29 #include <unistd.h>
30 #include <unordered_map>
31 #include <unordered_set>
32 #include <utility>
33 #include <vector>
34 #include <zstd.h>
35
36 using namespace std;
37 using namespace std::chrono;
38
39 #define DEFAULT_DBPATH "/var/lib/mlocate/plocate.db"
40
41 const char *dbpath = DEFAULT_DBPATH;
42 bool ignore_case = false;
43 bool only_count = false;
44 bool print_nul = false;
45 bool use_debug = false;
46 bool patterns_are_regex = false;
47 bool use_extended_regex = false;
48 int64_t limit_matches = numeric_limits<int64_t>::max();
49 int64_t limit_left = numeric_limits<int64_t>::max();
50
51 steady_clock::time_point start;
52
53 void apply_limit()
54 {
55         if (--limit_left > 0) {
56                 return;
57         }
58         dprintf("Done in %.1f ms, found %" PRId64 " matches.\n",
59                 1e3 * duration<float>(steady_clock::now() - start).count(), limit_matches);
60         if (only_count) {
61                 printf("%ld\n", limit_matches);
62         }
63         exit(0);
64 }
65
66 class Serializer {
67 public:
68         ~Serializer() { assert(limit_left <= 0 || pending.empty()); }
69         void print(uint64_t seq, uint64_t skip, const string msg);
70
71 private:
72         uint64_t next_seq = 0;
73         struct Element {
74                 uint64_t seq, skip;
75                 string msg;
76
77                 bool operator<(const Element &other) const
78                 {
79                         return seq > other.seq;
80                 }
81         };
82         priority_queue<Element> pending;
83 };
84
85 void Serializer::print(uint64_t seq, uint64_t skip, const string msg)
86 {
87         if (only_count) {
88                 if (!msg.empty()) {
89                         apply_limit();
90                 }
91                 return;
92         }
93
94         if (next_seq != seq) {
95                 pending.push(Element{ seq, skip, move(msg) });
96                 return;
97         }
98
99         if (!msg.empty()) {
100                 if (print_nul) {
101                         printf("%s%c", msg.c_str(), 0);
102                 } else {
103                         printf("%s\n", msg.c_str());
104                 }
105                 apply_limit();
106         }
107         next_seq += skip;
108
109         // See if any delayed prints can now be dealt with.
110         while (!pending.empty() && pending.top().seq == next_seq) {
111                 if (!pending.top().msg.empty()) {
112                         if (print_nul) {
113                                 printf("%s%c", pending.top().msg.c_str(), 0);
114                         } else {
115                                 printf("%s\n", pending.top().msg.c_str());
116                         }
117                         apply_limit();
118                 }
119                 next_seq += pending.top().skip;
120                 pending.pop();
121         }
122 }
123
124 struct Needle {
125         enum { STRSTR,
126                REGEX,
127                GLOB } type;
128         string str;  // Filled in no matter what.
129         regex_t re;  // For REGEX.
130 };
131
132 bool matches(const Needle &needle, const char *haystack)
133 {
134         if (needle.type == Needle::STRSTR) {
135                 return strstr(haystack, needle.str.c_str()) != nullptr;
136         } else if (needle.type == Needle::GLOB) {
137                 int flags = ignore_case ? FNM_CASEFOLD : 0;
138                 return fnmatch(needle.str.c_str(), haystack, flags) == 0;
139         } else {
140                 assert(needle.type == Needle::REGEX);
141                 return regexec(&needle.re, haystack, /*nmatch=*/0, /*pmatch=*/nullptr, /*flags=*/0) == 0;
142         }
143 }
144
145 class AccessRXCache {
146 public:
147         AccessRXCache(IOUringEngine *engine)
148                 : engine(engine) {}
149         void check_access(const char *filename, bool allow_async, function<void(bool)> cb);
150
151 private:
152         unordered_map<string, bool> cache;
153         struct PendingStat {
154                 string filename;
155                 function<void(bool)> cb;
156         };
157         map<string, vector<PendingStat>> pending_stats;
158         IOUringEngine *engine;
159 };
160
161 void AccessRXCache::check_access(const char *filename, bool allow_async, function<void(bool)> cb)
162 {
163         if (engine == nullptr || !engine->get_supports_stat()) {
164                 allow_async = false;
165         }
166
167         for (const char *end = strchr(filename + 1, '/'); end != nullptr; end = strchr(end + 1, '/')) {
168                 string parent_path(filename, end - filename);  // string_view from C++20.
169                 auto cache_it = cache.find(parent_path);
170                 if (cache_it != cache.end()) {
171                         // Found in the cache.
172                         if (!cache_it->second) {
173                                 cb(false);
174                                 return;
175                         }
176                         continue;
177                 }
178
179                 if (!allow_async) {
180                         bool ok = access(parent_path.c_str(), R_OK | X_OK) == 0;
181                         cache.emplace(parent_path, ok);
182                         if (!ok) {
183                                 cb(false);
184                                 return;
185                         }
186                         continue;
187                 }
188
189                 // We want to call access(), but it could block on I/O. io_uring doesn't support
190                 // access(), but we can do a dummy asynchonous statx() to populate the kernel's cache,
191                 // which nearly always makes the next access() instantaneous.
192
193                 // See if there's already a pending stat that matches this,
194                 // or is a subdirectory.
195                 auto it = pending_stats.lower_bound(parent_path);
196                 if (it != pending_stats.end() && it->first.size() >= parent_path.size() &&
197                     it->first.compare(0, parent_path.size(), parent_path) == 0) {
198                         it->second.emplace_back(PendingStat{ filename, move(cb) });
199                 } else {
200                         it = pending_stats.emplace(filename, vector<PendingStat>{}).first;
201                         engine->submit_stat(filename, [this, it, filename{ strdup(filename) }, cb{ move(cb) }] {
202                                 // The stat returned, so now do the actual access() calls.
203                                 // All of them should be in cache, so don't fire off new statx()
204                                 // calls during that check.
205                                 check_access(filename, /*allow_async=*/false, move(cb));
206                                 free(filename);
207
208                                 // Call all others that waited for the same stat() to finish.
209                                 // They may fire off new stat() calls if needed.
210                                 vector<PendingStat> pending = move(it->second);
211                                 pending_stats.erase(it);
212                                 for (PendingStat &ps : pending) {
213                                         check_access(ps.filename.c_str(), /*allow_async=*/true, move(ps.cb));
214                                 }
215                         });
216                 }
217                 return;  // The rest will happen in async context.
218         }
219
220         // Passed all checks.
221         cb(true);
222 }
223
224 class Corpus {
225 public:
226         Corpus(int fd, IOUringEngine *engine);
227         ~Corpus();
228         void find_trigram(uint32_t trgm, function<void(const Trigram *trgmptr, size_t len)> cb);
229         void get_compressed_filename_block(uint32_t docid, function<void(string_view)> cb) const;
230         size_t get_num_filename_blocks() const;
231         off_t offset_for_block(uint32_t docid) const
232         {
233                 return hdr.filename_index_offset_bytes + docid * sizeof(uint64_t);
234         }
235
236 public:
237         const int fd;
238         IOUringEngine *const engine;
239
240         Header hdr;
241 };
242
243 Corpus::Corpus(int fd, IOUringEngine *engine)
244         : fd(fd), engine(engine)
245 {
246         // Enable to test cold-cache behavior (except for access()).
247         if (false) {
248                 off_t len = lseek(fd, 0, SEEK_END);
249                 if (len == -1) {
250                         perror("lseek");
251                         exit(1);
252                 }
253                 posix_fadvise(fd, 0, len, POSIX_FADV_DONTNEED);
254         }
255
256         complete_pread(fd, &hdr, sizeof(hdr), /*offset=*/0);
257         if (memcmp(hdr.magic, "\0plocate", 8) != 0) {
258                 fprintf(stderr, "plocate.db is corrupt or an old version; please rebuild it.\n");
259                 exit(1);
260         }
261         if (hdr.version != 0) {
262                 fprintf(stderr, "plocate.db has version %u, expected 0; please rebuild it.\n", hdr.version);
263                 exit(1);
264         }
265 }
266
267 Corpus::~Corpus()
268 {
269         close(fd);
270 }
271
272 void Corpus::find_trigram(uint32_t trgm, function<void(const Trigram *trgmptr, size_t len)> cb)
273 {
274         uint32_t bucket = hash_trigram(trgm, hdr.hashtable_size);
275         engine->submit_read(fd, sizeof(Trigram) * (hdr.extra_ht_slots + 2), hdr.hash_table_offset_bytes + sizeof(Trigram) * bucket, [this, trgm, cb{ move(cb) }](string_view s) {
276                 const Trigram *trgmptr = reinterpret_cast<const Trigram *>(s.data());
277                 for (unsigned i = 0; i < hdr.extra_ht_slots + 1; ++i) {
278                         if (trgmptr[i].trgm == trgm) {
279                                 cb(trgmptr + i, trgmptr[i + 1].offset - trgmptr[i].offset);
280                                 return;
281                         }
282                 }
283
284                 // Not found.
285                 cb(nullptr, 0);
286         });
287 }
288
289 void Corpus::get_compressed_filename_block(uint32_t docid, function<void(string_view)> cb) const
290 {
291         // Read the file offset from this docid and the next one.
292         // This is always allowed, since we have a sentinel block at the end.
293         engine->submit_read(fd, sizeof(uint64_t) * 2, offset_for_block(docid), [this, cb{ move(cb) }](string_view s) {
294                 const uint64_t *ptr = reinterpret_cast<const uint64_t *>(s.data());
295                 off_t offset = ptr[0];
296                 size_t len = ptr[1] - ptr[0];
297                 engine->submit_read(fd, len, offset, cb);
298         });
299 }
300
301 size_t Corpus::get_num_filename_blocks() const
302 {
303         return hdr.num_docids;
304 }
305
306 void scan_file_block(const vector<Needle> &needles, string_view compressed,
307                      AccessRXCache *access_rx_cache, uint64_t seq, Serializer *serializer,
308                      size_t *matched)
309 {
310         unsigned long long uncompressed_len = ZSTD_getFrameContentSize(compressed.data(), compressed.size());
311         if (uncompressed_len == ZSTD_CONTENTSIZE_UNKNOWN || uncompressed_len == ZSTD_CONTENTSIZE_ERROR) {
312                 fprintf(stderr, "ZSTD_getFrameContentSize() failed\n");
313                 exit(1);
314         }
315
316         string block;
317         block.resize(uncompressed_len + 1);
318
319         size_t err = ZSTD_decompress(&block[0], block.size(), compressed.data(),
320                                      compressed.size());
321         if (ZSTD_isError(err)) {
322                 fprintf(stderr, "ZSTD_decompress(): %s\n", ZSTD_getErrorName(err));
323                 exit(1);
324         }
325         block[block.size() - 1] = '\0';
326
327         auto test_candidate = [&](const char *filename, uint64_t local_seq, uint64_t next_seq) {
328                 access_rx_cache->check_access(filename, /*allow_async=*/true, [matched, serializer, local_seq, next_seq, filename{ strdup(filename) }](bool ok) {
329                         if (ok) {
330                                 ++*matched;
331                                 serializer->print(local_seq, next_seq - local_seq, filename);
332                         } else {
333                                 serializer->print(local_seq, next_seq - local_seq, "");
334                         }
335                         free(filename);
336                 });
337         };
338
339         // We need to know the next sequence number before inserting into Serializer,
340         // so always buffer one candidate.
341         const char *pending_candidate = nullptr;
342
343         uint64_t local_seq = seq << 32;
344         for (const char *filename = block.data();
345              filename != block.data() + block.size();
346              filename += strlen(filename) + 1) {
347                 bool found = true;
348                 for (const Needle &needle : needles) {
349                         if (!matches(needle, filename)) {
350                                 found = false;
351                                 break;
352                         }
353                 }
354                 if (found) {
355                         if (pending_candidate != nullptr) {
356                                 test_candidate(pending_candidate, local_seq, local_seq + 1);
357                                 ++local_seq;
358                         }
359                         pending_candidate = filename;
360                 }
361         }
362         if (pending_candidate == nullptr) {
363                 serializer->print(seq << 32, 1ULL << 32, "");
364         } else {
365                 test_candidate(pending_candidate, local_seq, (seq + 1) << 32);
366         }
367 }
368
369 size_t scan_docids(const vector<Needle> &needles, const vector<uint32_t> &docids, const Corpus &corpus, IOUringEngine *engine)
370 {
371         Serializer docids_in_order;
372         AccessRXCache access_rx_cache(engine);
373         uint64_t matched = 0;
374         for (size_t i = 0; i < docids.size(); ++i) {
375                 uint32_t docid = docids[i];
376                 corpus.get_compressed_filename_block(docid, [i, &matched, &needles, &access_rx_cache, &docids_in_order](string_view compressed) {
377                         scan_file_block(needles, compressed, &access_rx_cache, i, &docids_in_order, &matched);
378                 });
379         }
380         engine->finish();
381         return matched;
382 }
383
384 // We do this sequentially, as it's faster than scattering
385 // a lot of I/O through io_uring and hoping the kernel will
386 // coalesce it plus readahead for us.
387 uint64_t scan_all_docids(const vector<Needle> &needles, int fd, const Corpus &corpus)
388 {
389         AccessRXCache access_rx_cache(nullptr);
390         Serializer serializer;  // Mostly dummy; handles only the limit.
391         uint32_t num_blocks = corpus.get_num_filename_blocks();
392         unique_ptr<uint64_t[]> offsets(new uint64_t[num_blocks + 1]);
393         complete_pread(fd, offsets.get(), (num_blocks + 1) * sizeof(uint64_t), corpus.offset_for_block(0));
394         string compressed;
395         uint64_t matched = 0;
396         for (uint32_t io_docid = 0; io_docid < num_blocks; io_docid += 32) {
397                 uint32_t last_docid = std::min(io_docid + 32, num_blocks);
398                 size_t io_len = offsets[last_docid] - offsets[io_docid];
399                 if (compressed.size() < io_len) {
400                         compressed.resize(io_len);
401                 }
402                 complete_pread(fd, &compressed[0], io_len, offsets[io_docid]);
403
404                 for (uint32_t docid = io_docid; docid < last_docid; ++docid) {
405                         size_t relative_offset = offsets[docid] - offsets[io_docid];
406                         size_t len = offsets[docid + 1] - offsets[docid];
407                         scan_file_block(needles, { &compressed[relative_offset], len }, &access_rx_cache, docid, &serializer, &matched);
408                 }
409         }
410         return matched;
411 }
412
413 // Takes the given posting list, unions it into the parts of the trigram disjunction
414 // already read; if the list is complete, intersects with “cur_candidates”.
415 //
416 // Returns true if the search should be aborted (we are done).
417 bool new_posting_list_read(TrigramDisjunction *td, vector<uint32_t> decoded, vector<uint32_t> *cur_candidates, vector<uint32_t> *tmp)
418 {
419         if (td->docids.empty()) {
420                 td->docids = move(decoded);
421         } else {
422                 tmp->clear();
423                 set_union(decoded.begin(), decoded.end(), td->docids.begin(), td->docids.end(), back_inserter(*tmp));
424                 swap(*tmp, td->docids);
425         }
426         if (--td->remaining_trigrams_to_read > 0) {
427                 // Need to wait for more.
428                 if (ignore_case) {
429                         dprintf("  ... %u reads left in OR group %u (%zu docids in list)\n",
430                                 td->remaining_trigrams_to_read, td->index, td->docids.size());
431                 }
432                 return false;
433         }
434         if (cur_candidates->empty()) {
435                 if (ignore_case) {
436                         dprintf("  ... all reads done for OR group %u (%zu docids)\n",
437                                 td->index, td->docids.size());
438                 }
439                 *cur_candidates = move(td->docids);
440         } else {
441                 tmp->clear();
442                 set_intersection(cur_candidates->begin(), cur_candidates->end(),
443                                  td->docids.begin(), td->docids.end(),
444                                  back_inserter(*tmp));
445                 swap(*cur_candidates, *tmp);
446                 if (ignore_case) {
447                         if (cur_candidates->empty()) {
448                                 dprintf("  ... all reads done for OR group %u (%zu docids), intersected (none left, search is done)\n",
449                                         td->index, td->docids.size());
450                                 return true;
451                         } else {
452                                 dprintf("  ... all reads done for OR group %u (%zu docids), intersected (%zu left)\n",
453                                         td->index, td->docids.size(), cur_candidates->size());
454                         }
455                 }
456         }
457         return false;
458 }
459
460 void do_search_file(const vector<Needle> &needles, const char *filename)
461 {
462         int fd = open(filename, O_RDONLY);
463         if (fd == -1) {
464                 perror(filename);
465                 exit(1);
466         }
467
468         // Drop privileges.
469         if (setgid(getgid()) != 0) {
470                 perror("setgid");
471                 exit(EXIT_FAILURE);
472         }
473
474         start = steady_clock::now();
475         if (access("/", R_OK | X_OK)) {
476                 // We can't find anything, no need to bother...
477                 return;
478         }
479
480         IOUringEngine engine(/*slop_bytes=*/16);  // 16 slop bytes as described in turbopfor.h.
481         Corpus corpus(fd, &engine);
482         dprintf("Corpus init done after %.1f ms.\n", 1e3 * duration<float>(steady_clock::now() - start).count());
483
484         vector<TrigramDisjunction> trigram_groups;
485         if (patterns_are_regex) {
486                 // We could parse the regex to find trigrams that have to be there
487                 // (there are actually known algorithms to deal with disjunctions
488                 // and such, too), but for now, we just go brute force.
489                 // Using locate with regexes is pretty niche.
490         } else {
491                 for (const Needle &needle : needles) {
492                         parse_trigrams(needle.str, ignore_case, &trigram_groups);
493                 }
494         }
495
496         unique_sort(
497                 &trigram_groups,
498                 [](const TrigramDisjunction &a, const TrigramDisjunction &b) { return a.trigram_alternatives < b.trigram_alternatives; },
499                 [](const TrigramDisjunction &a, const TrigramDisjunction &b) { return a.trigram_alternatives == b.trigram_alternatives; });
500
501         // Give them names for debugging.
502         unsigned td_index = 0;
503         for (TrigramDisjunction &td : trigram_groups) {
504                 td.index = td_index++;
505         }
506
507         // Collect which trigrams we need to look up in the hash table.
508         unordered_map<uint32_t, vector<TrigramDisjunction *>> trigrams_to_lookup;
509         for (TrigramDisjunction &td : trigram_groups) {
510                 for (uint32_t trgm : td.trigram_alternatives) {
511                         trigrams_to_lookup[trgm].push_back(&td);
512                 }
513         }
514         if (trigrams_to_lookup.empty()) {
515                 // Too short for trigram matching. Apply brute force.
516                 // (We could have searched through all trigrams that matched
517                 // the pattern and done a union of them, but that's a lot of
518                 // work for fairly unclear gain.)
519                 uint64_t matched = scan_all_docids(needles, fd, corpus);
520                 if (only_count) {
521                         printf("%" PRId64 "\n", matched);
522                 }
523                 return;
524         }
525
526         // Look them all up on disk.
527         for (auto &[trgm, trigram_groups] : trigrams_to_lookup) {
528                 corpus.find_trigram(trgm, [trgm{ trgm }, trigram_groups{ &trigram_groups }](const Trigram *trgmptr, size_t len) {
529                         if (trgmptr == nullptr) {
530                                 dprintf("trigram %s isn't found\n", print_trigram(trgm).c_str());
531                                 for (TrigramDisjunction *td : *trigram_groups) {
532                                         --td->remaining_trigrams_to_read;
533                                         if (td->remaining_trigrams_to_read == 0 && td->read_trigrams.empty()) {
534                                                 dprintf("zero matches in %s, so we are done\n", print_td(*td).c_str());
535                                                 if (only_count) {
536                                                         printf("0\n");
537                                                 }
538                                                 exit(0);
539                                         }
540                                 }
541                                 return;
542                         }
543                         for (TrigramDisjunction *td : *trigram_groups) {
544                                 --td->remaining_trigrams_to_read;
545                                 td->max_num_docids += trgmptr->num_docids;
546                                 td->read_trigrams.emplace_back(*trgmptr, len);
547                         }
548                 });
549         }
550         engine.finish();
551         dprintf("Hashtable lookups done after %.1f ms.\n", 1e3 * duration<float>(steady_clock::now() - start).count());
552
553         for (TrigramDisjunction &td : trigram_groups) {
554                 // Reset for reads.
555                 td.remaining_trigrams_to_read = td.read_trigrams.size();
556
557                 if (ignore_case) {  // If case-sensitive, they'll all be pretty obvious single-entry groups.
558                         dprintf("OR group %u (max_num_docids=%u): %s\n", td.index, td.max_num_docids, print_td(td).c_str());
559                 }
560         }
561
562         // TODO: For case-insensitive (ie. more than one alternative in each),
563         // prioritize the ones with fewer seeks?
564         sort(trigram_groups.begin(), trigram_groups.end(),
565              [&](const TrigramDisjunction &a, const TrigramDisjunction &b) {
566                      return a.max_num_docids < b.max_num_docids;
567              });
568
569         unordered_map<uint32_t, vector<TrigramDisjunction *>> uses_trigram;
570         for (TrigramDisjunction &td : trigram_groups) {
571                 for (uint32_t trgm : td.trigram_alternatives) {
572                         uses_trigram[trgm].push_back(&td);
573                 }
574         }
575
576         unordered_set<uint32_t> trigrams_submitted_read;
577         vector<uint32_t> cur_candidates, tmp, decoded;
578         bool done = false;
579         for (TrigramDisjunction &td : trigram_groups) {
580                 if (!cur_candidates.empty() && td.max_num_docids > cur_candidates.size() * 100) {
581                         dprintf("%s has up to %u entries, ignoring the rest (will "
582                                 "weed out false positives later)\n",
583                                 print_td(td).c_str(), td.max_num_docids);
584                         break;
585                 }
586
587                 for (auto &[trgmptr, len] : td.read_trigrams) {
588                         if (trigrams_submitted_read.count(trgmptr.trgm) != 0) {
589                                 continue;
590                         }
591                         trigrams_submitted_read.insert(trgmptr.trgm);
592                         // Only stay a certain amount ahead, so that we don't spend I/O
593                         // on reading the latter, large posting lists. We are unlikely
594                         // to need them anyway, even if they should come in first.
595                         if (engine.get_waiting_reads() >= 5) {
596                                 engine.finish();
597                                 if (done)
598                                         break;
599                         }
600                         engine.submit_read(fd, len, trgmptr.offset, [trgmptr{ trgmptr }, len{ len }, &done, &cur_candidates, &tmp, &decoded, &uses_trigram](string_view s) {
601                                 if (done)
602                                         return;
603
604                                 uint32_t trgm __attribute__((unused)) = trgmptr.trgm;
605                                 const unsigned char *pldata = reinterpret_cast<const unsigned char *>(s.data());
606                                 size_t num = trgmptr.num_docids;
607                                 decoded.resize(num);
608                                 decode_pfor_delta1_128(pldata, num, /*interleaved=*/true, &decoded[0]);
609
610                                 assert(uses_trigram.count(trgm) != 0);
611                                 bool was_empty = cur_candidates.empty();
612                                 if (ignore_case) {
613                                         dprintf("trigram %s (%zu bytes) decoded to %zu entries\n", print_trigram(trgm).c_str(), len, num);
614                                 }
615
616                                 for (TrigramDisjunction *td : uses_trigram[trgm]) {
617                                         done |= new_posting_list_read(td, decoded, &cur_candidates, &tmp);
618                                         if (done)
619                                                 break;
620                                 }
621                                 if (!ignore_case) {
622                                         if (was_empty) {
623                                                 dprintf("trigram %s (%zu bytes) decoded to %zu entries\n", print_trigram(trgm).c_str(), len, num);
624                                         } else if (cur_candidates.empty()) {
625                                                 dprintf("trigram %s (%zu bytes) decoded to %zu entries (none left, search is done)\n", print_trigram(trgm).c_str(), len, num);
626                                         } else {
627                                                 dprintf("trigram %s (%zu bytes) decoded to %zu entries (%zu left)\n", print_trigram(trgm).c_str(), len, num, cur_candidates.size());
628                                         }
629                                 }
630                         });
631                 }
632         }
633         engine.finish();
634         if (done) {
635                 return;
636         }
637         dprintf("Intersection done after %.1f ms. Doing final verification and printing:\n",
638                 1e3 * duration<float>(steady_clock::now() - start).count());
639
640         uint64_t matched = scan_docids(needles, cur_candidates, corpus, &engine);
641         dprintf("Done in %.1f ms, found %" PRId64 " matches.\n",
642                 1e3 * duration<float>(steady_clock::now() - start).count(), matched);
643
644         if (only_count) {
645                 printf("%" PRId64 "\n", matched);
646         }
647 }
648
649 string unescape_glob_to_plain_string(const string &needle)
650 {
651         string unescaped;
652         for (size_t i = 0; i < needle.size(); i += read_unigram(needle, i).second) {
653                 uint32_t ch = read_unigram(needle, i).first;
654                 assert(ch != WILDCARD_UNIGRAM);
655                 if (ch == PREMATURE_END_UNIGRAM) {
656                         fprintf(stderr, "Pattern '%s' ended prematurely\n", needle.c_str());
657                         exit(1);
658                 }
659                 unescaped.push_back(ch);
660         }
661         return unescaped;
662 }
663
664 regex_t compile_regex(const string &needle)
665 {
666         regex_t re;
667         int flags = REG_NOSUB;
668         if (ignore_case) {
669                 flags |= REG_ICASE;
670         }
671         if (use_extended_regex) {
672                 flags |= REG_EXTENDED;
673         }
674         int err = regcomp(&re, needle.c_str(), flags);
675         if (err != 0) {
676                 char errbuf[256];
677                 regerror(err, &re, errbuf, sizeof(errbuf));
678                 fprintf(stderr, "Error when compiling regex '%s': %s\n", needle.c_str(), errbuf);
679                 exit(1);
680         }
681         return re;
682 }
683
684 void usage()
685 {
686         printf(
687                 "Usage: plocate [OPTION]... PATTERN...\n"
688                 "\n"
689                 "  -c, --count            print number of matches instead of the matches\n"
690                 "  -d, --database DBPATH  search for files in DBPATH\n"
691                 "                         (default is " DEFAULT_DBPATH ")\n"
692                 "  -i, --ignore-case      search case-insensitively\n"
693                 "  -l, --limit LIMIT      stop after LIMIT matches\n"
694                 "  -0, --null             delimit matches by NUL instead of newline\n"
695                 "  -r, --regexp           interpret patterns as basic regexps (slow)\n"
696                 "      --regex            interpret patterns as extended regexps (slow)\n"
697                 "      --help             print this help\n"
698                 "      --version          print version information\n");
699 }
700
701 void version()
702 {
703         printf("plocate %s\n", PLOCATE_VERSION);
704         printf("Copyright 2020 Steinar H. Gunderson\n");
705         printf("License GPLv2+: GNU GPL version 2 or later <https://gnu.org/licenses/gpl.html>.\n");
706         printf("This is free software: you are free to change and redistribute it.\n");
707         printf("There is NO WARRANTY, to the extent permitted by law.\n");
708         exit(0);
709 }
710
711 int main(int argc, char **argv)
712 {
713         constexpr int EXTENDED_REGEX = 1000;
714         static const struct option long_options[] = {
715                 { "help", no_argument, 0, 'h' },
716                 { "count", no_argument, 0, 'c' },
717                 { "database", required_argument, 0, 'd' },
718                 { "ignore-case", no_argument, 0, 'i' },
719                 { "limit", required_argument, 0, 'l' },
720                 { "null", no_argument, 0, '0' },
721                 { "version", no_argument, 0, 'V' },
722                 { "regexp", no_argument, 0, 'r' },
723                 { "regex", no_argument, 0, EXTENDED_REGEX },
724                 { "debug", no_argument, 0, 'D' },  // Not documented.
725                 { 0, 0, 0, 0 }
726         };
727
728         setlocale(LC_ALL, "");
729         for (;;) {
730                 int option_index = 0;
731                 int c = getopt_long(argc, argv, "cd:hil:n:0VD", long_options, &option_index);
732                 if (c == -1) {
733                         break;
734                 }
735                 switch (c) {
736                 case 'c':
737                         only_count = true;
738                         break;
739                 case 'd':
740                         dbpath = strdup(optarg);
741                         break;
742                 case 'h':
743                         usage();
744                         exit(0);
745                 case 'i':
746                         ignore_case = true;
747                         break;
748                 case 'l':
749                 case 'n':
750                         limit_matches = limit_left = atoll(optarg);
751                         if (limit_matches <= 0) {
752                                 fprintf(stderr, "Error: limit must be a strictly positive number.\n");
753                                 exit(1);
754                         }
755                         break;
756                 case '0':
757                         print_nul = true;
758                         break;
759                 case 'r':
760                         patterns_are_regex = true;
761                         break;
762                 case EXTENDED_REGEX:
763                         patterns_are_regex = true;
764                         use_extended_regex = true;
765                         break;
766                 case 'D':
767                         use_debug = true;
768                         break;
769                 case 'V':
770                         version();
771                         break;
772                 default:
773                         exit(1);
774                 }
775         }
776
777         if (use_debug) {
778                 // Debug information would leak information about which files exist,
779                 // so drop setgid before we open the file; one would either need to run
780                 // as root, or use a locally-built file.
781                 if (setgid(getgid()) != 0) {
782                         perror("setgid");
783                         exit(EXIT_FAILURE);
784                 }
785         }
786
787         vector<Needle> needles;
788         for (int i = optind; i < argc; ++i) {
789                 Needle needle;
790                 needle.str = argv[i];
791
792                 // See if there are any wildcard characters, which indicates we should treat it
793                 // as an (anchored) glob.
794                 bool any_wildcard = false;
795                 for (size_t i = 0; i < needle.str.size(); i += read_unigram(needle.str, i).second) {
796                         if (read_unigram(needle.str, i).first == WILDCARD_UNIGRAM) {
797                                 any_wildcard = true;
798                                 break;
799                         }
800                 }
801
802                 if (patterns_are_regex) {
803                         needle.type = Needle::REGEX;
804                         needle.re = compile_regex(needle.str);
805                 } else if (any_wildcard) {
806                         needle.type = Needle::GLOB;
807                 } else if (ignore_case) {
808                         // strcasestr() doesn't handle locales correctly (even though LSB
809                         // claims it should), but somehow, fnmatch() does, and it's about
810                         // the same speed as using a regex.
811                         needle.type = Needle::GLOB;
812                         needle.str = "*" + needle.str + "*";
813                 } else {
814                         needle.type = Needle::STRSTR;
815                         needle.str = unescape_glob_to_plain_string(needle.str);
816                 }
817                 needles.push_back(move(needle));
818         }
819         if (needles.empty()) {
820                 fprintf(stderr, "plocate: no pattern to search for specified\n");
821                 exit(0);
822         }
823         do_search_file(needles, dbpath);
824 }