Fcalva/chibi-scheme
1
0
Fork 0
mirror of https://github.com/ashinn/chibi-scheme.git synced 2025-05-19 21:59:17 +02:00
Code Issues Projects Releases Packages Wiki Activity Actions

breaking down and adding polar notation

Browse source
This commit is contained in:
Alex Shinn 2012-07-10 23:19:57 +09:00
parent c063dadd2e
commit 2d62671df8
1 changed files with 26 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

27
sexp.c
View file

@ -2207,6 +2207,27 @@ sexp sexp_read_complex_tail (sexp ctx, sexp in, sexp real) {
sexp_gc_release1(ctx);
return sexp_complex_normalize(res);
}
#if SEXP_USE_MATH
sexp sexp_read_polar_tail (sexp ctx, sexp in, sexp magnitude) {
sexp_gc_var2(res, theta);
sexp_gc_preserve2(ctx, res, theta);
theta = sexp_read_number(ctx, in, 10);
if (sexp_exceptionp(theta)) {
res = theta;
} else if (sexp_complexp(theta) || !sexp_numberp(theta)) {
res = sexp_read_error(ctx, "invalid polar numeric syntax", theta, in);
} else {
res = sexp_make_complex(ctx, SEXP_ZERO, SEXP_ZERO);
sexp_complex_real(res) = sexp_cos(ctx, NULL, 1, theta);
sexp_complex_real(res) = sexp_mul(ctx, magnitude, sexp_complex_real(res));
sexp_complex_imag(res) = sexp_sin(ctx, NULL, 1, theta);
sexp_complex_imag(res) = sexp_mul(ctx, magnitude, sexp_complex_imag(res));
}
sexp_gc_release2(ctx);
return sexp_complex_normalize(res);
}
#endif
#endif
sexp sexp_read_float_tail (sexp ctx, sexp in, double whole, int negp) {
@ -2386,9 +2407,13 @@ sexp sexp_read_number (sexp ctx, sexp in, int base) {
sexp_gc_release2(ctx);
return res;
#if SEXP_USE_COMPLEX
} else if (c=='i' || c=='I' || c=='+' || c=='-') {
} else if (c=='i' || c=='I' || c=='+' || c=='-' || c=='@') {
if (base != 10)
return sexp_read_error(ctx, "found non-base 10 complex", SEXP_NULL, in);
#if SEXP_USE_MATH
if (c=='@')
return sexp_read_polar_tail(ctx, in, sexp_make_fixnum(negativep ? -val : val));
#endif
sexp_push_char(ctx, c, in);
return sexp_read_complex_tail(ctx, in, sexp_make_fixnum(negativep ? -val : val));
#endif