
proc goblinImport {object fileName format} {

    global goblinRandGeometry goblinRandLength
    

    # Native GOBLIN Format
	
	if {[string equal $format "gob"]} {
        goblin read $object $fileName
		return
    }
	
	
	# DIMACS Formats

    if {[string equal $format "edge"]} {
		set file [open "$fileName" r]
		set n 0

		while {[gets $file thisLine] >= 0} {
		    if {[scan $thisLine "p edge %d %d" n m]==2} {
        		if {$n<=0 || $m<=0} {
                    close $file
		            error "Insufficient problem dimensions"
		        }
                set goblinRandGeometry 0
                set goblinRandLength 0
        		goblin sparse graph $object $n
                $object set ucap 1
                $object set lcap 0
                $object set length 1
                $object set demand 1
    		}

            set length 1
		    if {[scan $thisLine "e %d %d %f" u v length]>=2} {
        		if {$u>$n || $u<1 || $v>$n || $v<1} {
                    close $file
					$object delete
		            error "Node index exeeds problem dimension"
		        }
        
        		set a [$object arc insert [expr $u-1] [expr $v-1]]
		        $object arc $a set length $length
		    }
		}

		close $file

        if {$n==0} {error "Missing problem line"}

		return
	}

    if {[string equal $format "max"]} {
		global goblinSource goblinTarget

		set file [open "$fileName" r]
		set n 0
		set goblinSource *
		set goblinTarget *

		while {[gets $file thisLine] >= 0} {
		    if {[scan $thisLine "p max %d %d" n m]==2} {
        		if {$n<=0 || $m<=0} {
                    close $file
		            error "Insufficient problem dimensions"
		        }
                set goblinRandGeometry 0
        		goblin sparse digraph $object $n
        		$object set demand 0
    		}
    
		    if {[scan $thisLine "n %d %s" u type]==2} {
        		if {$u>$n || $u<1} {
                    close $file
					$object delete
		            error "Node index exeeds problem dimension"
		        }

        		if {[string equal $type "s"]} {set goblinSource [expr $u-1]}
		        if {[string equal $type "t"]} {set goblinTarget [expr $u-1]}
    		}
    
		    if {[scan $thisLine "a %d %d %f" u v cap]==3} {
        		if {$u>$n || $u<1 || $v>$n || $v<1} {
                    close $file
					$object delete
		            error "Node index exeeds problem dimension"
		        }
        
        		set a [$object arc insert [expr $u-1] [expr $v-1]]
		        $object arc $a set ucap $cap
		    }
		}

		close $file

        if {$n==0} {error "Missing problem line"}

		if {[string equal $goblinSource "*"] || [string equal $goblinTarget "*"]} {
			$object delete
    		error "Missing source and/or target node"
		}

		return
    }

    if {[string equal $format "min"]} {
		set file [open "$fileName" r]
		set n 0

		while {[gets $file thisLine] >= 0} {
		    if {[scan $thisLine "p min %d %d" n m]==2} {
        		if {$n<=0 || $m<=0} {
                    close $file
		            error "Insufficient problem dimensions"
		        }
                set goblinRandGeometry 0
        		goblin sparse digraph $object $n
        		$object set demand 0
    		}
    
		    if {[scan $thisLine "n %d %f" u demand]==2} {
        		if {$u>$n || $u<1} {
                    close $file
					$object delete
		            error "Node index exeeds problem dimension"
		        }
        
        		$object node [expr $u-1] set demand [expr -$demand]
    		}
    
		    if {[scan $thisLine "a %d %d %f %f %f" u v lcap ucap length]==5} {
        		if {$u>$n || $u<1 || $v>$n || $v<1} {
                    close $file
					$object delete
		            error "Node index exeeds problem dimension"
		        }
        
        		set a [$object arc insert [expr $u-1] [expr $v-1]]
		        $object arc $a set ucap $ucap
		        $object arc $a set lcap $lcap
		        $object arc $a set length $length
		    }
		}

		close $file

        if {$n==0} {error "Missing problem line"}

		return
    }

    if {[string equal $format "asn"]} {
		set file [open "$fileName" r]
		set n 0

		while {[gets $file thisLine] >= 0} {
		    if {[scan $thisLine "p asn %d %d" n m]==2} {
        		if {$n<=0 || $m<=0} {
                    close $file
		            error "Insufficient problem dimensions"
		        }
                set goblinRandGeometry 0
        		goblin sparse digraph $object $n
        		$object set demand 1
    		}
    
		    if {[scan $thisLine "n %d %s" u type]==2} {
        		if {$u>$n || $u<1} {
                    close $file
					$object delete
		            error "Node index exeeds problem dimension"
		        }
        
        		if {[string equal $type "s"]} {
                    $object node [expr $u-1] set demand -1
                } else {
                    $object node [expr $u-1] set demand 1
                }
    		}
    
		    if {[scan $thisLine "a %d %d %f" u v length]==3} {
        		if {$u>$n || $u<1 || $v>$n || $v<1} {
                    close $file
					$object delete
		            error "Node index exeeds problem dimension"
		        }
        
        		set a [$object arc insert [expr $u-1] [expr $v-1]]
		        $object arc $a set length $length
		    }
		}

		close $file

        if {$n==0} {error "Missing problem line"}

		return
	}

    if {[string equal $format "geom"]} {
		set file [open "$fileName" r]
		set n 0

		while {[gets $file thisLine] >= 0} {
		    if {[scan $thisLine "p geom %d %d" n dim]==2} {
        		if {$n<=0} {
                    close $file
		            error "Insufficient problem dimension"
		        }
        		if {$dim!=2} {
                    close $file
		            error "GOBLIN allows only 2-dim. representations"
		        }
				goblin dense graph $object $n
				$object set demand 1
				
				global goblinMethGeometry
				set goblinMethGeometry 2
    		}
    
		    if {[scan $thisLine "n %d %f %f" u x y]==3} {
        		if {$u>$n || $u<1} {
                    close $file
					$object delete
		            error "Node index exeeds problem dimension"
		        }
        
        		$object node [expr $u-1] set cx $x
        		$object node [expr $u-1] set cy $y
    		}
		}

		close $file

        if {$n==0} {error "Missing problem line"}

		return
    }


	# TSPLIB Formats

    if {[string equal $format "tsp"]} {
		set file [open "$fileName" r]

        set display_data ""
        while {[gets $file thisLine]>=0 && \
                [scan $thisLine {%[A-_] : %s} spec value]==2} {
            if {[string equal $spec "TYPE"]} {set type $value}
            if {[string equal $spec "DIMENSION"]} {set n $value}
            if {[string equal $spec "EDGE_WEIGHT_TYPE"]} {set geometry $value}
            if {[string equal $spec "EDGE_WEIGHT_FORMAT"]} {set format $value}
            if {[string equal $spec "DISPLAY_DATA_TYPE"]} {set display_data $value}
        }

        if {![info exists type]} {
            close $file
            error "Missing TYPE specifier"
        }
        if {![string equal $type "TSP"]} {
            close $file
            error "TYPE must be TSP"
        }
        if {![info exists n]} {
            close $file
            error "Missing DIMENSION specifier"
        }
        if {![info exists geometry]} {
            close $file
            error "Missing EDGE_WEIGHT_TYPE specifier"
        }

        set goblinRandGeometry 0
        goblin dense graph $object $n
        $object set demand 2
        
        global goblinMethGeometry
        if {[string equal $geometry "EXPLICIT"]} {
            set goblinMethGeometry 0

            if {![info exists format]} {
                close $file
                $object delete
                error "Missing EDGE_WEIGHT_FORMAT specifier"
            }

            if {![string equal [string trim $thisLine] "EDGE_WEIGHT_SECTION"]} {
                close $file
                $object delete
                error "Missing EDGE_WEIGHT_SECTION specifier"
            }
            
            set goblinMethGeometry 0
            gets $file thisLine
            set thisLine [string trim $thisLine]
            
            if {[string equal $format "LOWER_DIAG_ROW"]} {
                set iLower {0}
                set iUpper {$n}
                set jLower {0}
                set jUpper {$i+1}
            } else {
                if {[string equal $format "UPPER_DIAG_ROW"]} {
                    set iLower {0}
                    set iUpper {$n}
                    set jLower {$i}
                    set jUpper {$n}
                } else {
                    if {[string equal $format "UPPER_ROW"]} {
                        set iLower {0}
                        set iUpper {$n-1}
                        set jLower {$i+1}
                        set jUpper {$n}
                    } else {
                        if {[string equal $format "FULL_MATRIX"]} {
                            set iLower {0}
                            set iUpper {$n}
                            set jLower {0}
                            set jUpper {$n}
                        } else {
                            close $file
                            $object delete
                            error "EDGE_WEIGHT_FORMAT: $format not supported in GOBLIN"
                        }
                    }
                }
            }
            
            for {set i [expr $iLower]} {$i<[expr $iUpper]} {incr i} {
                for {set j [expr $jLower]} {$j<[expr $jUpper]} {incr j} {
                    if {[scan $thisLine {%[0-9]} thisWeight]==0} {
                        close $file
                        $object delete
                        error "Parse error in EDGE_WEIGHT_SECTION"
                    }
                        
                    $object arc [$object adjacency $i $j] set length $thisWeight
                        
                    set firstBlank [string first " " $thisLine]
                    incr firstBlank
                    if {$firstBlank>0} {
                        set thisLine [string range $thisLine $firstBlank end]
                        set thisLine [string trim $thisLine]
                    } else {
                        gets $file thisLine
                        set thisLine [string trim $thisLine]
                    }
                }
            }
            
            for {set i 0} {$i<$n} {incr i} {
                $object arc [$object adjacency $i $i] set length *
            }

            if {[string equal $display_data "TWOD_DISPLAY"]} {
                if {![string equal [string trim $thisLine] "DISPLAY_DATA_SECTION"]} {
                    close $file
                    $object delete
                    error "Missing DISPLAY_DATA_SECTION specifier"
                }
				
                while {[gets $file thisLine] >= 0 && \
                        ![string equal [string trim $thisLine] "EOF"]} {
		            if {[scan $thisLine "%d %f %f" u x y]==3} {
            		    $object node [expr $u-1] set cx $x
            		    $object node [expr $u-1] set cy $y
    	    	    }
		        }
            }
        } else {
            if {[string equal $geometry "EUC_2D"]} {
                set goblinMethGeometry 2
            } else {
                if {[string equal $geometry "MAX_2D"]} {
                    set goblinMethGeometry 3
                } else {
                    if {[string equal $geometry "MAN_2D"]} {
                        set goblinMethGeometry 1
                    } else {
                        if {[string equal $geometry "GEO"]} {
                            set goblinMethGeometry 4
                        } else {
                            close $file
                            $object delete
                            error "Format EDGE_WEIGHT_TYPE: $geometry is not supported"
                        }
                    }
                }
            }

            if {![string equal [string trim $thisLine] "NODE_COORD_SECTION"]} {
                close $file
                $object delete
                error "Missing NODE_COORD_SECTION specifier"
            }
				
            while {[gets $file thisLine] >= 0 && \
                    ![string equal [string trim $thisLine] "EOF"]} {
		        if {[scan $thisLine "%d %f %f" u x y]==3} {
            		$object node [expr $u-1] set cx $x
            		$object node [expr $u-1] set cy $y
    	    	}
		    }
        }

		close $file
		return
    }
	
    if {[string equal $format "atsp"]} {
		set file [open "$fileName" r]

        set display_data ""
        while {[gets $file thisLine]>=0 && \
                [scan $thisLine {%[A-_] : %s} spec value]==2} {
            if {[string equal $spec "TYPE"]} {set type $value}
            if {[string equal $spec "DIMENSION"]} {set n $value}
            if {[string equal $spec "EDGE_WEIGHT_TYPE"]} {set geometry $value}
            if {[string equal $spec "EDGE_WEIGHT_FORMAT"]} {set format $value}
            if {[string equal $spec "DISPLAY_DATA_TYPE"]} {set display_data $value}
        }

        if {![info exists type]} {
            close $file
            error "Missing TYPE specifier"
        }
        if {![string equal $type "ATSP"]} {
            close $file
            error "TYPE must be ATSP"
        }
        if {![info exists n]} {
            close $file
            error "Missing DIMENSION specifier"
        }
        if {![info exists geometry]} {
            close $file
            error "Missing EDGE_WEIGHT_TYPE specifier"
        }

        set goblinRandGeometry 0
        goblin dense digraph $object $n
        $object set demand 0
        
        global goblinMethGeometry
        if {[string equal $geometry "EXPLICIT"]} {
            set goblinMethGeometry 0

            if {![info exists format]} {
                close $file
                $object delete
                error "Missing EDGE_WEIGHT_FORMAT specifier"
            }

            if {![string equal [string trim $thisLine] "EDGE_WEIGHT_SECTION"]} {
                close $file
                $object delete
                error "Missing EDGE_WEIGHT_SECTION specifier"
            }
            
            set goblinMethGeometry 0
            gets $file thisLine
            set thisLine [string trim $thisLine]
            
            if {[string equal $format "FULL_MATRIX"]} {
                set iLower {0}
                set iUpper {$n}
                set jLower {0}
                set jUpper {$n}
            } else {
                close $file
                $object delete
                error "EDGE_WEIGHT_FORMAT: $format not supported in GOBLIN"
            }
            
            for {set i [expr $iLower]} {$i<[expr $iUpper]} {incr i} {
                for {set j [expr $jLower]} {$j<[expr $jUpper]} {incr j} {
                    if {[scan $thisLine {%[0-9]} thisWeight]==0} {
                        close $file
                        $object delete
                        error "Parse error in EDGE_WEIGHT_SECTION"
                    }
                        
                    $object arc [$object adjacency $i $j] set length $thisWeight
                        
                    set firstBlank [string first " " $thisLine]
                    incr firstBlank
                    if {$firstBlank>0} {
                        set thisLine [string range $thisLine $firstBlank end]
                        set thisLine [string trim $thisLine]
                    } else {
                        gets $file thisLine
                        set thisLine [string trim $thisLine]
                    }
                }
            }
            
            for {set i 0} {$i<$n} {incr i} {
                $object arc [$object adjacency $i $i] set length *
            }

            if {[string equal $display_data "TWOD_DISPLAY"]} {
                if {![string equal [string trim $thisLine] "DISPLAY_DATA_SECTION"]} {
                    close $file
                    $object delete
                    error "Missing DISPLAY_DATA_SECTION specifier"
                }
				
                while {[gets $file thisLine] >= 0 && \
                        ![string equal [string trim $thisLine] "EOF"]} {
		            if {[scan $thisLine "%d %f %f" u x y]==3} {
            		    $object node [expr $u-1] set cx $x
            		    $object node [expr $u-1] set cy $y
    	    	    }
		        }
            }
        } else {
            close $file
            $object delete
            error "Format EDGE_WEIGHT_TYPE: $geometry is not supported"
        }

		close $file
		return
    }
	
    if {[string equal $format "hcp"]} {
		set file [open "$fileName" r]

        while {[gets $file thisLine]>=0 && \
                [scan $thisLine {%[A-_] : %s} spec value]==2} {
            if {[string equal $spec "TYPE"]} {set type $value}
            if {[string equal $spec "DIMENSION"]} {set n $value}
            if {[string equal $spec "EDGE_DATA_FORMAT"]} {set format $value}
        }

        if {![info exists type]} {
            close $file
            error "Missing TYPE specifier"
        }
        if {![string equal $type "HCP"]} {
            close $file
            error "TYPE must be HCP"
        }
        if {![info exists n]} {
            close $file
            error "Missing DIMENSION specifier"
        }
        if {![info exists format]} {
            close $file
            error "Missing EDGE_DATA_FORMAT specifier"
        }

        if {![string equal [string trim $thisLine] "EDGE_DATA_SECTION"]} {
            close $file
            error "Missing EDGE_DATA_SECTION specifier"
        }

        set goblinRandGeometry 0
        goblin sparse graph $object $n
        $object set demand 2
        
        global goblinMethGeometry
        set goblinMethGeometry 0
            
        if {[string equal $format "EDGE_LIST"]} {
            while {[gets $file thisLine]>=0 && \
                [scan $thisLine {%d %d} head tail]==2} {
                $object arc insert [expr $head-1] [expr $tail-1]
            }
        } else {
            close $file
            $object delete
            error "Format EDGE_DATA_TYPE: $format is not supported"
        }

		close $file
		return
    }
	
	error "Unknown Format: $format"
}
                    
