Searched refs:energy_1 (Results 1 – 1 of 1) sorted by relevance
16 PetscReal energy_0 = 0, energy_1 = 0; in main() local82 …for (PetscInt p = 0; p < Np; ++p) energy_1 += wq[p] * (PetscSqr(coords[p * 2 + 0]) + PetscSqr(coor… in main()85 …SELF, "Energy = %20.12e error = %20.12e\n", (double)energy_0, (double)((energy_1 - energy_0) / ene… in main()