# -*- coding: utf-8 -*-
"""原创离线合成样例；不连接模型、不读取账户。输出位于脚本旁。
已有输入保留，修改输入后可以重跑。失败边界由 --negative 单独演示。
"""
from pathlib import Path
import sys, json, csv, math, io, hashlib, html, base64, re
sys.stdout.reconfigure(encoding='utf-8')
sys.stderr.reconfigure(encoding='utf-8')
BASE = Path(__file__).resolve().parent
NEG = '--negative' in sys.argv
def write(name, value):
    p=BASE/name; p.parent.mkdir(parents=True,exist_ok=True)
    if isinstance(value,bytes): p.write_bytes(value)
    else: p.write_text(value,encoding='utf-8',newline='')
def seed(name,value):
    if not (BASE/name).exists(): write(name,value)
def jwrite(name,value): write(name,json.dumps(value,ensure_ascii=False,indent=2))
def jseed(name,value):
    if not (BASE/name).exists(): jwrite(name,value)
def load(name): return json.loads((BASE/name).read_text(encoding='utf-8'))
def table(name,heads,rows):
    f=io.StringIO(newline='');w=csv.writer(f);w.writerow(heads);w.writerows(rows);write(name,f.getvalue())
def page(name,title,body):
    write(name,'<!doctype html><html lang="zh-CN"><meta charset="utf-8"><meta name="viewport" content="width=device-width,initial-scale=1"><title>'+html.escape(title)+'</title><style>body{max-width:1000px;margin:2rem auto;padding:0 1rem;font:17px/1.7 system-ui}table{border-collapse:collapse}td,th{border:1px solid #aaa;padding:.5rem}pre{white-space:pre-wrap;overflow-wrap:anywhere}img,svg{max-width:100%;height:auto}a{color:#125caa}</style><h1>'+html.escape(title)+'</h1>'+body+'</html>')
def emit(data): print(json.dumps(data,ensure_ascii=False))

seed('vectors.csv','id,x,y\nq,1,0\na,10,1\nb,.7,.7\n')
vs={r['id']:[float(r['x']),float(r['y'])] for r in csv.DictReader((BASE/'vectors.csv').open())};q=vs.pop('q')
if NEG: vs['zero']=[0,0]
out=[]
for key,v in vs.items():
    nq=math.sqrt(sum(x*x for x in q));nv=math.sqrt(sum(x*x for x in v))
    if nq==0 or nv==0: raise ValueError('零向量不能求余弦')
    cosine=sum(x*y for x,y in zip(q,v))/(nq*nv);distance=math.dist(q,v);dot=sum(x*y for x,y in zip(q,v));normalized_sq=sum((x/nq-y/nv)**2 for x,y in zip(q,v));out.append([key,cosine,distance,dot,normalized_sq,2*(1-cosine)])
table('distance-rankings.csv',['id','cosine','euclidean','dot','normalized_squared_l2','two_one_minus_cos'],out)
emit({'cosine_desc':[r[0] for r in sorted(out,key=lambda r:-r[1])],'euclidean_asc':[r[0] for r in sorted(out,key=lambda r:r[2])],'rows':out})
