计算器 abacus 2 发布

[zhcosin]

计算器 abacus 2 今日发布了，下载网址 http://www.oschina.net/code/snippet_736932_13725，技术文档网址：http://my.oschina.net/zhcosin/blog?catalog=254400， abacus 是一个简易的、可扩展的、开源的小型计算器，支持四则混合运算和众多的数学函数，如三角函数和组合数等，支持符号常量，第2版又添加了对用户自定义函数的支持，对变量参与运算的支持也正在紧锣密鼓的开发中。
详情请参阅：http://my.oschina.net/zhcosin/blog/101082

[赵4老师]

仅供参考：

/*---------------------------------------
函数型计算器(VC++6.0,Win32 Console)程序由 yu_hua 于2007-07-27设计完成
功能：
目前提供了10多个常用数学函数:
    ⑴正弦sin
    ⑵余弦cos
    ⑶正切tan
    ⑷开平方sqrt
    ⑸反正弦arcsin
    ⑹反余弦arccos
    ⑺反正切arctan
    ⑻常用对数lg
    ⑼自然对数ln
    ⑽ｅ指数exp
    ⑾乘幂函数∧
用法：
如果要求2的32次幂，可以打入2^32<回车>
如果要求30度角的正切可键入tan(Pi/6)<回车>
注意不能打入：tan(30)<Enter>
如果要求1.23弧度的正弦，有几种方法都有效：
sin(1.23)<Enter>
sin 1.23 <Enter>
sin1.23  <Enter>
如果验证正余弦的平方和公式,可打入sin(1.23)^2+cos(1.23)^2 <Enter>或sin1.23^2+cos1.23^2 <Enter>
此外两函数表达式连在一起,自动理解为相乘如：sin1.23cos0.77+cos1.23sin0.77就等价于sin(1.23)*cos(0.77)+cos(1.23)*sin(0.77)
当然你还可以依据三角变换，再用sin(1.23+0.77)也即sin2验证一下。
本计算器充分考虑了运算符的优先级因此诸如：2+3*4^2 实际上相当于：2+(3*(4*4))
另外函数名前面如果是数字,那么自动认为二者相乘.
同理，如果某数的右侧是左括号，则自动认为该数与括弧项之间隐含一乘号。
如：3sin1.2^2+5cos2.1^2 相当于3*sin2(1.2)+5*cos2(2.1)
又如：4(3-2(sqrt5-1)+ln2)+lg5 相当于4*(3-2*(√5 -1)+loge(2))+log10(5)
此外，本计算器提供了圆周率 Pi键入字母时不区分大小写,以方便使用。
----------------------------------------*/
#include <iostream>
#include <iomanip>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <cmath>
#include <stdio.h>
#include <string.h>
#include <windows.h>
using namespace std;
const char Tab=0x9;
const int  DIGIT=1;
const int MAXLEN=16384;
char s[MAXLEN],*endss;
int pcs=15;
double fun(double x,char op[],int *iop) {
    while (op[*iop-1]<32) //本行使得函数嵌套调用时不必加括号,如 arc sin(sin(1.234)) 只需键入arc sin sin 1.234<Enter>
        switch (op[*iop-1]) {
        case  7: x=sin(x);  (*iop)--;break;
        case  8: x=cos(x);  (*iop)--;break;
        case  9: x=tan(x);  (*iop)--;break;
        case 10: x=sqrt(x); (*iop)--;break;
        case 11: x=asin(x); (*iop)--;break;
        case 12: x=acos(x); (*iop)--;break;
        case 13: x=atan(x); (*iop)--;break;
        case 14: x=log10(x);(*iop)--;break;
        case 15: x=log(x);  (*iop)--;break;
        case 16: x=exp(x);  (*iop)--;break;
        }
    return x;
}
double calc(char *expr,char **addr) {
    static deep; //递归深度
    static char *fname[]={ "sin","cos","tan","sqrt","arcsin","arccos","arctan","lg","ln","exp",NULL};
    double ST[10]={0.0}; //数字栈
    char op[10]={'+'}; //运算符栈
    char c,*rexp,*pp,*pf;
    int ist=1,iop=1,last;
    if (!deep) {
        pp=pf=expr;
        do {
            c = *pp++;
            if (c!=' '&& c!=Tab)
                *pf++ = c;
        } while (c!='\0');
    }
    pp=expr;
    if ((c=*pp)=='-'||c=='+') {
        op[0] = c;
        pp++;
    }
    last = !DIGIT;
    while ((c=*pp)!='\0') {
        if (c=='(') {//左圆括弧
            deep++;
            ST[ist++]=calc(++pp,addr);
            deep--;
            ST[ist-1]=fun(ST[ist-1],op,&iop);
            pp = *addr;
            last = DIGIT;
            if (*pp == '('||isalpha(*pp) && strnicmp(pp,"Pi",2)) {//目的是：当右圆括弧的右恻为左圆括弧或函数名字时，默认其为乘法
                op[iop++]='*';
                last = !DIGIT;
                c = op[--iop];
                goto operate ;
            }
        }
        else if (c==')') {//右圆括弧
            pp++;
            break;
        } else if (isalpha(c)) {
            if (!strnicmp(pp,"Pi",2)) {
                if (last==DIGIT) {
                    cout<< "π左侧遇）" <<endl;exit(1);
                }
                ST[ist++]=3.14159265358979323846264338328;
                ST[ist-1]=fun(ST[ist-1],op,&iop);
                pp += 2;
                last = DIGIT;
                if (!strnicmp(pp,"Pi",2)) {
                    cout<< "两个π相连" <<endl;exit(2);
                }
                if (*pp=='(') {
                    cout<< "π右侧遇（" <<endl;exit(3);
                }
            } else {
                for(int i=0; (pf=fname[i])!=NULL; i++)
                    if (!strnicmp(pp,pf,strlen(pf)))break;
                if (pf!=NULL) {
                    op[iop++] = 07+i;
                    pp += strlen(pf);
                } else {
                    cout<< "陌生函数名" <<endl;exit(4);
                }
            }
        } else if (c=='+'||c=='-'||c=='*'||c=='/'||c=='^') {
            char cc;
            if (last != DIGIT) {
                cout<< "运算符粘连" <<endl;exit(5);
            }
            pp++;
            if (c=='+'||c=='-') {
                do {
                    cc = op[--iop];
                    --ist;
                    switch (cc) {
                    case '+':  ST[ist-1] += ST[ist];break;
                    case '-':  ST[ist-1] -= ST[ist];break;
                    case '*':  ST[ist-1] *= ST[ist];break;
                    case '/':  ST[ist-1] /= ST[ist];break;
                    case '^':  ST[ist-1] = pow(ST[ist-1],ST[ist]);break;
                    }
                } while (iop);
                op[iop++] = c;
            } else if (c=='*'||c=='/') {
operate:        cc = op[iop-1];
                if (cc=='+'||cc=='-') {
                    op[iop++] = c;
                } else {
                    --ist;
                    op[iop-1] = c;
                    switch (cc) {
                    case '*':  ST[ist-1] *= ST[ist];break;
                    case '/':  ST[ist-1] /= ST[ist];break;
                    case '^':  ST[ist-1] = pow(ST[ist-1],ST[ist]);break;
                    }
                }
            } else {
                cc = op[iop-1];
                if (cc=='^') {
                    cout<< "乘幂符连用" <<endl;exit(6);
                }
                op[iop++] = c;
            }
            last = !DIGIT;
        } else {
            if (last == DIGIT) {
                cout<< "两数字粘连" <<endl;exit(7);
            }
            ST[ist++]=strtod(pp,&rexp);
            ST[ist-1]=fun(ST[ist-1],op,&iop);
            if (pp == rexp) {
                cout<< "非法字符" <<endl;exit(8);
            }
            pp = rexp;
            last = DIGIT;
            if (*pp == '('||isalpha(*pp)) {
                op[iop++]='*';
                last = !DIGIT;
                c = op[--iop];
                goto operate ;
            }
        }
    }
    *addr=pp;
    if (iop>=ist) {
        cout<< "表达式有误" <<endl;exit(9);
    }
    while (iop) {
        --ist;
        switch (op[--iop]) {
        case '+':  ST[ist-1] += ST[ist];break;
        case '-':  ST[ist-1] -= ST[ist];break;
        case '*':  ST[ist-1] *= ST[ist];break;
        case '/':  ST[ist-1] /= ST[ist];break;
        case '^':  ST[ist-1] = pow(ST[ist-1],ST[ist]);break;
        }
    }
    return ST[0];
}
int main(int argc,char **argv) {
    if (argc<=1) {
        if (GetConsoleOutputCP()!=936) system("chcp 936>NUL");//中文代码页
        cout << "计算函数表达式的值。"<<endl<<"支持(),+,-,*,/,^,Pi,sin,cos,tan,sqrt,arcsin,arccos,arctan,lg,ln,exp"<<endl;
        while (1) {
            cout << "请输入表达式：";
            gets(s);
            if (s[0]==0) break;//
            cout << s <<"=";
            cout << setprecision(15) << calc(s,&endss) << endl;
        }
    } else {
        strncpy(s,argv[1],MAXLEN-1);s[MAXLEN-1]=0;
        if (argc>=3) {
            pcs=atoi(argv[2]);
            if (pcs<0||15<pcs) pcs=15;
            printf("%.*lf\n",pcs,calc(s,&endss));
        } else {
            printf("%.15lg\n",calc(s,&endss));
        }
    }
    return 0;
}

