while (op_top >= 0) { char op = op_stack[op_top]; op_top--; double b = num_stack[num_top]; num_top--; double a = num_stack[num_top]; num_top--; double result = 0; if (op == '+') { result = a + b; } else if (op == '-') { result = a - b; } else if (op == '*') { result = a * b; } else if (op == '/') { result = a / b; } else if (op == '%') { result = (int)a % (int)b; } else if (op == 's') { result = sqrt(a); } else if (op == '^') { result = pow(a, b); } num_top++; num_stack[num_top] = result; } return num_stack[0]; }详细解释一下这个程序

这是一个计算器程序的核心部分，能够实现对表达式的计算。程序使用两个栈来实现计算，一个是操作数栈，一个是操作符栈。先将表达式中的数字和操作符分别压入栈中，然后依次取出操作符和对应的操作数进行计算，再将计算结果压入操作数栈中，直到操作符栈为空，最后将操作数栈中的结果返回。

具体来说，程序的主要流程如下：

1. 初始化两个栈：操作数栈和操作符栈
2. 从左往右扫描表达式中的每个字符，如果是数字则压入操作数栈中，如果是操作符则判断优先级，如果当前操作符的优先级小于或等于操作符栈顶元素的优先级，则从操作数栈中取出两个操作数，从操作符栈中取出一个操作符进行计算，将计算结果压入操作数栈中，重复此过程直到当前操作符的优先级大于操作符栈顶元素的优先级，最后将当前操作符压入操作符栈中
3. 如果表达式扫描完毕，但操作符栈中还有元素，则依次取出操作符和对应的操作数进行计算，直到操作符栈为空
4. 返回操作数栈中的唯一元素，即为计算结果

在计算时，程序支持加减乘除取模、开方和幂运算。其中加减乘除取模的计算方式比较简单，直接进行相应的运算即可，而开方和幂运算需要调用相关的数学函数来实现。