Hello, i have tried . Here, is my code>
real z = 5; // Stretching factor for exponential density
// Exponentially stretched borders: fine near ends (i.e., corners)
border bottom(t=0,1){ x=(exp(z*t)-1)/(exp(z)-1); y=0; label=1; };
border right(t=0,1){ x=1; y=(exp(z*t)-1)/(exp(z)-1); label=2; };
border top(t=1,0){ x=(exp(z*t)-1)/(exp(z)-1); y=1; label=3; };
border left(t=1,0){ x=0; y=(exp(z*t)-1)/(exp(z)-1); label=4; };
int n = 10; // Higher n = finer mesh
mesh Th = buildmesh(
bottom(n*4)
+ right(n*10)
+ top(n*2)
+ left(n*5)
);
plot(Th, wait=1, cmm="Uniform boundary-refined mesh");
savemesh(Th, "boundary_refined.msh");
// Refine only near the boundary (optional)
real tol = 0.05;
func isBoundaryRefined = (x < tol) || (x > 1 - tol) || (y < tol) || (y > 1 - tol);
mesh Th2 = adaptmesh(Th, isBoundaryRefined, IsMetric=1, nbvx=100000, err=1e-3);
plot(Th2, wait=1, cmm="Boundary + corner-refined");
savemesh(Th2, "boundary_corner_refined.msh");
But exactly, i am not getting that i want. Please help.