[benjaminweber]

下载后看了看代码，有如下的意见： 1. 你的parser是自己写的，语法被隐藏到代码中了。以后扩展语法会很不方便。建议用Bison/ANTLR重写。 2. 你的cmd也是自己写的，勇气可嘉，但是功能不好。建议使用tcl的shell环境，查找cppTcl库，并动态扩展数学功能。 3. 你的计算器不能读入文件，遇到大计算只能每次自己输入，很麻烦。建议支持读取用户代码。 4. 没有与其他软件的接口。比如可以将代码转成Matlab格式或支持直接读入Matlab代码。

[zhcosin]

引用 2 楼 benjaminweber 的回复:

下载后看了看代码，有如下的意见： 1. 你的parser是自己写的，语法被隐藏到代码中了。以后扩展语法会很不方便。建议用Bison/ANTLR重写。 2. 你的cmd也是自己写的，勇气可嘉，但是功能不好。建议使用tcl的shell环境，查找cppTcl库，并动态扩展数学功能。 3. 你的计算器不能读入文件，遇到大计算只能每次自己输入，很麻烦。建议支持读取用户代码。 ……

恩，谢谢建议。 1. 语法部分的确如你所说，语法被隐藏到代码中了，以后采用语法制导的语法分析。 2. cmd 那个命令循环，功能很弱，呵呵，没花功夫在这个上面，谢谢你提供的工具。 3. 恩，这个功能，可以添加，:-)

