Mercurial > repos > public > sbplib

comparison +multiblock/+domain/Annulus.m @ 1257:4d472d020ccf feature/FMMlabb

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Assert that the inner radius is positive
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 12 Feb 2020 13:27:27 +0100
parents 479d5672fecd
children 09a5783a3d37
comparison
equal deleted inserted replaced
1256:479d5672fecd 1257:4d472d020ccf
12 default_arg('c_inner', [0; 0]); 12 default_arg('c_inner', [0; 0]);
13 default_arg('r_outer', 1) 13 default_arg('r_outer', 1)
14 default_arg('c_outer', [0; 0]); 14 default_arg('c_outer', [0; 0]);
15 % Assert that the problem is well-defined 15 % Assert that the problem is well-defined
16 d = norm(c_outer-c_inner,2); 16 d = norm(c_outer-c_inner,2);
17 assert(r_inner > 0, 'Inner radius must be greater than zero');
17 assert(r_outer > d+r_inner, 'Inner disk not contained in outer disk'); 18 assert(r_outer > d+r_inner, 'Inner disk not contained in outer disk');
18 19
19 cir_out_A = parametrization.Curve.circle(c_outer,r_outer,[-pi/2 pi/2]); 20 cir_out_A = parametrization.Curve.circle(c_outer,r_outer,[-pi/2 pi/2]);
20 cir_in_A = parametrization.Curve.circle(c_inner,r_inner,[pi/2 -pi/2]); 21 cir_in_A = parametrization.Curve.circle(c_inner,r_inner,[pi/2 -pi/2]);
21 22