Home My Page Projects Code Snippets Project Openings diderot
Summary Activity Tracker Tasks SCM

SCM Repository

[diderot] Diff of /branches/charisee/src/compiler/ein/mkoperators.sml
ViewVC logotype

Diff of /branches/charisee/src/compiler/ein/mkoperators.sml

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 2842, Tue Dec 2 01:38:14 2014 UTC revision 2843, Mon Dec 8 01:27:25 2014 UTC
# Line 153  Line 153 
153              E.Prod[ E.Epsilon(0, 1, 2), E.Tensor(0, [E.V 1]),  E.Tensor(1, [E.V 2 ]) ])              E.Prod[ E.Epsilon(0, 1, 2), E.Tensor(0, [E.V 1]),  E.Tensor(1, [E.V 2 ]) ])
154      }      }
155    
156        (*2-d cross product Eps_{ij}U_i V_j*)
157        val crossProduct2 = E.EIN{
158            params = [E.TEN(1,[2]), E.TEN(1,[2])],
159            index= [],
160            body=E.Sum([(E. V 0,0,1),(E.V 1,0,1)],
161            E.Prod[ E.Eps2(0, 1), E.Tensor(0, [E.V 0]),  E.Tensor(1, [E.V 1]) ])
162            }
163    
164    
165      (* Identiy: <\delta_{i j}>_{i j}  *)      (* Identiy: <\delta_{i j}>_{i j}  *)
166      fun identity(dim) =E.EIN{      fun identity(dim) =E.EIN{
167          params = [],          params = [],
# Line 298  Line 307 
307          end          end
308    
309      (*(F_y/dx - F_x/dy )k*)      (*(F_y/dx - F_x/dy )k*)
310      val curl2d=E.EIN{      (*val curl2d=E.EIN{
311          params = [E.FLD 2],          params = [E.FLD 2],
312          index = [],          index = [],
313          body = E.Sub(E.Apply(E.Partial([E.C 0]), E.Field(0,[E.C 1])),          body = E.Sub(E.Apply(E.Partial([E.C 0]), E.Field(0,[E.C 1])),
314                       E.Apply(E.Partial([E.C 1]), E.Field(0,[E.C 0])))                       E.Apply(E.Partial([E.C 1]), E.Field(0,[E.C 0])))
315        }*)
316    
317        (*\EinExp{\sum_{ij}\mathcal{E}_{ij} \frac{ F_j}{\partial x_i}*)
318            val curl2d=E.EIN{
319            params = [E.FLD 2],
320            index = [],
321            body = E.Sum([(E.V 0,0,1), (E.V 1,0,1)],
322                    E.Prod[E.Eps2(0, 1),
323                    E.Apply( E.Partial([E.V 0]), E.Field(0,[E.V 1]))])
324    
325      }      }
326    
327    
328      val curl3d=E.EIN{      val curl3d=E.EIN{
329          params = [E.TEN(1,[3])],          params = [E.TEN(1,[3])],
330          index = [3],          index = [3],

Legend:
Removed from v.2842  
changed lines
  Added in v.2843

root@smlnj-gforge.cs.uchicago.edu
ViewVC Help
Powered by ViewVC 1.0.0