/*
 * Licensed to the Apache Software Foundation (ASF) under one
 * or more contributor license agreements.  See the NOTICE file
 * distributed with this work for additional information
 * regarding copyright ownership.  The ASF licenses this file
 * to you under the Apache License, Version 2.0 (the
 * "License"); you may not use this file except in compliance
 * with the License.  You may obtain a copy of the License at
 *
 *   http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing,
 * software distributed under the License is distributed on an
 * "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY
 * KIND, either express or implied.  See the License for the
 * specific language governing permissions and limitations
 * under the License.
 */
package hivemall.knn.similarity;

import hivemall.model.FeatureValue;
import hivemall.utils.hadoop.HiveUtils;

import java.util.Arrays;
import java.util.HashMap;
import java.util.List;
import java.util.Map;

import org.apache.hadoop.hive.ql.exec.Description;
import org.apache.hadoop.hive.ql.exec.UDFArgumentException;
import org.apache.hadoop.hive.ql.metadata.HiveException;
import org.apache.hadoop.hive.ql.udf.UDFType;
import org.apache.hadoop.hive.ql.udf.generic.GenericUDF;
import org.apache.hadoop.hive.serde2.objectinspector.ListObjectInspector;
import org.apache.hadoop.hive.serde2.objectinspector.ObjectInspector;
import org.apache.hadoop.hive.serde2.objectinspector.primitive.PrimitiveObjectInspectorFactory;
import org.apache.hadoop.io.FloatWritable;

//@formatter:off
@Description(name = "cosine_similarity",
        value = "_FUNC_(ftvec1, ftvec2) - Returns a cosine similarity of the given two vectors",
        extended = "WITH docs as (\n" + 
                "  select 1 as docid, array('apple:1.0', 'orange:2.0', 'banana:1.0', 'kuwi:0') as features\n" + 
                "  union all\n" + 
                "  select 2 as docid, array('apple:1.0', 'orange:0', 'banana:2.0', 'kuwi:1.0') as features\n" + 
                "  union all\n" + 
                "  select 3 as docid, array('apple:2.0', 'orange:0', 'banana:2.0', 'kuwi:1.0') as features\n" + 
                ") \n" + 
                "select\n" + 
                "  l.docid as doc1,\n" + 
                "  r.docid as doc2,\n" + 
                "  cosine_similarity(l.features, r.features) as similarity\n" + 
                "from \n" + 
                "  docs l\n" + 
                "  CROSS JOIN docs r\n" + 
                "where\n" + 
                "  l.docid != r.docid\n" + 
                "order by \n" + 
                "  doc1 asc,\n" + 
                "  similarity desc;\n" + 
                "\n" + 
                "doc1    doc2    similarity\n" + 
                "1       3       0.5443311\n" + 
                "1       2       0.5\n" + 
                "2       3       0.9525793\n" + 
                "2       1       0.5\n" + 
                "3       2       0.9525793\n" + 
                "3       1       0.5443311")
//@formatter:on
@UDFType(deterministic = true, stateful = false)
public final class CosineSimilarityUDF extends GenericUDF {

    private ListObjectInspector arg0ListOI, arg1ListOI;

    @Override
    public ObjectInspector initialize(ObjectInspector[] argOIs) throws UDFArgumentException {
        if (argOIs.length != 2) {
            throw new UDFArgumentException("cosine_similarity takes 2 arguments");
        }
        this.arg0ListOI = HiveUtils.asListOI(argOIs, 0);
        this.arg1ListOI = HiveUtils.asListOI(argOIs, 1);

        return PrimitiveObjectInspectorFactory.writableFloatObjectInspector;
    }

    @Override
    public FloatWritable evaluate(DeferredObject[] arguments) throws HiveException {
        List<String> ftvec1 = HiveUtils.asStringList(arguments[0], arg0ListOI);
        List<String> ftvec2 = HiveUtils.asStringList(arguments[1], arg1ListOI);
        float similarity = cosineSimilarity(ftvec1, ftvec2);
        return new FloatWritable(similarity);
    }

    public static float cosineSimilarity(final List<String> ftvec1, final List<String> ftvec2) {
        if (ftvec1 == null || ftvec2 == null) {
            return 0.f;
        }

        final FeatureValue probe = new FeatureValue();
        final Map<String, Float> map1 = new HashMap<String, Float>(ftvec1.size() * 2 + 1);
        double score1 = 0.d;
        for (String ft : ftvec1) {
            FeatureValue.parseFeatureAsString(ft, probe);
            float v = probe.getValueAsFloat();
            score1 += (v * v);
            String f = probe.getFeature();
            map1.put(f, v);
        }
        double l1norm1 = Math.sqrt(score1);

        float dotp = 0.f;
        double score2 = 0.d;
        for (String ft : ftvec2) {
            FeatureValue.parseFeatureAsString(ft, probe);
            float v2 = probe.getValueAsFloat();
            score2 += (v2 * v2);
            String f2 = probe.getFeature();
            Float v1 = map1.get(f2);
            if (v1 != null) {
                dotp += (v1.floatValue() * v2);
            }
        }
        double l1norm2 = Math.sqrt(score2);

        final double denom = l1norm1 * l1norm2;
        if (denom <= 0.f) {
            return 0.f;
        } else {
            return (float) (dotp / denom);
        }
    }

    @Override
    public String getDisplayString(String[] children) {
        return "cosine_similarity(" + Arrays.toString(children) + ")";
    }

}