extern/kalker
1
0
Fork 0
mirror of https://github.com/PaddiM8/kalker.git synced 2024-12-13 01:50:41 +01:00
Code Issues Packages Projects Releases Wiki Activity

Integration with complex numbers

Browse Source
This commit is contained in:
bakk 2021-05-23 02:21:35 +02:00
parent ab6e758726
commit b0a1d3581b
2 changed files with 30 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

39
kalk/src/calculus.rs
View File

@ -80,26 +80,45 @@ fn simpsons_rule(
let mut result = KalkNum::default(); let mut result = KalkNum::default();
const N: i32 = 900; const N: i32 = 900;
let a = interpreter::eval_expr(context, a_expr, "")?.to_f64(); let a = interpreter::eval_expr(context, a_expr, "")?;
let b = interpreter::eval_expr(context, b_expr, "")?.to_f64(); let b = interpreter::eval_expr(context, b_expr, "")?;
let h = (b - a) / N as f64; let h = (b.sub_without_unit(a.clone())).div_without_unit(KalkNum::from(N));
for i in 0..=N { for i in 0..=N {
let variable_value = a
.clone()
.add_without_unit(KalkNum::from(i).mul_without_unit(h.clone()));
context.symbol_table.set(Stmt::VarDecl( context.symbol_table.set(Stmt::VarDecl(
Identifier::from_full_name(integration_variable), Identifier::from_full_name(integration_variable),
Box::new(Expr::Literal(a + i as f64 * h)), if variable_value.has_imaginary() {
Box::new(Expr::Binary(
Box::new(Expr::Literal(variable_value.to_f64())),
TokenKind::Plus,
Box::new(Expr::Binary(
Box::new(Expr::Literal(variable_value.imaginary_to_f64())),
TokenKind::Star,
Box::new(Expr::Var(Identifier::from_full_name("i"))),
)),
))
} else {
Box::new(Expr::Literal(variable_value.to_f64()))
},
)); ));
let factor = match i { let factor = KalkNum::from(match i {
0 | N => 1, 0 | N => 1,
_ if i % 3 == 0 => 2, _ if i % 3 == 0 => 2,
_ => 3, _ => 3,
}; });
// factor * f(x_n) // factor * f(x_n)
result.value += factor as f64 * interpreter::eval_expr(context, expr, "")?.value; let mul = factor.mul_without_unit(interpreter::eval_expr(context, expr, "")?);
result.value += mul.value;
result.imaginary_value += mul.imaginary_value;
} }
result.value *= (3f64 * h) / 8f64; Ok(result.mul_without_unit(KalkNum::new_with_imaginary(
3f64 / 8f64 * h.value,
Ok(result) &h.unit,
3f64 / 8f64 * h.imaginary_value,
)))
} }

2
kalk/src/parser.rs
View File

@ -111,7 +111,7 @@ pub enum CalcError {
impl ToString for CalcError { impl ToString for CalcError {
fn to_string(&self) -> String { fn to_string(&self) -> String {
match self match self {
CalcError::ExpectedDx => format!("Expected eg. dx, to specify for which variable the operation is being done to. Example with integration: ∫(0, 1, x dx) or ∫(0, 1, x, dx). You may need to put parenthesis around the expression before dx/dy/du/etc."), CalcError::ExpectedDx => format!("Expected eg. dx, to specify for which variable the operation is being done to. Example with integration: ∫(0, 1, x dx) or ∫(0, 1, x, dx). You may need to put parenthesis around the expression before dx/dy/du/etc."),
CalcError::IncorrectAmountOfArguments(expected, func, got) => format!( CalcError::IncorrectAmountOfArguments(expected, func, got) => format!(
"Expected {} arguments for function {}, but got {}.", "Expected {} arguments for function {}, but got {}.",