GNU bug report logs - #23553
Error in running some loop in set! geometry

Previous Next

To add a comment to this bug, you must first unarchive it, by sending
a message to control AT debbugs.gnu.org, with unarchive 23553 in the body.
You can then email your comments to 23553 AT debbugs.gnu.org in the normal way.

Toggle the display of automated, internal messages from the tracker.

Message #5 received at submit <at> debbugs.gnu.org (full text, mbox):

From: Partha Saha <psaha089 <at> gmail.com>
To: bug-guile <at> gnu.org
Subject: Error in running some loop in set! geometry
Date: Fri, 13 May 2016 10:34:28 +0000 (UTC)

Dear Prof. Steve and all the meep users

I am new to meep and Scheme coding and I think I have a problem making a
loop in Scheme but I don't know what I have done wrong.

I am trying to make a 1D slotted nanobeam slab but finding some error in
running the loop. My ctl file is as follows:

; Some parameters to describe the geometry:                                     
(define-param eps 12) ; dielectric constant of waveguide                        
(define-param w 1.1627906977) ; width of waveguide                         
             
(define-param r_center 0.4203935599) ; radius of holes at the center
(define-param r_end 0.3631484794) ; radius of holes at the end
(define-param h 0.3935599284) ; height of slab  
(define-param w_slot 0.0894454383) ; width of slot                         
            
;(define-param d 1.4) ; defect spacing (ordinary spacing = 1) but in my case
no defect hence disabled                   
(define-param N 30) ; number of holes on either side of defect      
(define-param i-max 20) ; number of holes on either side of tapered region  
(define-param di 1); change in i   
(define-param dmirr 1); change in mirr  
(define-param mirr-max 10) ; number of holes on either side of mirror holes
region 
          

; The cell dimensions                                                           
(define-param sy 6) ; size of cell in y direction (perpendicular to wvg.)    
(define-param sz 3) ; size of cell in z direction (perpendicular to wvg.)     
(define-param pad 2) ; padding between last hole and PML edge                   
(define-param dpml 1) ; PML thickness 

(define sx (* 2 (+ pad dpml N))) ; size of cell in x direction



(set! geometry-lattice (make lattice (size sx sy sz)))



(set! geometry   
(append geometry; combine lists of objects  
   (list (make block (center 0 0 0) (size infinity w h)
               (material (make dielectric (epsilon eps))))

          (make block (center 0 0 0) (size (+ (* 2 N) 2) w_slot infinity)
(material air)) ; air slot in the waveguide

          (do ((i 1 (+ i di))) ((>= i i-max))  (make cylinder (center (+ 0.5
(* (+ i -1) 1)) 0 0) (radius (+ r_center (/ (* (* (+ i -1) (+ i -1)) (-
r_end r_center)) (* (+ i-max -1) (+ i-max -1))))) (height infinity)
(material air)))

          (do ((i 1 (+ i di))) ((>= i i-max))  (make cylinder (center (+
-0.5 (* (+ i -1) -1)) 0 0) (radius (+ r_center (/ (* (* (+ i -1) (+ i -1))
(- r_end r_center)) (* (+ i-max -1) (+ i-max -1))))) (height infinity)
(material air)))

          (do ((mirr 1 (+ mirr dmirr))) ((>= mirr mirr-max)) (make cylinder
(center (+ i-max mirr -0.5) 0 0) (radius r_end) (height infinity) (material
air)))

          (do ((mirr 1 (+ mirr dmirr))) ((>= mirr mirr-max)) (make cylinder
(center (- 0.5 i-max mirr) 0 0) (radius r_end) (height infinity) (material
air)))
        

);end of list
);end of append
);end of geometry




(set! pml-layers (list (make pml (thickness dpml))))

(set-param! resolution 20)


(define-param fcen 0.3653594) ; pulse center frequency
(define-param df 0.2)  ; pulse width (in frequency)
(define-param nfreq 500) ; number of frequencies at which to compute flux


(set! sources (list
               (make source
                 (src (make gaussian-src (frequency fcen) (fwidth df)))
                 (component Ey)
                 (center (+ dpml (* -0.5 sx)) 0 0)
                 (size 0 w 0))))

(set! symmetries (list (make mirror-sym (direction Y) (phase -1))))


(run-sources+ 300
              (at-beginning output-epsilon)
              (after-sources (harminv Ez (vector3 (+ N 0.2) 0 0) fcen df)))



The error message is
ERROR: In procedure cdr:
ERROR: In procedure cdr: Wrong type (expecting pair): #<unspecified>

The full error is
Backtrace:
In ice-9/boot-9.scm:
 157: 11 [catch #t #<catch-closure 249aa60> ...]
In unknown file:
   ?: 10 [apply-smob/1 #<catch-closure 249aa60>]
In ice-9/eval.scm:
 432: 9 [eval # #]
 432: 8 [eval # #]
In unknown file:
   ?: 7 [primitive-load "nanobeam_edited_with_do_loop.ctl"]
In ice-9/eval.scm:
 387: 6 [eval # ()]
 432: 5 [eval # #]
 432: 4 [eval # #]
 467: 3 [eval # #]
In unknown file:
   ?: 2 [make-structure 3 #(66.0 6.0 3.0) #(0.0 0.0 0.0) 20 #t 1.0e-4 100000
#f ...]
In ice-9/eval.scm:
 387: 1 [eval # #]
In unknown file:
   ?: 0 [cdr #<unspecified>]

ERROR: In procedure cdr:
ERROR: In procedure cdr: Wrong type (expecting pair): #<unspecified>

Can anyone help me in finding the error?

Message #10 received at 23553-done <at> debbugs.gnu.org (full text, mbox):

From: Andy Wingo <wingo <at> pobox.com>
To: Partha Saha <psaha089 <at> gmail.com>
Cc: 23553-done <at> debbugs.gnu.org
Subject: Re: bug#23553: Error in running some loop in set! geometry
Date: Tue, 28 Feb 2017 16:25:30 +0100

On Fri 13 May 2016 12:34, Partha Saha <psaha089 <at> gmail.com> writes:

> Dear Prof. Steve and all the meep users

This address is for bugs in Guile itself, not meep.  Closing this bug as
out of scope.  Happy hacking with Guile, though :)

Andy

Previous Next

GNU bug tracking system
Copyright (C) 1999 Darren O. Benham, 1997,2003 nCipher Corporation Ltd, 1994-97 Ian Jackson.