// zero. This detects things like eps.i.j.k * p.j * p.k = 0.
if (local_dummy_indices.size() >= 2) {
lst dummy_syms;
- for (int i=0; i<local_dummy_indices.size(); i++)
+ for (exvector::size_type i=0; i<local_dummy_indices.size(); i++)
dummy_syms.append(local_dummy_indices[i].op(0));
if (r.symmetrize(dummy_syms).is_zero()) {
free_indices.clear();
// Yes, construct list of all dummy index symbols
lst dummy_syms;
- for (int i=0; i<dummy_indices.size(); i++)
+ for (exvector::size_type i=0; i<dummy_indices.size(); i++)
dummy_syms.append(dummy_indices[i].op(0));
// Chop the sum into terms and symmetrize each one over the dummy
// (a/b)^-x -> {sym((b/a)^x), 1}
return (new lst(replace_with_symbol(power(n_basis.op(1) / n_basis.op(0), -n_exponent), sym_lst, repl_lst), _ex1))->setflag(status_flags::dynallocated);
}
-
- } else { // n_exponent not numeric
-
- // (a/b)^x -> {sym((a/b)^x, 1}
- return (new lst(replace_with_symbol(power(n_basis.op(0) / n_basis.op(1), n_exponent), sym_lst, repl_lst), _ex1))->setflag(status_flags::dynallocated);
}
}
+
+ // (a/b)^x -> {sym((a/b)^x, 1}
+ return (new lst(replace_with_symbol(power(n_basis.op(0) / n_basis.op(1), n_exponent), sym_lst, repl_lst), _ex1))->setflag(status_flags::dynallocated);
}
// Contraction of two epsilon tensors is a determinant
bool variance = is_a<varidx>(self->op(1));
matrix M(num, num);
- for (int i=0; i<num; i++) {
- for (int j=0; j<num; j++) {
+ for (unsigned i=0; i<num; i++) {
+ for (unsigned j=0; j<num; j++) {
if (minkowski)
M(i, j) = lorentz_g(self->op(i+1), other->op(j+1), pos_sig);
else if (variance)