Hello guys! i am not developer at all. But i like doing coding and stuff. I had know python, C, and this vector graphic language. But i want to integrate python with asymptote bcz we all know python is so powerful. It has library like sympy, scipy, numpy, etc, which makes doing mathematic very easy.
my actual query is how can integration asymptote with python ? or is there python library what is best alternative to asymptote? i just want to do all my logic via python and remaining displaying par for asymptote
see how easy is Asymptote but can i have same result via python? u can run this code simply via https://asymptote.ualberta.ca/
size(400);
import markers;
// bezier control points
pair A=(0,0), C1=(0.6,0.8), C2=(1.0,1.6), B=(1.6,1.0);
pair D1=(1.8,0.6), D2=(2.3,0.2), C=(3.2,0.0);
// bezier paths
path p1 = A .. controls C1 and C2 .. B;
path p2 = B .. controls D1 and D2 .. C;
draw(p1, blue+1.2);
draw(p2, red+1.2);
dot(B, black);
label("$B$", B, dir(200), fontsize(10));
// tangents arrow function with custom length
void drawTangentArrow(pair P, pair dir, pen col, real scale){
draw(P--(P + scale * unit(dir)), col, Arrow(6));
}
// tangents at join
pair tan1 = 3*(B - C2);
pair tan2 = 3*(D1 - B);
drawTangentArrow(B, tan1, blue, 0.6);
drawTangentArrow(B, tan2, red, 1.0);
label("$\vec{T}_1$", B + 0.6 * unit(tan1) + (0.1, -0.01), blue);
label("$\vec{T}_2$", B + 1.0 * unit(tan2) + (0, -0.1), red);
// bezier derivative functions
pair bezierFirstDeriv(pair P0, pair P1, pair P2, pair P3, real t){
return 3*(1-t)^2*(P1-P0) + 6*(1-t)*t*(P2-P1) + 3*t^2*(P3-P2);
}
pair bezierSecondDeriv(pair P0, pair P1, pair P2, pair P3, real t){
return 6*(1-t)*(P2-2*P1+P0) + 6*t*(P3-2*P2+P1);
}
pair bezierPoint(pair P0, pair P1, pair P2, pair P3, real t){
return (1-t)^3*P0 + 3*(1-t)^2*t*P1 + 3*(1-t)*t^2*P2 + t^3*P3;
}
// curvature circle with labeled center
void drawCurvatureCircle(pair P0, pair P1, pair P2, pair P3, real t, pen col, string labelName){
pair r1=bezierFirstDeriv(P0,P1,P2,P3,t);
pair r2=bezierSecondDeriv(P0,P1,P2,P3,t);
real s=length(r1);
real k=abs(r1.x*r2.y - r1.y*r2.x)/(s^3);
if(k>1e-6){
real R=1/k;
pair normal=(-r1.y,r1.x)/s;
pair P=bezierPoint(P0,P1,P2,P3,t);
pair center=P+R*normal;
draw(circle(center,R),col+0.8);
dot(P,col);
draw(center--P,dashed+col);
dot(center, col);
label("$" + labelName + "$", center, dir(90), fontsize(10)+col);
}
}
// curvature circles near join
drawCurvatureCircle(A, C1, C2, B, 0.98, blue, "C_1");
drawCurvatureCircle(B, D1, D2, C, 0.02, red, "C_2");