[赵4老师]

由 yu_hua 于2007-07-27设计完成 不是我写的，我只是拿来用。供参考而已。

[cbzjzsb123]

[rickys2080]

我也写过！！！

[zhcosin]

引用 1 楼 zhao4zhong1 的回复:

仅供参考： C/C++ code?1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586878889……

你给的这个小程序有点意思，不过如下： 1. 你那有两个小小的编译错误。 2. 你支持的那些函数怎么都是只带一个参数的？不能支持多个参数的函数么？

[c_rrb]

没有那边的帐号，楼主最好发表一份到CSDN来，谢谢啊

[zhcosin]

添加了对变量的支持。

[zhcosin]

以后想实现一种简单的编程语言，这样就可以定义有控制流程的函数了，从而可以完成更加复杂的计算任务，正在啃编译原理.....

[ppsharp]

给个CSDN的链接就好了

[zhcosin]

引用 11 楼 ppsharp 的回复:

给个CSDN的链接就好了

http://download.csdn.net/detail/zhcosin/5008980

[FancyMouse]

这程序，用来练手是一回事。真的要拿来用的话我还是用bc。

[zhcosin]

引用 13 楼 FancyMouse 的回复:

这程序，用来练手是一回事。真的要拿来用的话我还是用bc。

那是自然，这玩意目前还只是一个实验品，但我想等到我为它添加了可编程的功能，它的实用性当可以大大增加，当然了，还有一些非核心功能也是需要逐个完善的，比如对于读取配置文件，支持命令行的编辑和历史记录等。