Adding point integral class
A suggested interface for the point_integral feature can be found at:
lp:~martinal/ufc/pending3_point_integral_addition
Any suggestions?
Blueprint information
- Status:
- Started
- Approver:
- None
- Priority:
- Essential
- Drafter:
- None
- Direction:
- Needs approval
- Assignee:
- None
- Definition:
- Drafting
- Series goal:
- None
- Implementation:
-
Beta Available
- Milestone target:
- None
- Started by
- Marie Rognes
- Completed by
Whiteboard
MER: The suggested interface looks good -- I've reimplemented this (easier b/c of flurry of ufc changes since original branch). In addition to the generic point, I suggest we add a function for local vertices as below. We can generate significantly more efficient code for that, and it is a common use case.
/// Tabulate the tensor for the contribution from a local vertex
virtual void point_integral:: tabulate_
MER: I've added the above, with tabulate_tensor both for a set of arbitrary points (relative to a cell) and for a local vertex index now.