Implement data-blind scalar quantization - #16030
Conversation
Allow callers to disable centering at the format level, which also disables writing of float vectors since they are no longer needed. Includes a path to handle of mix of centered and uncentered segments as input. In this case the uncentered/no float vectors will be dequantized and requantized but this case should be relatively uncommon. Includes OSQ changes to allow a zero vector for COSINE if the vector is not a unit vector. Maybe fix this in upstream callers?
|
This PR has not had activity in the past 2 weeks, labeling it as stale. If the PR is waiting for review, notify the dev@lucene.apache.org list. Thank you for your contribution! |
|
This PR has not had activity in the past 2 weeks, labeling it as stale. If the PR is waiting for review, notify the dev@lucene.apache.org list. Thank you for your contribution! |
|
Sorry for getting to this PR so late, but I realized that data-blind quantization would be super helpful for another use-case too: in an attempt to support multiple HNSW graphs backed by the same vector storage (#14758), we need to be able to effectively de-duplicate vectors across fields. Lucene currently supports a centroid-centered quantization, which might cause the same raw vector to be quantized differently across different fields. If it is data-blind (i.e. centered on the zero vector), quantized bytes can be sanely de-duplicated. See #16506. |
|
+1, I love that we are making progress on data-blind quantization. This should work well when incoming vectors are isotropic (all dimensions behave the same i.e. the histograms of their per-dimension values approximate little baby gaussians with mean 0) and variance as required to be on unit sphere at that dimensionality. But I think your average trained model won't just produce isotropic embeddings? For such cases (probably the common case? not sure), we have pre-conditioning / random Hadamard rotation (another PR in flight for this? -- yes #16092!) which should (usually? there are adversaries (intentional or otherwise) for any rotation matrix right?) scrub anisotropic vectors. Sort of like the record and record players in GEB -- thank you Kurt Gödel!). These are all experimental vector codecs but I hope they eventually become default. If we always pre-condition then we almost always can do data blind quantization that is just as good as non-data-blind quantization. Does this PR make any effort / at least javadocs to explain that you should ensure your incoming vectors are isotropic? Or to spot check if they really seem to be isotropic? luceneutil has all sorts of smell detection ("smelling pipeline" a recent genai model called it!) now to detect all sorts of problems your otherwise very-opaque-to-humans vectors might have. |
|
I'll try to revive this over the next week or two.
IIRC the base quantizer should behave pretty well if the distribution is uncentered as well (this is what the lower/upper interval are for and why the dot product is unsigned), but if it's not ~Gaussian I would expect the error rate to be high. RE: rotation -- I think we might want this to provide this but it should be layered just above Lucene. My concern is that if rotation is delegated to the segments during search costs will be very high -- I would expect O(1-2usecs) for a heavily optimized FWHT on good hardware, and you would have to multiply this cost by the number of segments. It'll look fine in a benchmark that's merged to one segment but poor in practice.
This is not documented, but is generally a constraint for most quantizers. It would be easy to document but I'm not sure what we would do if we detected your vectors don't quantize well. If we had an auto-quantization setting of some kind at that point you would just fall back to float32 or float16. In general I'm annoyed by the amount of code I have to copy here, especially since all of the flat formats on the read side are just a fixed stride index read. Maybe the layer of abstraction should be closer to the vector layer ( |
|
This PR has not had activity in the past 2 weeks, labeling it as stale. If the PR is waiting for review, notify the dev@lucene.apache.org list. Thank you for your contribution! |
|
Ok, I've revived this branch, renaming to Lucene106 and merging in main and all of the fixes that entails. PTAL. |
kaivalnp
left a comment
There was a problem hiding this comment.
Thanks @mccullocht, this enables some great improvements in indexing vectors! (faster merges / smaller indexes)
I didn't get through the whole PR, but leaving some initial thoughts..
| - 2 * score; | ||
| // Ensure that 'score' (the squared euclidean distance) is non-negative. The computed value | ||
| // may be negative as a result of quantization loss. | ||
| return 1 / (1f + Math.max(score, 0f)); |
There was a problem hiding this comment.
Can we use existing utility functions for this? Same for dot product.
There was a problem hiding this comment.
This file has been deleted.
|
|
||
| @Override | ||
| public FlatFieldVectorsWriter<?> addField(FieldInfo fieldInfo) throws IOException { | ||
| if (fieldInfo.getVectorEncoding().equals(VectorEncoding.FLOAT32)) { |
There was a problem hiding this comment.
#16473 recently added scalar-quantization support for FP16 vectors, should we include that too?
There was a problem hiding this comment.
Moving everything back into Lucene104 means we get this feature "for free"
| * <li><b>int</b> the field number | ||
| * <li><b>int</b> the vector encoding ordinal | ||
| * <li><b>int</b> the vector similarity ordinal | ||
| * <li><b>vint</b> the vector dimensions | ||
| * <li><b>vlong</b> the offset to the vector data in the .veq file | ||
| * <li><b>vlong</b> the length of the vector data in the .veq file | ||
| * <li><b>vint</b> the number of vectors | ||
| * <li><b>vint</b> the wire number for ScalarEncoding | ||
| * <li><b>[float]</b> the centroid | ||
| * <li><b>float</b> the centroid square magnitude | ||
| * <li>The sparse vector information, if required, mapping vector ordinal to doc ID |
There was a problem hiding this comment.
This appears to be the same as the Lucene104 codec?
Is there no change to the on-disk representation? If so, could / should we add the data-blind option to the Lucene104 class instead of creating a new format?
There was a problem hiding this comment.
This is largely the same. There's no guidance on this and I probably could merge it back into the Lucene104 codec as a version bump. I have mixed feelings about new version vs codec internal version bump because it's less obvious with the codec internal version bump that I am breaking you on upgrade.
There was a problem hiding this comment.
Pushed everything back into the Lucene104 codec. Disabling centering results in an internal VERSION bump that omits the centroid and derived data and replaces it with zeros if needed. Writers with centering enable (default) won't see anything different at all, writers with centering disabled will experience the usually one way door.
| */ | ||
|
|
||
| /** | ||
| * Lucene 10.5 scalar quantized vector format, extending 10.4 with data-blind mode ({@code |
| fi.isDataBlind() ? null : rawVectorsReader.getFloatVectorValues(field); | ||
|
|
||
| if (rawFloatVectorValues == null) { | ||
| // Data-blind mode: full-precision float vectors were never stored. Reconstruct floats from |
There was a problem hiding this comment.
[DISCUSS] I think Lucene should disallow operations on data it dropped (i.e. throw an error when trying to retrieve floats / rescore using floats that were not stored in the index, instead of using a lossy value).
Does Lucene have other use-cases that support indexing some data (and specific operations to search it), but not retain / support returning the original data?
One example: if a user wrote some segments using enableCentering = false, then re-opens the index and attempts to merge with enableCentering = true, the dequantized vectors would be used instead of the original ones, leading to further loss of information?
There was a problem hiding this comment.
There might be some use for this in some very high dimensional data sets -- I may want to quantize the original vectors at a high bit rate (4 or 8) bits and discard them, and also provide 1 or 2 bit quantization with centering that are used for graph navigation.
In your hypothetical situation yes you would absolutely lose precision, how much depends on how aggressive the initial quantization is.
There was a problem hiding this comment.
It's not part of the automated indexing flow but IIUC there is already support for truncating the float vectors that works in a comparable (but not identical) way -- if the float vector files mysteriously end up empty then we reconstruct the values from the quantized data. At 4 or 8 bits this is almost always pretty reasonable.
| if (centroid == null) return false; | ||
| for (float v : centroid) { | ||
| if (v != 0f) return false; | ||
| } |
There was a problem hiding this comment.
Q: what if the centroid happens to be zero, even if enableCentering = true was used?
Would this cause the returned FloatVectorValues to be built from dequantized bytes, even when the full floats were available?
Should enableCentering be persisted in the index, instead of inferring from the centroid?
There was a problem hiding this comment.
If the center is already a zero vector that is geometrically equivalent to our data blind configuration. Centering wouldn't subtract anything from the vector, and for angular similarity center_dot would be zero.
mccullocht
left a comment
There was a problem hiding this comment.
Hey @kaivalnp thanks for the review, I know this is very large. I will address the rest of your comments over the next couple of days.
| * <li><b>int</b> the field number | ||
| * <li><b>int</b> the vector encoding ordinal | ||
| * <li><b>int</b> the vector similarity ordinal | ||
| * <li><b>vint</b> the vector dimensions | ||
| * <li><b>vlong</b> the offset to the vector data in the .veq file | ||
| * <li><b>vlong</b> the length of the vector data in the .veq file | ||
| * <li><b>vint</b> the number of vectors | ||
| * <li><b>vint</b> the wire number for ScalarEncoding | ||
| * <li><b>[float]</b> the centroid | ||
| * <li><b>float</b> the centroid square magnitude | ||
| * <li>The sparse vector information, if required, mapping vector ordinal to doc ID |
There was a problem hiding this comment.
This is largely the same. There's no guidance on this and I probably could merge it back into the Lucene104 codec as a version bump. I have mixed feelings about new version vs codec internal version bump because it's less obvious with the codec internal version bump that I am breaking you on upgrade.
| fi.isDataBlind() ? null : rawVectorsReader.getFloatVectorValues(field); | ||
|
|
||
| if (rawFloatVectorValues == null) { | ||
| // Data-blind mode: full-precision float vectors were never stored. Reconstruct floats from |
There was a problem hiding this comment.
There might be some use for this in some very high dimensional data sets -- I may want to quantize the original vectors at a high bit rate (4 or 8) bits and discard them, and also provide 1 or 2 bit quantization with centering that are used for graph navigation.
In your hypothetical situation yes you would absolutely lose precision, how much depends on how aggressive the initial quantization is.
| if (centroid == null) return false; | ||
| for (float v : centroid) { | ||
| if (v != 0f) return false; | ||
| } |
There was a problem hiding this comment.
If the center is already a zero vector that is geometrically equivalent to our data blind configuration. Centering wouldn't subtract anything from the vector, and for angular similarity center_dot would be zero.
Add an option to the quantization format to enable or disable centering (enabled by default). When centering is disabled we also stop writing the float vectors which can lead to significant storage savings. Special handling is included during merges -- we check that all of the input is in the same encoding, and handle transcoding if some of the input is float vectors.
Large portions of this change were generated using claude code. I reviewed, tweaked, and tested the code before putting it up for review.
This change appears as a new parameter on the existing
Lucene104ScalarQuantizedVectorsFormatto enable centering (default is true). Disabling centering causes an internal format version bump triggering the new behavior -- disabling computation and storage of the centroid vector and related parameters.luceneutil results -- 1M cohere vectors, 8 bit quantization.
before:
after
The harness extrapolates vector size from the input size so believe the on-disk index_size number -- this is about 4x smaller. Force merge is faster since we don't have to re-quantize vectors on merge. Recall is very similar but YMMV.
See also #16029