Dear Matt and Meep users,
I have read your(matt) post in meep discuss about anisotropy material(
http://www.mail-archive.com/[email protected]/msg01405.html ),
which have resolved the setting of epsilon-tensor very well.
However when I used your method to set the D-conductivity tensor, there are
some problems.
My code of definition as follows:
*(define-class dielectric material-type
(define-property epsilon no-default 'number)
(define-property epsilon-tensor (matrix3x3 (vector3 1 0 0)
(vector3 0 1 0)
(vector3 0 0 1)) 'matrix3x3)
(define-derived-property scaled-epsilon-tensor 'matrix3x3
(lambda (object)
(matrix3x3* (object-property-value
object 'epsilon)
(object-property-value object 'epsilon-tensor))))
(define-property D-conductivity no-default 'number)
(define-property D-conductivity-tensor (matrix3x3 (vector3 1
0 0) (vector3 0 1 0)
(vector3 0 0 1)) 'matrix3x3)
(define-derived-property scaled-D-conductivity-tensor
'matrix3x3
(lambda (object)
(matrix3x3* (object-property-value object
'D-conductivity)
(object-property-value object 'D-conductivity-tensor)))))
(define MO1 (make dielectric (epsilon 1) (epsilon-tensor (matrix3x3
(vector3 4.8 0 0) (vector3 0 4.8 0) (vector3 0 0 5.2)))
(D-conductivity 0)
(D-conductivity-tensor (matrix3x3 (vector3 4.8 0 0) (vector3 0 4.8 0)
(vector3 0 0 5.2)))))*
The output of Meep is :
meep> MO1
*((dielectric material-type) (scaled-epsilon-tensor . #(#(4.8 0 0) #(0 4.8
0) #(0 0 5.2))) (epsilon . 1) (epsilon-tensor . #(#(4.8 0 0) #(0 4.8 0) #(0
0 5.2))) (D-conductivity . 0) (D-conductivity-tensor . #(#(4.8 0 0) #(0 4.8
0) #(0 0 5.2))) (scaled-D-conductivity-tensor . #(#(0 0 0) #(0 0 0) #(0 0
0))))
1. From the above, although no error happens, the **D-conductivity-tensorand
**scaled-D-conductivity-tensor are not equal, which is a problem to be
solved.
2 When I started a simulation using the above definition of MO1, the
output of meep :
-----------
Initializing structure...
Working in 2D dimensions.
meep: unknown material type
cylinder, center = (0,0,0)
radius 4, height 1e+20, axis (0, 0, 1)
cylinder, center = (0,0,0)
radius 3, height 1e+20, axis (0, 0, 1)
cylinder, center = (0,0,0)
radius 2, height 1e+20, axis (0, 0, 1)
cylinder, center = (0,0,0)
radius 1.2, height 1e+20, axis (0, 0, 1)
dielectric constant epsilon = 1
cylinder, center = (0,0,0)
radius 1, height 1e+20, axis (0, 0, 1)
Could you and someone else gives me some suggestions, thanks ahead!
Huang
**
*
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
